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FOREWORD 


This  report  describes  the  theoretical  basis  behind  the  low  frequency  propagation  and 
scattering  models  used  in  Personal  Computer  Shallow  Water  Acoustics  Tool-set 
(PC  SWAT  7.0).  PC  SWAT  is  a  user-friendly  sonar  simulation  developed  by  Dr.  Sammelmann. 
It  is  used  widely  throughout  the  Department  of  Defense. 
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1.0  INTRODUCTION 


This  report  describes  the  theoretical  basis  for  the  low  frequency  propagation  and 
scattering  models  used  by  PC  SWAT  7.0.  These  models  are  based  on  the  normal  mode 
description  of  propagation  of  sound  in  a  range  independent  waveguide. 

In  addition  to  the  Introduction  and  Reference  sections,  this  report  is  composed  of 
thirty  sections.  Section  2  gives  a  brief  derivation  of  the  acoustic  wave  equation.  Section  3 
provides  a  brief  review  of  the  theory  of  Sturm-Louville  operators.  Sections  4, 5,  and  6  provide  a 
description  of  the  construction  of  the  Greens’  Function  in  a  horizontally  stratified  waveguide. 
Section  7  provides  a  derivation  of  the  normalization  of  the  depth  functions  in  a  waveguide  in 
terms  of  the  boundary  conditions  at  the  top  and  bottom  of  the  waveguide.  Section  8  describes  the 
construction  of  the  normal  mode  description  of  propagation  in  a  waveguide  with  a  piece-wise 
constant  sound  speed  and  density  profile  and  a  rigid  basement.  Section  9  describes  the  extension 
of  the  methods  of  Section  8  to  include  a  homogeneous  half-space  as  the  basement  of  the 
waveguide.  Section  10  extends  the  normal  mode  description  to  include  piece- wise  linear  sound 
speed  profiles.  Sections  1 1  and  12  provide  examples  of  the  normal  mode  description  in  an 
isovelocity  waveguide  with  a  rigid  and  homogeneous  half-space  basement.  Sections  13  and  14 
describe  the  effects  of  small-scale  roughness  of  the  coherent  (mean)  field.  Section  15  provides 
the  derivation  of  the  time  domain  description  of  the  propagation  of  normal  modes  based  on  a 
saddle  point  analysis  of  the  propagation  of  a  normal  mode.  Sections  16, 17,  and  1 8  describe 
representations  of  cylindrical  and  spherical  wave  functions,  and  the  transformations  between 
these  two  sets  of  solutions  of  the  ffee-field  Helmholtz  Equation.  Section  19  describes  the 
expansion  of  the  normal  mode  contributions  to  propagation  in  terms  of  a  spherical  basis  set  about 
the  source  and  receive  points.  Section  20  describes  the  use  of  the  multipole  expansion  of  the 
Greens’  Function  in  the  inclusion  of  the  effects  of  the  directivity  of  the  source  and  the  receiver  on 
the  acoustic  field.  Section  21  describes  the  effects  of  rotating  the  local  coordinate  system  on  the 
spherical  expansions  of  the  Greens’  Function.  Section  22  presents  a  review  of  the 
Helmholtz  Integral  Equation  in  preparation  for  the  introduction  of  the  spherical  T-matrix 
approach  to  scattering  from  a  target.  Sections  23  and  24  derive  the  spherical  T-matrix  approach 
for  scattering  from  rigid  and  elastic  solids.  Section  25  describes  the  scattering  from  an  elastic 
target  in  a  waveguide  based  on  the  preceding  section  on  the  spherical  T-matrix.  Sections  26, 27, 
and  28  describe  the  normal  mode  and  time  domain  representation  of  the  scattering  from  a  rough 
interface  and  volume  inhomogeneities.  Section  29  describes  the  plane  wave  approximation  of 
the  computation  of  the  signal-to-noise  ratio  using  the  normal  mode  method.  Section  30  describes 
the  global  matrix  method  of  solving  for  the  normal  modes,  and  Section  31  presents  some 
examples  computed  using  PC  SWAT  7.0. 
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2.0  ACOUSTIC  WAVE  EQUATION 


Conservation  of  mass  and  momentum  are  given  by  the  hydrodynamic  equations  of 
motion1,2 : 


•fp  +  V.(A?)  =  0 
ot 

(2.1) 

d  _  _ 

p — v  +  pv  •  Vv  =  — VP+  pF 
dt 

(2.2) 

where  p  is  the  density  of  the  fluid,  v  is  the  velocity  of  the  fluid,  P  is  the  pressure  of  the  fluid,  and 
F  is  the  external  force  per  unit  mass  (for  example,  gravity).  The  acoustic  wave  equation  is 
defined  by  linearizing  the  above  hydrodynamic  equations  by  making  the  following  change  of 
variables  in  the  above  equations. 

p  =  p0+Sp 

(2.3) 

V  =  v0  +  Sv 

(2.4) 

p=p0+p 

(2.5) 

The  linearized  equations  of  motion  are: 

•f- Sp + v  •  (4?v0) + v  •  (p0Sv)  =  o 

at 

— Sv  +  v0  •  V<£  +  Sv  •  Vv0  =  Vp  +  VP0 
dt  A  Po 

If  one  neglects  the  thermal  conductivity  of  the  fluid,  and  diffusion  of  the  components  of 
the  fluid,  the  acoustic  propagation  of  a  wave  can  be  treated  as  an  adiabatic  process,  where  the 
density  and  pressure  fluctuations  are  proportional: 

Sp  =  c~2p  (2.8) 

In  the  above  equation,  c  is  the  speed  of  sound  in  the  fluid  given  by  the  partial  derivative  of  the 
pressure  in  the  fluid  with  respect  to  density  at  constant  entropy: 

c2=A(<%=  0)  (2.9) 

dp 


(2.6) 

(2.7) 
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Assume  that  in  the  absence  of  the  wave  the  medium  is  at  rest.  In  addition,  neglect 
gravitational  effects  on  the  wave,  that  is,  the  gradient  of  the  background  pressure  field  in 
Equation  2.7.  In  this  case,  the  linearized  equations  of  motion  are: 


d<£  X7 


(2.10) 


dp 

dt 


p0c2V  •  dv 


(2.11) 


Eliminating  the  velocity  of  the  acoustic  wave  in  the  above  pair  of  first  order  differential 
equations,  one  arrives  at  the  following  second  order  differential  equation  in  terms  of  the  pressure 
of  the  acoustic  wave  for  a  stationary  medium: 


Vp+|(-Lr^)  =  0 

ot  p0c  dt 


(2.12) 


In  the  case  of  a  monochromatic  wave  with  em  time  dependence,  the  acoustic  wave 
equation  takes  on  the  following  form: 

P$m  —  Vp  +  k2p  =  0  (2.13) 

Po 


where 


c 

is  the  acoustic  wavenumber  of  the  wave. 


(2.14) 
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3.  STURM-LOUVILLE  OPERATORS 


This  section  provides  a  brief  overview  of  the  theory  of  Sturm-Louville  operators.  An 
overview  of  the  Sturm-Louville  Equation  can  be  found  in  Morse  and  Feshbach3,4. 

The  Sturm-Louville  Differential  Equation  is  a  linear,  second  order  differential  equation  of 
the  following  form: 

--^-p(x)~y/(x)  +  (q(x)  -  lw(x))y/(x)  =  0  (3.1) 

ax  ax 


The  differential  operator 


L  =  i-7^7  +  (q(x)  ~  Mx))} 

w(  x)  ax  ax 


(3.2) 


is  a  self-adjoint  operator  on  the  space  of  smooth  functions  C°°(a>b)  on  which  the  following  inner 
product  is  defined: 

h 

<  y/\  y/'>=  J^(jc)^'(*)w(*)<& 


THEOREM  1 

This  theorem  establishes  the  criterion  for  a  pair  of  solutions  of  the  Sturm-Louville 
Equation  to  be  linearly  independent. 

If  the  functions  p(x)  and  q(x)  are  continuous  function  in  the  interval  (a,b)  and  y/^x)  and  y/2(x) 
are  a  pair  of  solutions  of  the  differential  equation: 


^2  ^ 

-77  ¥(*)  +  P(x)—V(x)  +  q(x)y/(x)  =  0 

ax  ax 

for  which  the  Wronskian  (W)  is  defined  below: 

2)  =  yr,(x) ^ ~  V2{x) 

ax  ax 


(3.3) 


(3.4) 
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is  non-zero  for  a  point  in  the  interval  (a,b).  Then  every  solution  of  this  differential  equation  in 
this  interval  is  a  linear  combination  of  these  two  functions. 


THEOREM  2 

This  theorem  yields  information  about  the  Wronskian  of  a  pair  of  solutions  of  the 
Sturm-Louville  Equation. 

If  (x)  and  y/2(x)  are  solutions  of  the  Sturm-Louville  Differential  Equation: 
d  d 

-—/>(*)— IK*)  +  (<?(*)  -  hv(x))y/{x)  =  0  (3.5) 

Then  the  Wronskian  of  these  two  functions  is  of  the  following  form: 

fr(r„ri)  =  0.6) 

dx  dx 


THEOREM  3 

This  theorem  states  that  two  eigenfunctions  of  the  Sturm-Louville  Equation  are 
necessarily  orthogonal  if  they  have  different  eigenvalues.  Note,  the  converse  statement  is  not 
true. 

Let  \j/(x)  and  \f/'{x)  be  two  eigenfunctions  of  the  Sturm-Louville  Equation: 
d  d 

POO— IK*)  +  (q(x)  -  Xw(x))y/{x)  =  0  (3.7) 

with  eigenvalues  A  &  A' ,  respectively.  Then  the  inner  product 

h 

<  y/ 1  y/'>-  jy/T x)\y' (x)w(x)dx  =  0 

a 

vanishes. 
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THEOREM  4 

This  theorem  describes  the  orthonormality  and  completeness  of  the  spectrum  of  the 
Sturm-Louville  Equation.  Equations  3.9  and  3.10  describe  the  orthonormality  of  a  suitable  basis 
of  eigenfunctions  of  the  Sturm-Louville  Equation.  Equation  3.11  describes  the  completeness 
condition  for  this  basis. 

Suppose  the  eigenvectors  of  the  Sturm-Louville  Equation 

+  (<?U)  “  Xw(xj)y/(x)  =  0  (3.8) 

ax  ax 

consist  of  a  discrete  set  of  eigenvectors  y/m(x)  with  a  discrete  spectrum  and  a  continuous  set  of 
eigenvectors  ^(*)  with  a  continuous  spectrum.  Then  one  can  normalize  a  complete  set  of 
eigenvectors,  such  that  the  following  orthonormality  conditions  are  valid: 

<¥m\Vn  >=  ^ 

<¥x\  Wx>=8(A-X)  (3*10) 


and  the  following  completeness  condition  is  valid. 


IX VWJJ)  +  \¥xix)y/x{x')dX  =  ^S(x  ~  *’) 


(3.11) 


THEOREM  5 

This  theorem  describes  the  construction  of  the  Green’s  Function  for  the 
Sturm-Louville  Equation  from  a  complete  orthonormal  basis  of  the  spectrum  of  the  equation. 

Let  G(x,x':  X)  be  a  solution  of  the  inhomogeneous  Sturm-Louville  Equation: 


-  —  p(x) — G(x,  x':  A)  +  (q(x)  -  Aw(r))G(r ,  x' :  X)  —  -8 ( x  -  x' ) 
dx  dx 

(3.12) 

subject  to  the  boundary  conditions: 

axG{a,x':  X)  +  fa—G{a,x':  X)  =  0 
dx 

(3.13) 

a2G(b,x':  X)  +  J32—G(b,x':  X)  =  0 
dx 

(3.14) 
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Let  y/m(x)  be  a  complete,  orthonormal  basis  of  the  discrete  spectrum,  and  if/(x,A)  be  a 

complete,  orthonormal  basis  of  the  continuous  spectrum  as  described  in  Theorem  4.  Then  the 
solution  of  the  above  differential  equation  for  the  Green’s  Function  of  the 
Sturm-Louville  Equation  is  of  the  form: 


G(x,x':A)  = 


YVMVmW)  ,  1 
^  (A-Am)  2m 


J  (A  -A1) 


(3.15) 


in  terms  of  this  basis. 


THEOREM  6: 

This  theorem  describes  an  alternative  solution  for  the  Green’s  Function  of  the 
Sturm-Louville  Equation,  which  does  not  require  the  construction  of  an  orthonormal  basis  of  the 
homogeneous  equation. 

Let  G(x,x':  A)  be  a  solution  of  the  inhomogeneous  Sturm-Louville  Equation: 

-  p(x)  ^j-  G(x,  x':A)  +  (q(x)  -  Aw(x))G(x,  x':  A)  = -S(x  -  x')  (3.1 6) 

subject  to  the  boundary  conditions. 

«i G(a,x':  A)  +  px-^-G{a,x':  A)  =  0  (3.17) 

a2G(b,x ':  A)  +  P2  ^-G{b,x':  A)  =  0  (3. 18) 


Let  %{x)  and  y/2(  x)  be  a  pair  of  solutions  of  the  differential  equation: 


“  +  (q(x)  -  Awix))y/Xx)  =  0 

(3.19) 

and  the  boundary  condition: 

aiW\(a)  +  fii~Vi(a) =  o 

ax 

(3.20) 

<x*V2(b)  +  p2-^-y/2(b)  =  0 
dx 

(3.21) 
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at  x  =  a,  and  x  =  b,  respectively.  Then  the  Green’s  Function  can  be  represented  in  the  form: 


G(x,x':  X)  =  - 


^iQ<)^2(x>) 

W(x') 


(3.22) 


in  terms  of  the  functions  y/x(x)  and  y/2(x) ,  where 


W(x)  =  y i(x) 


dy2{x) 

dx 


dx 


^2  00 


is  the  Wronskian  of  these  two  functions. 


(3.23) 
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4.  GREEN’S  FUNCTION  FOR  A  RANGE  INDEPENDENT  OCEAN 


This  section  gives  a  brief  overview  of  the  construction  of  the  Green’s  Function  for  a 
range  independent  waveguide.  One  must  adopt  a  cylindrical  coordinate  system  where  the  upper 
surface  is  given  by  the  plane  z  =0,  and  the  z-axis  is  oriented  downward.  The  source  is  located  at 
a 

depth  zs  at  a  horizontal  range  r  =  0. 

The  Green’s  Function  for  the  acoustic  wave  equation  in  a  constant  density  ocean  is  a 
solution  of  the  differential  equation: 

(-jT+“-7'  +  -T^T+'7T +k2)G(r,rs)  =  --S(r-rs)S(z-zs)S(<p-0s)  (4.1) 

dr  r  dr  r  d<p  dz  r 

subject  to  appropriate  boundary  conditions. 

Make  the  simplifying  assumption  that  the  speed  of  sound  is  a  function  of  the  depth 
coordinate  only.  Make  the  further  simplifying  assumption  that  the  source  is  omni-directional  and 
located  a  horizontal  range  r  =  0.  In  this  case,  Green’s  Function  simplifies  to  the  following 
differential  equation: 

(-77  +  +  k^2  )G(r’  z’z*)  =  “-  S{r)8{z  -  zs )  (4.2) 

dr  r  dr  dz  2m 

Here 


=  0  ,<p,<ps,z,zs)d</> 


(4.3) 


is  the  relationship  of  the  three-dimensional  Green's  Function  of  Equation  4.1  to  the 
two-dimensional  Green’s  Function  of  Equation  4.2. 

Neglecting  the  inhomogeneous  term  on  the  right  hand  side  of  Equation  4.2,  the  left  hand 
side  of  the  differential  equation  is  separable.  Define  the  one-dimensional  Green’s  Functions  for 
the  radial  and  depth  coordinate  as  solutions  of  the  following  differential  equations  for  a  given 
separation  parameter  A,: 

4tG>  :  X)  +  ~Gr(r :  X)  +  XGr(r :  X)  =  ~S(r)  (4.4) 

dr  r  dr  2m 
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—£ Gz(z,zs  :X)  +  (k(z)2  -X)G.(z,zs :  X)  =  -5(z - z ,) 


(4.5) 


The  two-dimensional  Green’s  Function  G(r,z,zs) can  be  written  as  the  spectral  integral 
G(r, z, z, )  =  -^  jGr (r :  X)G, (z, zt  :  X)dX  (4.6) 


of  the  radial  and  depth  Green’s  Function  over  an  appropriate  Contour  C.  The  Contour  C  can  be 
either  chosen  to  enclose  all  the  poles  of  the  radial  Green’s  Function  or  the  depth  Green’s 
Function  in  a  positive  sense,  that  is,  the  Contour  C  is  chosen  such  that  either  the  spectral  integral 


Gr(r:X)dX  =  —S(r) 

2m  l  Ijtr 


(4.7) 


over  the  radial  Green’s  Function  is  equal  to  the  Dirac  Delta  Function  with  respect  to  range,  or  the 
spectral  integral 


\Gz  (z, zs :  X)dX  =  S(z~zs) 


(4.8) 


over  the  depth  Green’s  Function  is  equal  to  the  Dirac  Delta  Function  with  respect  to  depth.  The 
fact  that  such  a  contour  may  be  chosen  is  due  to  the  fact  that  both  the  radial  and  depth  dependent 
Green’s  Functions  are  the  Green’s  Functions  of  a  Sturm-Louville  Differential  Equation. 


Now  one  may  make  the  change  of  variables: 
q2  =  X 

One  can  express  the  spectral  integral  for  the  Green’s  Function  as  the  integral: 
G(r, z,zs)  =  —  fGr(r  : q)Gz ( z , zs  :  q)qdq 

m  • 


(4.9) 


(4.10) 


»+/0 
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5.  RADIAL  GREEN’S  FUNCTION 


This  section  presents  a  derivation  of  the  radial  Green’s  Function. 

The  radial  Green’s  Function  for  a  range  independent  waveguide  satisfies  the  following 
differential  equation, 

^TGr(r:q)+~G,(r:q)  +  q2Or(r:q)  =  ~S(r)  (5.1) 

dr  r  dr  2  Ttr 

subject  to  the  boundary  conditions 

Km  r_*[r-j-Gr(r)]  =  --^- 
dr  in 

Hm„_  J~rd-GXr)-iqGXr))  =  0 

dr 

at  the  origin  and  at  infinity.  The  boundary  condition  at  the  origin  follows  from  integrating 
Equation  5.1  in  an  infinitesimal  region  about  the  origin.  The  boundary  condition  at  infinity 
follows  from  the  requirement  that  the  Green’s  Function  asymptotically  approaches  an  outgoing 
plane  wave. 

Equation  5.1  is  the  Bessel  Differential  Equation  of  Order  0.  Thus,  the  Green’s  Function 
is  of  the  form: 

G,(.r-.q)  =  AH"\qr)  +  BHf\qr) 

where  H^iqr)  and  H{02)(qr)  are  the  Hankel  Functions  of  the  first  and  second  kind.  The 

boundary  condition  at  infinity  allows  us  to  discard  the  term  with  the  Hankel  Function  of  the 
second  kind.  The  boundary  condition  at  the  origin  determines  the  coefficient  A.  The  radial 
Green’s  Function  is  given  by  the  following  expression: 
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6.  CONSTRUCTION  OF  THE  DEPTH-DEPENDENT  GREEN’S  FUNCTION 


This  section  describes  the  construction  of  the  depth  dependent  Green’s  Function  from 
solutions  of  the  homogeneous  differential  equation. 

The  depth  dependent  Green’s  Function  is  a  solution  of  the  following  differential  equation. 

LG2 (z  =  0, zt ; q)  s  P(z)^~  1  dG2(z  zs,q)  +  ^2^_q2 ^ ^ =  -£(z  -Zj)  (6. 1 ) 

dz  p(z)  dz 

The  Green’s  Function  is  subject  to  the  following  boundary  conditions. 

B,G2(z  =  0 =  fT(q)G1(z  =  =  0  (6.2) 

p{  0)  dz 

B2G2(z  =  d,zs;q)  =  f\q)G,(z  =  d,z,-,q)+^-^^^  =  0  (6.3) 

p(d)  dz 

The  Green’s  Function  satisfies  the  following  conditions  in  the  limit  z  approaches  the 
source  depth: 

lim£^0(G2(zJ  +  £,zs)~  G2(zs  - e,zs ))  =  0  (6.4) 

3g,(yg,».)-3^fr-8.»,)).  (6.5) 

dz  dz 

Let  p,(z)  and  p2(z)  be  a  pair  of  solutions  of  the  homogeneous  differential  equation  that 
satisfy  the  boundary  condition  at  the  surface  and  bottom  respectively,  that  is,  they  satisfy  the 
differential  equation 

LPx(z)  =  0  (6.6) 

Lp2(z)  =  0  (6.7) 

and  boundary  condition 

BlPl(z  =  0)  =  0  (6.8) 


B2P2(z  =  d)  =  0 
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The  depth  dependent  Green’s  Function  is  given  by  the  following  expression: 


G2(z,zs;q)  =  - 


li(z<)/?2(z>) 

W{zs-q) 


where 


is  the  Wronskian  of  the  two  solutions,  and 
z<  =  min(z,zs) 

z>  =  max(z,zt). 


(6.10) 


(6.11) 


(6.12) 

(6.13) 
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7.  NORMALIZATION  OF  DEPTH  FUNCTIONS 

This  section  describes  the  normalization  condition  of  the  depth  functions  and  the  depth 


dependent  Green’s  Function  following  the  discussion  by  Porter4. 

The  depth  function  is  a  solution  of  the  following  differential  equation: 

L{q)\ff{z)  =  p(z)j-(  )  ~  3 ^  +  {k\z)  -  q2)l//(z)  =  0  (7.1) 

dz  p{z )  dz 

and  the  following  boundary  conditions: 

B,(qMz  =  0)  =  /r(«)r(z  =  0)+  ^0))ay&=”'  =  °  P'2) 

BMv(z  =  <0  =  f\qw(z  =  =  0  (7,3) 

p(d)  dz 

In  general,  a  non-zero  solution  of  the  above  differential  equation  and  boundary  conditions 
only  exists  for  a  discrete  spectrum  of  eigenvalues. 

Let  pt(z)  and  p2(z)  be  a  pair  of  unforced  solutions  of  the  differential  equation: 

L(q)pjz)  -  0  (7.4) 

Uq)pt(z)  =  0  (7.5) 

and  boundary  condition  at  the  surface  (bottom),  respectively: 

B,(q)p,(z  =  0)  =  0  (7.6) 

BMPi^  =  d)=  0  (7.7) 


Here,  boundary  conditions  are  only  imposed  on  px{z)  and  p2(z)  at  one  end  of  the 
domain,  respectively.  This  relaxation  of  the  boundary  condition  at  one  end  of  the  domain  allows 
the  existence  of  a  continuous  spectrum  of  solutions.  The  functions  px{z)  and  p2(z)  are 
analogous  to  the  left-going  and  right-going  traveling  waves  on  a  one-dimensional  string.  The 
depth  functions  are  analogous  to  the  standing  waves  on  this  string  when  both  ends  of  the  string 
are  clamped. 
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Let 


tV(z,g)  =  Pl(z)^l-^lpAz)  (7.8) 

be  the  Wronskian  of  the  above  pair  of  solutions  of  this  differential  equation.  The  two  solutions 
are  said  to  be  linearly  independent  if  the  Wronskian  is  non-zero. 


Recall  the  depth  dependent  Green’s  Function  may  be  represented  in  the  form 


G2(z,zs;q)  =  - 


Pi(OP2(z>) 

W{xM\q) 


(7.9) 


in  terms  of  this  pair  of  solutions  and  the  Wronskian.  The  spectrum  of  normal  modes  is 
comprised  of  those  values  of  q  for  which  the  above  Wronskian  vanishes: 


W{z,\q)  =  0 


(7.10) 


Let  y/(z)  be  the  depth  function  at  eigenvalue  qm  at  which  the  above  Wronskian  vanishes. 
Without  loss  of  generality,  one  may  assume  p,(z)  and  p2(z)  form  a  one-parameter  family  of 
solutions  of  the  depth  equation  that  satisfy  the  limiting  conditions: 

KnW/>,(z)  =  !Kz)  (7.H) 

li«W  PiW  =  Mz)  (7.12) 

that  they  approach  the  depth  function  y/(z)  as  q  approaches  the  eigenvalue  qm . 


The  following  relationship: 


P2(z)L(<lJV(z)-(L(q)P2(z))V(z)  = 


(7.13) 


between  y/{z)  and  p2(z)  follows  from  the  fact  these  two  functions  are  solutions  of  the  depth 
equation  for  eigenvalue  qm  and  q,  respectively.  Integrating  the  above  quantity  over  the  water 
column,  one  arrives  at  the  expression 


((AW 


jjKf) 

dz 


dp2(z) 

dz 


y/(z))l  p(z))\ 


D 


L-=o 


D 

-(q2  -ql)  l p2(z)t//(z)/ p(z)dz 
0 


(7.14) 
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Similarly,  the  relationship: 

p,  (z)L(q)p2  ( z)-(L(q)pl  ( z))p2  (z)  =  ((p,  p2  (z))  /p(z))  =  0  (7. 1 5) 


proves  that  the  ratio 


W(z;q)/p(z)  =  (p,(z)%^-^p2(z))/p(z) 

dz  dz 

of  the  Wronskian  and  the  density  is  a  constant. 

Using  the  relationship: 

dz  dz  ^(z)  dz  p2(z)  dz 


and  the  boundary  conditions 

1  dyrjz  =  £>)  _  /8(gm) 
yr(z  =  £>)  dz  g8(gm) 

one  arrives  at  the  expression 


ft  (z>)  *z>) = (^2  ■  (z,)r(0) 

9z  dz  £  (?)  £  (O 


(7.16) 


(7.17) 


(7.18) 


(7.19) 


Since  the  depth  function  y/(z)  satisfies  the  boundary  condition  at  both  the  top  and  bottom 
surface,  one  has  the  following  expression  for  the  impedance  of  the  depth  function  at  the  top 
surface: 


1  3^(0),  f±gJ 

w(0)  dz  gT(qm) 


pi  0) 


(7.20) 


Utilizing  the  fact  that  the  ratio  of  the  Wronskian  and  density  is  a  constant  one  obtains  the 
following  relationship  for  the  impedance  of  the  function  p2(z)  at  the  top  surface: 


1  dp2(0)  _  1  dp,(0)  |  fmq)  _  fT(t ?)  Wjz^q) 

p2(  0)  dz  p,(0)  dz  p(0)  gTiq)F  p(z)p,(0)p2(0) 


(7.21) 
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Using  Equations  7.17, 7.20,  and  7.21  one  arrives  at  the  following  expression: 
dz  dz 

= (yy  -  +  OT* : «)  -  IK(r :  MO)/  p,  (0) 

Here,  use  of  the  vanishing  of  the  Wronskian 
w(z:qm)  =  0 

at  the  eigenvalue  qm  has  been  made  to  add  a  term  that  is  proportional  to  this  term  to 
Equation  7.22. 


(7.22) 


(7.23) 


Substituting  Equations  7.19  and  7.21  into  Equation  7.14,  one  arrives  at  the  following 
expression  for  the  normalization  of  the  depth  functions: 


(7.24) 


Divide  this  equation  by  (q2  —q2m)  and  take  the  limit  q  approaches  qm  to  obtain  the 
following  relationship  between  the  Wronskian  and  the  normalization  of  the  depth  functions. 


1  gfjJ 

p(z)  dz 


u  , 

=  2?m  j 


¥(z)W(z)  ,  a  ,fT(gJ 


dz-T-(^rr rM°)  +-f(- 


a?  gr(?m) 


a*  Vfoj 


)^(D)2 


(7.25) 


Without  loss  of  generality,  one  can  normalize  the  depth  functions  such  that  the  partial 
derivative  of  the  Wronskian  at  the  m’th  normal  mode  is  given  by  the  expression: 

(7.26) 

by  requiring  the  following  normalization  condition  on  the  depth  function: 

V£M£)  J  Vr(*->M0)»  +  _ .  ,, , 


(7.27) 
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In  the  case  of  a  pressure  release  surface  and  a  rigid  bottom,  the  partial  derivative  terms  in  the 
above  expression  vanishes  and  the  previous  expression  reduces  to  the  usual  normalization 
condition  in  the  case  the  spectrum  is  discrete. 
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8.  N  HOMOGENEOUS  FLUID  LAYERS  OVER  A  RIGID  BOTTOM 


This  section  describes  the  construction  of  the  Green’s  Function  for  N  fluid  layers  over  a 
rigid  bottom. 

The  Green’s  Function  is  given  by  the  following  spectral  integral: 


G(r,z,zs ) 


An 


jH(0')(qr)G!(z,zs  :q)qdq 


(8.1) 


where  G,(z,zs :  q)  is  the  depth  dependent  Green’s  Function.  This  spectral  integral  can  be 

represented  as  a  sum  of  residues  of  poles  of  the  depth  dependent  Green’s  Function  in  the  upper 
half  q-plane.  The  remainder  of  this  section  will  be  devoted  to  solving  the  characteristic  equation 
for  the  location  of  these  poles,  and  the  evaluation  of  the  depth  functions. 

Let  {p„,cn :  n  -  0,1  ..N  - 1}  denote  the  density  and  sound  speed  in  the  N  layers.  Let 
{zn  :n  =  0,1...  V}  denote  the  z-coordinate  of  the  N+l  interfaces.  Without  loss  of  generality,  one 
may  assume  z0=  0,  and  zN  -  D. 


The  depth  function  in  the  n’th  layer  is  required  to  satisfy  the  differential  equation: 


^-  +  (k2„~q2)  =  0  (8.2) 

az 

where  kn  =  col cn  is  the  wavenumber  in  the  n’th  layer.  The  functions: 

Fn+(z)  =  e+ih’(z~s”)  (8.3) 

F;(z)  =  e-rt”(— >  (8.4) 

denote  a  basis  for  the  downward  and  upward  going  waved  in  the  n’th  layer,  where 

h„(q)  =  i^[q^-k2n  (8.5) 

is  the  vertical  wavenumber  in  the  n’th  layer,  and  q  is  the  horizontal  wavenumber.  Here  one 
adopts  the  conventions  of  Ewing,  Jardetsky,  and  Press6  for  the  vertical  wavenumber;  that  is,  the 
imaginary  component  of  the  vertical  wavenumber  is  greater  than  zero  on  the  physical  sheet,  and 
the  branch  cut  is  the  hyperbola  on  which  the  imaginary  component  is  zero  in  the  complex  q- 
plane. 
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Let  the  following  sum  represent  the  depth  function  in  the  n’th  layer: 

Fn  (z)  =  AnFn+{z)  +  A~nF;(z)  (8.6) 

of  upward  and  downward-going  waves.  The  depth  function  is  required  to  satisfy  the  following 
boundary  conditions. 


^o(^o)  =  0 


F«  (2-1 )  =  Fn+i  (2-, ) ,  n  =  0,1  ..(N  -  2) 

1  g&J  1  ^-,(2^) 


Pn 


n+l 


<7z 


,  »  =  0,l..(tf-2) 


(8-7) 

(8.8) 

(8.9) 


dJjM  (zh 1 )  _  q 


(8.10) 


The  above  boundary  conditions  specify  the  upper  surface  to  be  a  pressure  release  surface, 
and  the  bottom  surface  a  rigid  surface.  The  remaining  boundary  conditions  are  continuity  of 
pressure  and  the  normal  velocity  between  adjacent  layers. 

Equation  8.7  implies  the  following  condition 

Ao  =  ~Ao  (8.11) 

on  the  depth  function  in  the  top  layer.  Equation  8.10  implies  the  following  condition 

A~N-\=+AN-\e*2,hs''dN''  (8.12) 

on  the  depth  function  in  the  bottom  layer.  The  remaining  boundary  conditions  imply  the 
relationship: 


*n+l 


MA 


where  An  is  the  column  vector 


A.  = 


J 


and  Mn  is  the  matrix 


M. 


K 

<2 

¥1 

Ml, 

(8.13) 


(8.14) 


(8.15) 
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M->  _  1  (  M'n  -^",1 

(8.16) 

Mnu  =  (1  +  Pn+X  h”  )e+iK<  /2 

Pn  K+\ 

(8.17) 

9 

M”2  =(1  Pn+i  h”  )e~h"d"  /  2 

P„  ^„+i 

(8.18) 

- 

M"  =  (1  /7"+1  hn  e+Mn)/2 

A.  ^,+1 

(8-19) 

M"n  =  (l+Pn+l  h"  )e~ih"d"  H 

Pn  K+\ 

(8.20) 

Here, 

dn  =  (Zn+ 1  "2J>  «  =  0,1...// -1 

(8.21) 

is  the  thickness  of  the  n’th  layer. 

The  characteristic  equation  for  the  normal  modes  is  the  transcendental  equation: 

dF„{zN)  =ihff  [e+/w»-,  _ A-_le~ihN-'dN-'^=  0 

(8.22) 

where  the  depth  function  in  the  bottom  layer  is  determined  by  the  following  relationship  between 
the  depth  function  in  the  top  layer  with  that  in  the  bottom  layer: 

An_\  —  Mn_2...M0A0 

(8.23) 

- 

a  -4+1l 

A)  “  ^0  J 

V  V 

(8.24) 

* 

The  above  relationship  determines  the  depth  functions  throughout  the  waveguide  up  to 
the  multiplicative  factor  A„ .  This  factor  is  determined  by  normalizing  the  depth  functions 
according  to  the  following  condition: 

Yiz)Fiz)<b= yf'(2)f-(z)<*=i 

o  P(Z )  ^  J  Pn 

(8.25) 
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Let  {qm}  be  a  complete  basis  for  the  normal  modes  in  the  complex  q-plane,  and 
F(z  '■  qm )  the  corresponding  depth  function.  Then  the  Green’s  Function  for  this  waveguide  is 
equal  to  the  following  summation: 


G(r,z,zs ) 


i 

4p(0 


■(in,) 


(8.26) 


where  the  depth  functions  satisfy  the  following  orthonormality  condition: 

DtF(z-(im)F(z-qn)±  =  gn 

oJ  />(*) 


(8.27) 


An  alternative  method  of  computing  the  Green’s  Function  is  to  directly  evaluate  the 
contour  integral  given  in  Equation  8. 1 .  This  is  the  approach  used  in  fast  field  programs  such  as 
“SAFARI”  and  “OASES”.  This  approach  requires  the  construction  of  the  depth  dependent 
Green’s  Function  off  mass  shell.  It  has  the  advantage  that  one  doesn’t  have  to  solve  the 
transcendental  equation  given  by  Equation  8.22  for  the  normal  modes. 

In  order  to  construct  the  off  mass  shell  solution  for  the  depth  dependent  Green’s  Function 
one  constructs  the  functions  pt(z)  and  p2  (z)  according  to  the  following  prescription.  First, 
define  the  quantities: 

z=mm(z,zs)  (8.28) 

z>=  max(z,  zs)  (8.29) 


Suppose  the  coordinates  z<  and  z>  are  located  in  the  n  and  n’  layers  respectively,  that  is, 
These  coordinates  satisfy  the  inequalities: 


Define  the  coefficients: 


r+fi 

Ur"-"-"- 

l-J 

II 

N-\  ^ 

j 

1 

N-2 

\7e  J 

(8.30) 

(8.31) 


(8.32) 


(8.33) 
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and  define  the  functions: 

p,(z)=a;f;(z)+a;f;(z) 

(8.34) 

(8.35) 

These  functions  satisfy  the  homogeneous  equation  for  the  depth  functions,  and  the 
boundary  condition  at  the  upper  and  lower  surface,  respectively.  The  depth  dependent 
Green’s  Function  is  given  by  the  following  expression  in  terms  of  these  two  functions: 


Gz(z,zs:q)  =  - 


W(px,p2)(zs) 


where 


W(pvp2)  =  px(z) 


dp2(z) 

dz 


dPjjz) 

dz 


Pl(z) 


(8.36) 


(8.37) 


is  the  Wronskian  of  this  pair  of  functions. 


As  a  means  of  comparing  the  above  contour  integral  representation  with  the  normal  mode 
representation  of  the  Green’s  Function,  the  normal  mode  representation  is  given  by  the  following 
expression  in  terms  of  the  depth  functions  of  the  normal  modes. 


G:{z,zs  :q)  = 


1  y  F(z:gm)F(zt  :gm) 
p{z^  (q2-q2m) 


(8.38) 


Substituting  the  above  expression  for  the  depth  dependent  Green’s  Function  into 
Equation  8.1,  one  arrives  at  the  normal  mode  expression  for  the  Green’s  Function  given  in 
Equation  8.26. 
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9.  N  HOMOGENEOUS  FLUID  LAYERS  OVER  A  HOMOGENEOUS  HALF-SPACE 


This  section  describes  the  modifications  of  the  previous  section,  when  the  rigid  bottom  of 
the  previous  section  is  replaced  by  a  homogeneous  half-space. 

One  of  the  major  differences  between  the  case  of  a  homogeneous  half-space  and  a  rigid 
bottom  is  that  in  the  case  of  a  half-space  there  are  generally  a  finite  number  of  normal  modes, 
where  as,  in  the  case  of  a  rigid  bottom  there  are  generally  an  infinite  number  of  normal  modes. 

In  addition,  the  integrand  of  the  equation 


1  -H»-/0 

G(r,  z,zs)  =  —  jH(0l)  ( qr)Gz  (z,  zs  :  q)qdq 
™  -~+/0 


(9.1) 


is  generally  not  an  even  function  of  the  vertical  wavenumber  of  the  homogeneous  half-space. 
This  results  in  a  branch  cut  contribution  of  the  spectral  integral  from  the  integration  of  the 
spectral  integral  around  the  branch  cut  of  the  vertical  wavenumber  of  the  homogeneous  half¬ 
space  in  the  complex  q-plane.  The  half-space  also  affects  the  normalization  of  the  normal 
modes. 


Let  {p„,cn :  n  =  0,1  ..N  - 1}  denote  the  density  and  sound  speed  in  the  N  layers.  Let 
{zn :  n  =  0,1  ...N}  denote  the  z-coordinate  of  the  N+l  interfaces.  Without  loss  of  generality,  one 
assumes  z0  =  0 ,  and  zN-D.  Let  {pN ,  cN }  denote  the  density  and  speed  of  sound  in  the 
homogeneous  half-space. 

The  depth  function  in  the  n’th  layer  is  required  to  satisfy  the  differential  equation: 

+  =  0  (9.2) 


where  kn  =  col  cn  is  the  wavenumber  in  the  n’th  layer.  The  functions: 
Fn+(z)  =  e+ih'^) 


(9.3) 


Fn-(z)  =  e-ih^)  (9.4) 

denote  a  basis  for  the  downward  and  upward  going  waved  in  the  n’th  layer,  where 

K(q)  =  i4ci2  ~kl  (9-5) 


is  the  vertical  wavenumber  in  the  n’th  layer,  and  q  is  the  horizontal  wavenumber. 
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One  may  again  represent  the  depth  function  in  the  n’th  layer  as  the  following  sum 
Fn(z)=A;F;(z)  +  A-nFn-(z)  (9.6) 

of  upward  and  downward  going  waves.  The  depth  function  is  required  to  satisfy  the  following 
boundary  conditions: 


F0(z0)  =  0 


Fn (*»♦. )  =  Fn+X (zn+) ),  «  =  0,1..(A^- 1) 


1  1  <+,(zB+l) 

Pn  P„+|  dz 


n  =  0,l..(Af-l) 


(9.7) 

(9.8) 

(9.9) 


Assuming  the  source  is  above  the  homogeneous  half-space,  one  requires  the  condition 
that  the  coefficient: 


a~n-°  (9.10) 

of  the  upward  moving  waves  in  the  homogeneous  half-space  vanish. 

Equation  9.7  implies  the  following  condition 

Ao  =  ~Ao  (9.11) 

on  the  depth  function  in  the  top  layer.  Boundary  conditions  shown  in  Equations  9.8  and  9.9 
imply  the  relationship. 


An+ 1  =  MnAn 

where  An  is  the  column  vector 


(9.12) 


A 


n 


f 

V 


(9.13) 


and  M n  is  the  matrix  given  by  Equations  8.15  through  8.20  in  the  previous  section. 

Equation  9.10  in  conjunction  with  Equations  9.8  and  9.9  imply  one  can  replace  the  boundary 
condition  for  the  rigid  bottom  by  the  following  boundary  condition: 


:«)  + 


8*  fa)  d 

Pn- i  dz 


fn-i(zn  :?)  =  0 


(9-14) 
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where  the  coefficients  fB(q)  and  gB(q)  are  given  by  the  following  expressions: 


/*(*)= -iM?) 

(9.15) 

Pn 

gB(q)  =  l 

(9.16) 

The  depth  functions  {F(z  :qm)} are  normalized  by  the  following  condition: 


J 


p(z)  2 qm  dq  gB(qmy 

The  Green’s  Function  is  given  by  the  following  residue  term  and  cut  contribution: 


(9.17) 


G(w,)  =  — £  H' "l(qmr)F(z  q. )F(z, :  q, )  +  Gt, (r.z.z, ) 
Here,  the  cut  contribution  is  given  by  the  following  expression: 
Gcul(r,z,zs)  =  -^  \H"\qr){G2(z,zs  :+hN)-G2(z,zs  :-hN)}hNdhN 


(9.18) 


(9.19) 


The  above  integral  is  an  integral  over  the  branch  cut  for  the  homogeneous  half-space, 
where  the  integrand  is  proportional  to  the  difference  in  the  depth  function  across  the  cut. 
Evaluation  of  the  above  integral  requires  the  construction  of  the  off  mass  shell  representation  of 
the  depth  dependent  Green’s  Function. 

In  order  to  construct  the  off  mass  shell  solution  for  the  depth  dependent  Green’s  Function 
one  must  construct  the  functions  px{z)  and  p2 (z) according  to  the  following  prescription.  First, 
define  the  quantities: 


z<=  min  (z,zs) 
z>=  max(z,zr) 


(9.20) 

(9.21) 


Suppose  the  coordinates  z<  and  z>  are  located  in  the  n  and  n’  layers  respectively,  that  is, 
these  coordinates  satisfy  the  inequalities: 


<  z<  <  z 


n+1 


*„•  <  *>  <  *„.+1 


(9.22) 

(9.23) 
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Define  the  coefficients: 


(a+) 

r+1l 

J 

...M0 

-i 

(a+] 

f1' 

nn' 

'A*  J 

,0, 

and  define  the  functions: 

P\  te) = K  f*  (z) + a;f;  (z) 
p2(z)  =  A;.F;.(z)  +  A-n.F;.(z ) 


(9.24) 

(9.25) 


(9.26) 

(9.27) 


These  functions  satisfy  the  homogeneous  equation  for  the  depth  functions,  and  the 
boundary  condition  at  the  upper  and  lower  surface,  respectively.  The  depth  dependent 
Green’s  Function  is  given  by  the  following  expression  in  terms  of  these  two  functions: 


G:(z,zs:q)  =  - 


PiMPz^) 

W(Pl>P2)(Zs) 


where 


W(.P\,p2)  =  px(z) 


dPi{z)  (friz) 


dz 


dz 


Pl(z) 


(9.28) 


(9.29) 


is  the  Wronskian  of  this  pair  of  functions. 
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10.  N  INHOMOGENEOUS  FLUID  LAYERS 


This  section  describes  the  modifications  of  the  previous  two  sections  in  the  case  the  index 
of  refraction  squared  is  a  piece-wise  linear  function  of  depth.  In  this  case,  the  complex 
exponential  functions  of  the  previous  sections  are  replaced  by  complex  Airy  functions.  In  the 
case  of  an  inhomogeneous  layer,  there  is  no  longer  a  natural  split  of  the  depth  functions  into 
upward  moving  and  downward  moving  components.  This  work  is  based  on  the  description  by 
Stickler7  on  the  computation  of  normal  modes  in  a  waveguide  with  a  piece-wise  linear  index  of 
refraction  squared. 

Let  {/?„ ,  cTn ,  cBn  :  n  =  0,1  ...(TV  - 1)}  denote  the  density,  sound  speed  at  the  top  of  the  layer, 
and  sound  speed  at  the  bottom  of  the  layer  for  the  N  inhomogeneous  layers,  where  the  index  of 


refraction  squared  is  the  following  linear  function  of  depth: 

c;2(z)  =  an-J3n(z-z„)  (10.1) 

«„=(cjr2  (10.2) 

fin=((cTnr2-(cTny2)/dtt  (10.3) 

Define  the  change  of  coordinates: 

z-*Z„(z)  (10.4) 

in  the  n’th  layer,  where  the  coordinate  function  Zn  (z)  is  required  to  satisfy  the  following 
constraint: 

=  -(CO1  /  c(z)2 -q1)  (10.5) 

dz 

The  coordinate  function  Zn  (z)  is  given  by  the  following  linear  function  of  z: 

Z„  (z)  =  -L\  {(O1  (a„  -  A  (z  -  z„ ))  -  g: 2 )  (10.6) 

where  the  coefficient  Ln  is  given  by  the  expression 

i.=(<»4A!)'''6  (10-7) 
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The  derivative  of  this  coordinate  function  is  given  by  the  expression: 

=*>2PA 


dZ„(z)  2  f2 


dz 

The  differential  equation  for  the  depth  function 
dFn{z) 


(10.8) 


dz 


+  (k2(z)-q2)Fn(z)  =  0 


(10.9) 


in  the  n’th  layer  is  transformed  into  the  following  differential  equation: 


dZtt 


-ZnF„ 


=  0 


(10.10) 


under  the  coordinate  transformation  z  ->  Zn .  Equation  10.10  is  the  Airy  Differential  Equation, 
which  implies  the  depth  function  is  given  by  the  following  linear  combination  of  Airy  Functions: 

Fn(z)  =  A+Ai(Zn(z))  +  A;Bi(Zn(z))  (10.11) 

Note  the  above  expression  for  the  depth  function  becomes  ill  conditioned  in  the  limit  the 
sound  speed  gradient  in  the  layer  vanishes,  that  is,  the  coefficient  pn  0  approaches  zero.  In 

this  case  the  coordinate  transformation  z  ->  Zn  becomes  ill  defined. 

The  ill-conditioned  nature  of  this  coordinate  transformation  for  infinitesimally  small 
sound  speed  gradients  can  lead  to  numerical  overflow  errors  in  the  computation  of  the  depth 
functions.  Therefore  it  is  necessary  to  place  a  minimum  limit  on  the  sound  speed  gradient  for 
which  the  above  linear  representation  of  the  index  of  refraction  squared  will  be  used.  If  the 
sound  speed  gradient  is  below  this  minimum,  the  sound  speed  profile  in  the  layer  is  considered  to 
be  a  constant,  and  the  depth  functions  of  the  previous  two  sections  are  used. 

Thus  one  introduces  the  following  basis  for  the  depth  functions  if  the  sound  speed 
gradient  is  sufficiently  small: 

Fn+(z)  =  e+ihj2~z") 

F-(z)  =  e~ihAz-^ 


(10.12) 

(10.13) 
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where  the  sound  speed 

cn=(cT„  +cb„)/2  (10.14) 

is  equal  to  the  average  of  the  speed  of  sound  at  the  top  and  bottom  surface,  and 

*„(<?)= <10-15) 

is  the  vertical  wavenumber.  Otherwise,  we  use  the  following  basis: 

(1016> 

F,-(z)  =  Bi(Z.(z))  (10.17) 

for  the  depth  functions,  provided  the  gradient  {ir!  is  sufficiently  large. 

Let 


1 (io.i8) 

dz  dz 

denote  the  Wronskian  of  the  above  basis  for  the  depth  functions.  In  the  case  of  vanishing 
gradient  (iso- velocity  layer),  this  Wronskian  is  given  by  the  expression: 


Wn{z:q)  =  -2ihn(q) 


(10.19) 


where  hn  (q)  is  the  vertical  wavenumber  of  this  iso-velocity  layer.  Otherwise,  the  Wronskian  is 
equal  to  the  expression: 


where  the  expression 

/4,(x)^w_^waw=i 

dx  dx  k 

is  used  for  the  Wronskian  of  the  Airy  Functions. 


(10.20) 


(10.21) 


The  depth  function  Fn  (z:q)  in  the  n’th  layer  may  be  represented  by  the  following  linear 
combination  of  the  functions  F„+  ( z :  q)  and  F~{z :  q) : 


Fn(z:q)  =  A:F;(z:q)  +  A;Fn-(z:q) 


(10.22) 
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The  depth  function  in  the  top  layer  satisfies  the  following  constraint  for  a  pressure  release 
surface: 


^o(z0)  =  0 


(10.23) 


In  addition,  the  depth  function  must  satisfy  the  constraint  of  continuity  of  pressure  and 
normal  velocity  at  the  interface  between  adjacent  layers: 

C^n+1  )  P"n+ 1 0'n+l  )  >  ^  —  0,1  ..(N  1)  (10.24) 


1  dF„(z„+i) 

Pn 


1  dF^,  |(zn+l) 
Pn*  1  dz 


n  =  0,l..(Ar-l) 


(10.25) 


In  the  case  of  a  rigid  basement,  the  depth  function  satisfies  the  boundary  condition: 


jjjM  (zv  )  _  q 
dz 


(10.26) 


at  the  bottom  surface.  In  the  case  of  a  homogeneous  half-space,  the  coefficient  for  an  upward 
moving  wave  in  the  basement  is  required  to  vanish  and  the  depth  function  is  required  to  satisfy 
the  constraints  of  continuity  of  pressure  and  normal  velocity  at  the  interface  with  the 
homogeneous  half-space: 


a~n=o 

^jV-I  ( ZN  )  =  Fn  (ZN  ) 

1  dFfl_l(ztl)_  1  dFN(zN) 

P  N—\n  dz  pN  dz 


(10.27) 

(10.28) 

(10.29) 


The  constraint  on  the  depth  function  by  the  presence  of  the  bottom  can  be  reformulated 
by  a  boundary  condition  of  the  following  form: 


Pn- i  dz 


(10.30) 


where  the  coefficients  fB(q )  and  gB(q )  are  given  by  the  following  expressions  in  the  case  of  a 
homogeneous  half-space: 


f\q ) 


zihM 

Pn 


(10.31) 


£*(<?)  =  1 


(10.32) 


10-4 


CSS/TR-02/10 


and  the  following  expressions  in  the  case  of  a  rigid  bottom: 


/'(*)  =  0  (10.33) 

gB(q)  =  l  (10.34) 


The  depth  functions  {F(z  :  <7m)}  are  normalized  by  the  following  condition: 


A+_l_  A(/lM)F(D-9_y  =1 


P(z) 


^  g  (qmY 


(10.35) 


Equation  10.23  implies  the  following  relationship  for  the  coefficients  of  the  depth 
function  in  the  top  layer: 


A~  —  —A+ 


(10.36) 


Equations  10.24  and  10.25  imply  the  following  relationship  between  the  coefficients  of 
adjacent  layers: 


An+ 1  =  MnAn 


Here  A„  is  the  coefficient  vector 


A.  = 


and  Mn  is  the  propagator  matrix 
Mn 


M_  = 


\M 2\  M 22  J 


whose  matrix  elements  are  given  by  the  following  relations 
1 


M ,T=  + 


M™=  + 


Mm  =- 
"*21 


^+1(*n+I) 

1 

1 

Wn+,{^) 


{K(zn+i) 

iFn(Zn+l) 


dFn+]  (Zn+ 1  ) 

Pn+\ 

dFn+(Zn«) 

dz 

Pn 

dz 

dFnAZn,\) 

Pn+\ 

dFn(Zn+\) 

dz 

Pn 

dz 

dFn+l  (Zn+\  ) 

Pn+\ 

dFn+(Zn+i) 

dz 

Pn 

dz 

Fn+l  (“ni  l  )} 

■^n+|(^n+l)} 
^n+1  (Zn+\ )} 


(10.37) 


(10.38) 


(10.39) 


(10.40) 

(10.41) 

(10.42) 
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M”=- 


\  f P~(-  \dFn+\(Zn+])  P n+l  dFn  (zn+ 1)  c,+  /_  ^ 

KAzn+i) 1 '  (  "+,)  dz  -  pn 


(10.43) 


Recall  that  the  expression: 


K(z-q)  =  K(z) 


dK(z) 

dz 


dK(z) 

dz 


FAz) 


(10.44) 


is  the  Wronskian  of  the  basis  for  the  depth  functions  in  the  n’th  layer. 

Using  Equations  10.36  and  10.37,  the  depth  functions  in  all  the  layers  are  determined  up 
to  a  common  multiplicative  factor.  However,  the  depth  functions  constructed  from 
Equations  10.36  and  10.37  do  not  necessarily  satisfy  the  boundary  condition  in  Equation  10.30  at 
the  bottom  layer.  This  equation  acts  as  the  characteristic  equation,  whose  solutions  determine  the 
location  of  the  normal  modes  in  the  complex  q-plane. 


Given  the  above  construction  of  the  coefficients  of  the  depth  function,  the  construction  of 
the  depth  dependent  Green’s  Function  and  the  spectral  integral  proceeds  in  analogy  to  the 
previous  two  sections. 
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11.  EXAMPLE:  ISO-VELOCITY  LAYER  OVER  A  RIGID  BOTTOM 


This  section  describes  an  iso- velocity  layer  over  a  rigid  bottom. 

The  characteristic  equation  for  the  normal  modes  in  this  case  is  given  by  the  equation: 
cos(h0d)  =  0  (11.1) 

where 

A„=i'V?!-*o2  <"-2) 

is  the  vertical  wavenumber.  The  solutions  of  this  equation  are  given  by  the  modes  vertical 
wavenumber: 


*,(»„)  =  (<■  + 1/2)| 

The  normalized  depth  function  is  given  by  the  following  relationship: 
F(z-.qm)  =  J^sin(h0(q„)d) 

The  Green’s  Function  is  given  by  the  following  modal  summation: 

G(r,  z,zs)  =  -+ 7  £  ( qmr )  sin(/z0  (qm  )z)  sin(/?0  (qm  )zs ) 

2  a  m 


(11.3) 


(11.4) 


(11.5) 
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12.  EXAMPLE:  PEKERIS  MODEL 


This  section  describes  the  Pekeris  Model,  that  is,  an  iso- velocity  layer  over  a 
homogeneous  half-space. 

The  normal  modes  for  the  Pekeris  Model  are  given  by  the  solutions  of  the  characteristic 
equation: 

mh0  cos (h0d)  -  ihx  sin (h0d)  =  0  (12.1) 

where  m  is  the  density  ratio. 


m  = 


Pi 


Po 


(12.2) 


The  depth  dependent  Green’s  Function  is  given  by  the  following  expression: 


G,(z,zs) 


sin(/t0z< )  (mh0  cos (/»0  (d  -  z>  ))  -  ihx  sin(ft0  (d  -  z>  ))) 
hQ  (mh0  cos(h0  (d))  -  ihx  sin (/?0  (d))) 


(12.3) 


The  depth  function  in  the  top  layer  is  given  by  the  following  expression: 
F(z :  q)  =  A*  sin (h0z) 

where  the  normalization  factor  is  given  by  the  expression. 

A+  mh0  2p(zs)h\ 

0  sin(V)  V  hxd(m2hl  -hx)  +  im(h20  -  h\ ) 


(12.4) 


(12.5) 


The  normal  mode  portion  of  the  Green’s  Function  is  given  by  the  following  sum  of 


residues: 
GMode(r,z,zs)  = 


(12.6) 


^  E  ^0°  (?/)  sin(/20z)  sin(hozs ) 


2"J 


r  _jA_y 

sin  (h0d) 


hxd(m2hl  ~h2)  +  im(h^  -h?) 


2  1.2  > 
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The  branch  cut  contribution  to  the  Green’s  Function  is  given  by  the  following  spectral 
integral: 


GCul(r,z,zs)  =  +^~  j ’hfdh^iqr) 
4  7t  J 


m  h 


2L2 


sm(h0z)sm(h0zs) 
cos(h0d)2  +/»,2  sin(h0d)2 


(12.7) 
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13.  EFFECT  OF  A  SMALL  WAVE-HEIGHT,  RANDOMLY  ROUGH,  TWO  FLUID 
INTERFACE  ON  THE  COHERENT  COMPONENT  OF  THE  REFLECTED  AND 

TRANSMITTED  FIELD 


This  section  describes  the  effects  of  a  randomly  rough  two-fluid  interface  with  small 
wave-height  on  the  coherent  component  of  the  reflected  and  transmitted  field8. 

Let  {/?, ,  c, }  and  {p2 ,  c2 }  denote  the  density  and  sound  speed  in  Regions  One  and  Two, 
where  the  plane  z  =  0  is  the  interface  between  Region  One  (z  >  0)  and  Region  Two  (z  <  0)  in  the 
unperturbed  case.  The  pressure  field  in  Regions  One  and  Two  are  solutions  of  the  differential 


equations: 

V2Pl+k2Pi=0  (13.1) 

V2p2+k22p2=0  (13-2) 

where  /?,  and  p2  are  the  pressure  fields  in  Regions  One  and  Two  respectively.  Here,  kt  and  k2 
are  the  wavenumber  in  Regions  One  and  Two  respectively.  The  pressure  field  is  subject  to  the 
following  boundary  conditions  at  the  interface  z  =0. 

P\  (^  =  0)  =  p2  (z  =  0)  (13.3) 

1  4?,(z  =  0)=  1  dp2(z  =  0 )  (134) 

/?!  dz  p2  dz 

Let 

Pi„c  =  exp[+/£0,x  •  r  -  iki0  ,z]  (13.5) 


be  the  incident  field,  where  k01  is  the  horizontal  wave- vector  of  the  incident  field  in  the  xy- 
plane,  and 

03.6) 

is  the  vertical  wavenumber  of  the  incident  field.  The  unperturbed  solution  for  the  pressure  field 
in  the  two  regions  is  given  by  the  following  expressions: 

p,  =  exp[+/£0  x  •  r  - ikm  .z]  +  R0  exp [+iko  x  •  r  +  ik^  .z]  (13.7) 
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Pi  =  T0  &xp[+ik0±  •?- ik20;z ] 
where 


(13.8) 


^10.:  ^  P\  ^20.:  ^  Pi 
k\0.:  !  P\  ■*"  ^20 ,r  ^  /^2 


(13.9) 


Tq -\  +  R0  — 


_ ,z !  P\ _ 

^20.1  ^  P2  "h  ^20,i  ^  Pi 


(13.10) 


are  the  reflection  and  transmission  coefficients  of  the  interface. 


Consider  the  case  that  the  interface  is  a  random  rough  surface  given  by  the  equation: 


z  =  a(r ) 

where  a(r)  is  a  random  variable  with  zero  mean  and  slope. 


(13.11) 


<#  >=0 


<  Vj/ar  >=  0 


(13.12) 

(13.13) 


The  pressure  field  is  required  to  satisfy  the  boundary  conditions: 
pl(z  =  a)  =  p2(z  =  a) 

1  3 p,  (z  =  or)  1  3p2  ( z  =  a ) 
p,  dn  p2  dn 


(13.14) 


(13.15) 


where 


|~  =  V„  =  .  1  (-V  xa»V  +  j-) 

'Jl  +  V1CfV±a  3z 


(13.16) 


is  the  derivative  in  the  direction  of  the  normal  to  the  rough  surface: 

1  ,  „ 
n  =  ...  . - - =  -  ( - , - ,1) 

yjT+Vjx^Vjx  a*  dy 


(13.17) 


Let  us  express  the  pressure  field  in  Regions  One  and  Two  as  the  sum  of  the  mean  field 
and  a  zero  mean  stochastic  field. 
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P,  =<  P\  >  +W1 
Pi  =<Pl  >+^2 


(13.18) 

(13.19) 


Expanding  the  fields  about  the  mean  surface,  one  arrives  at  the  following  boundary 
conditions  on  the  mean  surface  (z  =  0): 


„  .  .  «!  %>i  _  dP:  .  a‘  fe. 

Pl  IT  T"aT~P!  3z  2  9z 


(13.20) 


P. 

J_ 

P2 


3p,  92p,  a  93p,  f7  _  fy 

^-L+a-Ti7-  + - -T-- Vxa«Vp, 

9z  9z2  2  9z 


3p,  32p2  a1  9  p2  ^  ^ 


(13.21) 


Here,  one  has  only  kept  terms  up  to  quadratic  terms  in  the  random  variable.  Taking  the 
ensemble  average  of  the  above  equations,  one  obtains  the  following  equations  for  the  mean  field 
on  the  mean  surface  (z  =  0): 


<p,  >-<p2  >=FX 
1  9  <  p,  >  1  9  <  p2  > 


p,  3z  p2  9 z 


=  F, 


(13.22) 

(13.23) 


The  source  terms  on  the  right  hand  side  are  given  by  the  following  expressions  quadratic 
in  the  random  variable: 


Fx  =-<d 


F2  =  -  <  ai 


3w,  3w2 

<a2 

> 

f~92  <  p, 

>  3 2  <  p2> 

3z  9z 

2 

9z2 

3 z2 

1  32w, 

1  92w2 

<a2  > 

1  9 3  <  p,  > 

1  9 3  <  p2  > 

.Pi  dz2 

Pi  dz2  _ 

P 

2 

_P.  &3 

P2  . 

(13.24) 


<  V±a»  V| 


w,  w2 


Pi 


(13.25) 


Note,  the  source  terms  are  quadratic  in  the  random  variable.  Subtracting  Equations  13.24 
and  13.25  from  Equations  13.20  and  13.21,  one  arrives  at  the  following  boundary  conditions  for 
the  stochastic  field  on  the  mean  surface  (z  =  0). 
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w,  -  w2  =  G, 


(13.26) 


1  3w,  1  3w2  _ 

px  dz  p2  dz 


(13.27) 


The  source  terms  on  the  right  hand  side  of  the  boundary  conditions  are  given  by  the 
following  expressions  linear  in  the  random  variable: 


G,  =-d 


d<  p |  >  3  <  p2  > 

dz  dz 


(13.28) 


G2  =-(A 


1  32 


<P\  > 


L^i 


dz2 


1  3 2  <  jp2  > 
/?2  dz2 


+  V±ccV 


<Pi  > 
Pi 


<Pi> 

Pi 


(13.29) 


In  order  to  solve  Equations  13.22  and  13.23  for  the  mean  field,  one  must  eliminate  the 
stochastic  field  from  the  right-hand  side  of  these  equations.  In  order  to  do  this,  one  must  first 
solve  Equations  13.26  and  13.27  for  the  stochastic  fields  in  terms  of  the  mean  field  and  the 
random  wave  height.  One  may  solve  Equations  13.26  and  13.27  by  taking  the  Fourier  Transform 
of  these  equations.  This  enables  us  to  reduce  the  solution  of  these  equations  into  a  set  of 
algebraic  equations  for  the  spectral  strength  of  these  fields. 

Let  us  introduce  the  following  Fourier  decomposition  of  the  random  wave  height: 

=  (13.30) 


Assume  the  cross  correlation  function  of  this  random  variable  is  homogeneous.  Then  the 
cross  correlation  function  depends  only  upon  the  difference  of  the  two  points,  and  is  of  the 
following  form: 


<  a(r)a(r')  >=<  a2  >  cr(r  -  r')  =<  a2  >  —  jd2^LS(^)e+'^ti: 

2  TV 


(13.31) 


where  <a2  >  is  the  root  mean  square  wave  height  squared,  and  S(£)  is  the  power  spectrum  of 

the  surface.  The  following  ensemble  averages  follow  from  the  above  definition  of  the  power 
spectrum: 

<  a(r)a(rj)  >=<  a2  >  S(jj)e+i^'r  (13.32) 


< a(Z)a(Tl) >=  2k8 (£x -rjL)<a2  >S(jj) 


(13.33) 
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The  Fourier  decompositions  of  the  stochastic  fields  in  Regions  One  and  Two  have  the 


following  form: 

2  K  J 

(13.34) 

(13.35) 

The  functions  G,  and  G2  are  linear  functionals  of  the  random  variable  that  have  a 


Fourier  decomposition  of  the  following  form: 

G,  M = -L-  \d^L  -  u)f ,  tf  .17) 

(13.36) 

G2  (r)  =  ~~t  \d2Zx  \d2Ti±e+iL'pa(%  -  tj)g2  (&  7) 

(13.37) 

The  Fourier  coefficients  g{  (g,rj)  and  g2{^,rj)  have  the  following  form  in  terms  of  the 
Fourier  Transform  of  the  mean  field  and  its  derivatives: 
d  <  p\  >  d  <p2> 


(€>v)  =  - 


dz 


dz 


(7) 


(13.38) 


g2(£,7)  = 


1  d2  <  px  > 
~P\  d? 


1  92  <  P2  > 
Pi  te2 


(7)-7l*(I±-7±) 


<?l> 

Px 


±El  > 

Pi 


(7) 


(13.39) 


Substituting  Equations  13.34  through  13.37  into  Equations  13.26  and  13.27,  one  arrives 
at  the  following  algebraic  equations  for  the  Fourier  coefficients  of  the  stochastic  field: 

(£)  -  w2  (£)  =  ~~  [d2PM%  ~  7)g,  (& 7)  ( 1 3-4°) 

—  (£)  - — w2  (^)  =  ^-  \d lrixa(Z  -  7)g2  (£  7)  ( 1 3 .4 1 ) 

Pi  Pi  2*  J 

Solution  of  the  above  equations  leads  to  the  following  expression  for  the  Fourier 
coefficients  for  the  stochastic  field: 
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w, 


«)-f 

2ft  l£l:  I  P\  +  l%2 Pi 


^Ir. 


s,(fl=-pY«  a(4'n) — 

2ft  ^  +^2,JP2 


Pi 


gi(^v)+gi(^n) 


(13.42) 


(13.43) 


Substituting  the  above  expressions  for  the  stochastic  field  into  Equations  13.24  and  13.25 
for  the  source  terms  of  the  mean  field,  one  obtains  the  following  expressions: 


Pi  P\  j 


+tl,:)g2&V) 


<a  > 
(2n)2 

<oP  >  a2 

2  9z‘ 


p2fi  p 


s(f-n) 


AJPi+'Zi.JPi 


i<P\  >-<P2> ) 


(13.44) 


Fi  = 


<a  > 


(2x)2 

tutu 


\d2£L  j d\ 


S{g-rf) 


*4\.J  Pi  + ‘€2.  J  Pi 


P\  P2 


^l.:+^2.z)gl(^n)  + 


(  £2  Z2  ^ 
?U  »2,r 

P.  A 


g2(£*7) 


+ll#(l±-p±) 

--- L1 

y 

gii^rj) 

l  PIP2 

^Pl  Pi  J 

<  a2  >  93  f  <  /?,  >  <p2>^ 


2  dz3 


A 


Pi  ) 


(13.45) 


To  lowest  order,  the  mean  field  is  given  by  the  unperturbed  solution  given  below: 

<  P1.0  >=  exp[+/£01  •?] {exp[-/£]0  :z] + R0  exp[+/£10,z] }  ( 1 3.46) 

<Pi.o  >=T0exp[+ik01  • r-ik20:z ]  (13.47) 
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k0  x  is  the  horizontal  wave- vector  of  the  incident  field,  and 


is  the  vertical  wavenumber  in  Region  One,  and 


(13.48) 


is  the  vertical  wavenumber  in  Region  Two. 
coefficients  of  the  incident  wave: 


(13.49) 

R0  and  T0  are  the  reflection  and  transmission 


*10  ^20,z  !  Pi 

*io,z !  P\  *20, z  ^  Pi 


70=(1  +  R0)  = 


2*10 .1  /  P\ 


*10,2  ^  P\  *20,z  /  Pi 


(13.50) 

(13.51) 


The  Fourier  coefficients  of  the  unperturbed  solutions  are  given  by  the  following 
expressions: 


<  Ko  >  (€) =  -*o,±){exPH'*io,zz]+*o  exp[+*l0  ;z]}  (13.52) 

<  Pi,o  >  (£)  =  2 a£(£i  ~KJT0  exp [~ik20zz]  (13.53) 


Since  the  mean  field  and  stochastic  field  are  accurate  up  to  terms  quadratic  order  in  the 
random  wave  height,  one  can  substitute  the  unperturbed  solution  for  the  mean  field  in  Equations 
13.38  and  13.39  for  the  stochastic  field  source  terms.  This  substitution  leads  to  the  following 
simplification: 

g,(£,t 7)  =  2k8{xi1  -£0,x)2i(£)  (13.54) 

g2  =  2/r 8  (v x  -  k0L  )g2(£)  (13.55) 

l  (4)  =  *,  o,2  (1  -  )  ~  d  +  *o )  ( 1 3-56> 


g2(£)  =  (l  +  tf0) 


k2 


10,2  k  20,2  r  r  w  1  1>. 

- *0.1  •(£±-^0,1 )(--—) 

P\  Pi  P\  Pi  . 


(13.57) 


Similarly,  one  may  substitute  the  above  results  and  the  unperturbed  solution  into 
Equations  13.44  and  13.45  for  the  source  terms  of  the  mean  field  without  loss  of  accuracy. 
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<cr  > 


exp[/\x*r]|j2^j 


7  Sf„\k°/\n 

S\,z'Pi+*S2,Jp2l  Pi  Pi 


<  /y  >  - 

+— —  expt/V  •r](\  +  R0)(kl -kl) 


(13.58) 


<a  > 


exp[ik0L  •  r]  jd2<*±  — *0,l) — 

J  ^,JP\+lkJp2 


it  t  t2  t2  \ 

-^fL^,+l2,)gl(^)+  &<*) 

P1P2  l  Pi  Pi) 


,S  r  /*'(£., +&,)~/e,.f  1  1  V,*x 

+  Cx*(si-^o,i)  - — - £i(s)  + - g2(£) 

1  P\Pi  \P\  Pi) 


.<0T>  -  -...  *10.*  _.  *20..-  ,, 

^-exp[ifr0JL»r] - (l-/?0) - (1  +  P0) 

2  Pi  Pi 


(13.59) 


In  the  above  equations  one  notes  that  Equations  13.58  and  13.59  may  be  expressed  in  the 
following  form  proportional  to  the  root  mean  square  (rms)  wave  height  and  incident  field 
evaluated  at  the  unperturbed  surface. 


Ft  =<a2  >  p.nc(z  =  0)Fl 


(13.60) 


F2=<a2  >pinc(z  =  0)F2 


(13.61) 


*  «^l»7-7)-«u  +&.»«> 

2n  'Si..- 'Pi  +i$i,z  'Pi  L  Pi  Pi 


+-(!+/?„  Ml-kl) 


(13.62) 
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. j.  f^x  ^-y_x 

2^r  /  P\  +  z*^2,r  ^  Pi 


Mm 

P\Pl 


(ft.,  +<ft..)*l(ft)  + 


'ft!. 


ft  (ft) 

y 


.1 

-I— 

2 


y  7,3 


-•(ft.,+ft;...) 

P\Pi 

e 


r 


k~'°J  (l-R.)-^-(\  +  R0) 


V 


A 


Pi 


*,(£)+ 

\ 


1  1 


A 


VA  A 


&(f) 


The  generalized  reflection  (7?)  and  transmission  (7)  coefficients  are  defined  by  the 
following  equations: 


<A  >=(l  +  7?)p,n,(2  =  °) 

<Pi  >=Tpinc(z  =  0) 

±3<£1>_  =  -  R)p,Jz  =  0) 

P\  OZ  p, 

~^2±=_i“2>'”(z=0) 

A  Pi 

Substituting  the  above  expressions  into  the  boundary  conditions: 

< Pi  >-<P2  >=  Fw  =< a 2  > Fxpinc(z  =  0) 

±3<£l>_±3<^1>  =  f^  =<0?>PiPk(z  =  o) 

p,  oz  p2  dz 

one  arrives  at  the  following  equations  for  the  generalized  reflection  and  transmission 
coefficients: 

(1  +  R)-T  =<a2  >  Fx 

k,  k1 ,  , 

-  /— (1  -R)+ i—T  =<  a1  >F2 

P\  Pi 


(13.63) 

(13.64) 

(13.65) 

(13.66) 

(13.67) 

(13.68) 

(13.69) 

(13.70) 

(13.71) 
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Solving  linear  Equations  13.70  and  13.71,  one  arrives  at  the  following  expressions  for  the 
generalized  reflection  and  transmission  coefficients  for  the  coherent  component: 


R  = 

T  = 


kJ_P\  /  p2 

k\,z  l  P\  +  ^2,Z  I  Pi 

^KJPx 

KJ  A  +k2.J  Pi 


■+<a2  > 

K. 

.Pi 

—  <a 2  > 

.Pi 

F\  —iF2 


J/(*i  J  Pi  +k2,JPl) 
WlJPi  +klrJ  Pi) 


(13.72) 

(13.73) 


The  Fourier  transform  of  the  cross  correlation  function  of  the  scattered  field  is  given  by 
the  following  expectation  value: 


<  w, (£)£,*(£')  >=  2 k8{£l  -£/)  <  a2  >  5(£) 


itiMZVPi+M) 


kz'Pl+kz1  Pi 


(13.74) 


The  cross  correlation  function  of  the  scattered  field  is  given  by  the  following  expression: 

<  w,(r)w,*(r')  >= 


(13.75) 


2  n 


\d2k  <a2>  S{g) 


'Zi.M&lpi  +g2(t) 


kz1  Pi+kz1  Pi 


«Ppfi  •  (T  -  r')  +  i£(z  -  z')] 


Similarly,  the  cross  correlation  of  the  field  scattered  into  Region  Two  is  given  by  the 
following  expression. 


<  w2(r)w2’(r')  >= 

-  jd2{x  <  a*  >  S(f) 


2  K 


P+glik 


kz/p+krJPi 


exp[/|±  •  (r  -  ?) -  i£:(z -  z’)] 


(13.76) 


Suppose  the  spectrum  5(£±)  is  sharply  peaked  around  the  origin,  so  that  one  is  justified 
by  approximating  the  spectrum  as  the  delta  function: 

5(£x)  =  2;r52(£x)  (13.77) 

in  the  above  integrals.  In  this  case  one  arrives  at  the  following  asymptotic  expression  for  the 
coherent  reflection  coefficient  and  transmission  coefficient,  where 
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and 

Ao,  *  =  ^K,i*K,±-k2  (\3.79) 

are  the  vertical  wavenumber  in  Media  One  and  Two  in  the  case  of  a  flat  interface,  respectively. 

A 


^WA  k20,Jp2_+<a2  > 

klQ ,z  !  A  +  k20,z  I  Pi 


v20,z 

Pi 


Ao  iF20 


/(Ao,  ! Pi  "*■  Ao,z ! Pi) 


(13.80) 


T  =  - 


2*10  ,z!P\ 


k\0,z  I  Pi  "h  k20,z  I  Pi 


—  <a> 


MO.z 


A 


F  +  /F 

2  io  T  “20 


Kk\0,z  /  A  k20,z  I  Pi) 


(13.81) 


A>  = 


1 


1 


ikl0,zlPl+ik20,z/Pl 


kiO,zk20,zSlo(  )  i(klO  ,z  +  k20,z)§20 

Pi  P\ 


+  -(1  +  Ro)(Ko.z-kk:) 


(13.82) 


Ao  — 


Ao.z  I  P\  "h  ik20,z  I  Pi 


ikwAo,zg\o(kU),z  +  k20,z)-r—  +  (~L--TL)glO 


P\  Pi  P\  Pi 


) 


2  Pi  pi 


(13.83) 


S\0  ~  Ao,z  0  A) )  ikio,z  (1  +  A  ) 


(13.84) 


k2  k2 

~  /I  ,  D  \  f  10,2  *20,2  * 

£20  —  (1  ■*"  A)(  } 

A  P2 


(13.85) 


A 


Ap,z  Pi  k20,z  {_  Pi 
kiO,z  t  P\  +  Ao.z  f  Pi 


(13.86) 


After  some  tedious  algebra,  the  following  simplifications  maybe  made  to  the  above 
equations: 


R  =  R0(l-2k?0z<a2>) 


(13.87) 


T  =  T0(l  +  ^< a2  >(k]0:-k20;)2) 


(13.88) 
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The  reflection  and  transmission  coefficients  thus  given  agree  with  Equations  57  and  58  in 
Reference  18.  Equation  13.87  also  agrees  with  the  low  frequency  limit  of  the  Kirchhoff 
approximation,  which  gives  an  effective  reflections  coefficient  for  the  mean  field  of  the 
following  form: 

R  =  R0  exp[-2 k^  :  <  a2  >]  ( 1 3.89) 


Having  derived  the  effects  of  a  surface  with  small  wave-height  surface  roughness  on  the 
reflected  and  transmitted  field.  Let  us  reformulate  this  result  into  a  format  amenable  to 
calculation  in  a  waveguide. 


Let  pi0  and  p2Q  represent  the  solution  for  the  pressure  field  in  Media  One  and  Two  in  the 
absence  of  surface  roughness.  They  satisfy  the  boundary  conditions: 


P io  ~  Pv> 

1  _  j_  dp20 

Pi  dz  p2  dz 

at  the  interface  between  Media  One  and  Two,  and  the  differential  equations: 

=  o 

£±ZL  +  (kl-q')Pl,=  o 

dz 


(13.90) 

(13.91) 


(13.92) 

(13.93) 


Let  dp^  and  dp2  denote  the  difference 
%>x  =<  P\  >"PlO 
5p2  —<  p2  >  —  P20 


(13.94) 

(13.95) 


of  the  mean  field  in  Media  One  and  Two  in  the  presence  of  surface  roughness  and  the  flat  surface 
field  in  Media  One  and  Two.  They  satisfy  the  boundary  conditions 


2  ^  /-l  _  .  l  1  dpl0  N .  ^,1  ^  l U  _  ,  1,  ^  ^20 


dpi-dp^Fi  =<a  >(biPl0+b2—:^)=<ai>(biP20+b2 - — 

Pi  dz  A  dz 


) 


71^1  — =f2  =<«2  > (ci Ao + c2 ■7-~7~)  =< 0(2 > (ci P20 + c.  {  dPl° 


Pi  dz  p2  dz 


pi  dz 


pi  dz 


) 


(13.96) 

(13.97) 


13-12 


CSS/TR-02/10 


at  the  interface  between  Media  One  and  Two,  and  satisfy  the  same  differential  equations  as  the 
functions  pl0  and  p20 ,  respectively.  The  source  terms  Fl  and  F2  are  given  by  the  following 

expression: 


F,= 


<ar  > 
In 


s(L  -fcj 


(JuS, « X—  -  — )  -  +  ft.,  )*2  <f  > 

P2  Pi 


<<*  >,,2  ,2 


■K-KM 


10 


(13.98) 


Fi  = 


<Cp  > 

In  '  1  /#u/  Px+i^J  Pi 


Pi  Pi 


(  £2 


(^i,z  +£>,J£i(£)+ 


2  A 


Pi  Pi  j 


+li*di-4x) 


Pi  P2 


*,(£)  + 


A 


Pi) 


&(£) 


V 

/ 


<ar  > 


/j-2  /-2  )  ^  #10 

- \K  10 ,2—K.  20,z - 

2  pt  dz 


(13.99) 


where  one  has  neglected  higher  order  terms.  The  functions  g,  and  g2  are  given  by  the  following 
expressions  in  terms  of  the  functions  pl0  and  p20 : 


J,(f)  =  -(P,  =  -(P,  -  ft)—  (13.100) 

p,  dz  p2  arz 

?,  (O  =  (—  -— )  A.  -  *1  •  (fi  -  ^  X—  -  — )Ao  (13-101) 

Pi  Pi  Pi  Pi 

Substituting  the  above  expressions  into  Equations  13.98  and  13.99,  one  arrives  at  the 
following  expressions  for  the  coefficients  of  the  boundary  conditions  for  the  perturbed  fields: 


k=\K-k L> 


S( fl 


A) 


2k> 


■({,,  +f,-) 


’  A:2  A:2 

Pi  Pi 


•*i»(fi-*iX— -— ) 

Pi  Pi 


(13.102) 
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'  (P.  -P2)2  f  j2  e  5(£±-*i 


2#  P\Pi 


S\.J  P\+ZlJ  Pi 


(13.103) 


=~b F~/gV/X/o  ~^X---)i 

^  Sl,z  '  P\  +b2,r  '  Pi  P\  Pi  P\  Pi 


P\  Pi 


2^  /?,  /)2  G\,z'P\+bl,;'Pl  P\  Pi  P\  Pi 

(13.104) 


Ci  -  2(*‘.=  *^)  i  ^  *(^  _*i)> 

(13.105) 

Let  us  make  the  simplifying  assumption,  that  the  surface  roughness  spectrum  is  isotropic. 
In  this  case  one  can  reduce  the  above  expressions  into  the  following  relations: 


2  ,2  ~ 


b\  =  ^(ki  -&?)-(- —  —)  \q'dq' 


S(q'~q ) 


A  Pi  0  £.z  fpi  +tij  A 


-(6.,+£.z) 


(13.106) 


42=+/M!lwJ^ 

P\Pi  0  q  ^/A+^/A  "  '■ 

c-=-44>2>^j^tV 

A  A  0  £i,z  I  P\+  £2 ,z  I  Pi 


(13.107) 


(13.108) 


(13.109) 


Here  <7  =|  kx  |  is  the  magnitude  of  the  horizontal  wavenumber.  The  above  expressions  can  be 
substituted  into  Equations  13.96  and  13.97  for  the  boundary  conditions  of  the  perturbed  field. 
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14.  EFFECT  OF  SMALL  WAVE-HEIGHT  SURFACE  ROUGHNESS  ON 
PROPAGATION  OF  SOUND  IN  THE  OCEAN 


This  section  describes  the  effects  of  surface  and  bottom  roughness  on  propagation  of 
sound  in  an  otherwise  horizontally  stratified  waveguide. 


z  =  aT(x,y) 


(14.1) 


z  =  d+aB(x,y) 


(14.2) 


denote  a  realization  of  the  randomly  rough  surface  and  bottom  of  the  waveguide,  respectively. 
The  random  fields  aT  and  aB  are  homogeneous,  zero  mean  random  fields  with  cross  correlation 
functions  of  the  form: 


<  aT  ( r)aT (/*')  >=<  a2  ><JT(r-r')  =  <C‘T  >  jd2^±Sr (£x ) exp [+/£x  •  (r  -  r ')] 


(14.3) 


< aB ( r)aB (r‘)  >=< a2  ><JB{r-r')  =  <Ct-B-  — 

L71 


Jj2^x5s(^x)exp[+i|x  -F)]  (14.4) 


where  ST  (J;L )  and  SB(^±)  are  the  spectral  density  of  the  rms  wave  height  of  the  surface  and 
bottom,  respectively. 

In  the  case  of  a  flat  surface  and  bottom,  the  normal  mode  contribution  to  the  propagation 
of  sound  in  the  waveguide  is  of  the  following  form: 


/W  =-T-J— -Z# o)(q„r)F(z:q„)F(zs  :q„ ) 

4/>(z,)  „ 


(14.5) 


where  F(z  :qn)  is  the  depth  function  of  the  n’th  normal  mode.  The  depth  function  is  a  solution 
of  the  differential  equation: 


P(z)-T~7 Tjf(2 :  ?„)+«*)’ -gl)F(z :  ?,)  =  0 
dz  pyz)  dz 


(14.6) 
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subject  to  the  following  boundary  conditions: 


fr  (q.)F(i  =  0 :  q, ) + SlMlEk^llA  =  o 

p(  0)  dz 

(14.7) 

p(d)  dz 

(14.8) 

(14.9) 

gr(9)  =  0  +  *r.o(4)) 

(14.10) 

P(d) 

(14.11) 

gs(?)  =  (l  +  ^,o(9)) 

(14.12) 

Here, 

h(z :  q)  = +i-Jq*  -£2(z) 
is  the  vertical  wavenumber,  and 
*r,o(4)  =  -l 

R  (q)  -  W  i q) i pjd)  -  hg  (d :  q)  /  (cf) 


(14.13) 

(14.14) 

(14.15) 


h(d :  q)  /  /?(<*)  +  hB(d:q)/  pB  {d) 
are  the  reflection  coefficients  of  the  surface  and  bottom  in  the  case  of  a  flat  surface  and  bottom. 
Assume  the  sediment  is  modeled  as  a  fluid  layer,  where 

hB(z:q)=+iy[q^kJ(^)  (14.16) 


is  the  vertical  wavenumber  in  the  sediment,  kB(z)  is  the  wavenumber  in  the  sediment,  and 
pB(z)  is  the  density  of  the  sediment. 


The  effect  of  small  wave-height  surface  roughness  on  the  coherent  (mean)  component  of 
the  normal  modes  can  be  accounted  for  by  substituting  Equation  13.80  for  the  effective  reflection 
coefficient  of  the  rough  interface  into  Equation  14.9  through  14.12.  In  evaluating  the  parameters 
the  quantities  in  Section  13  are  defined  by  the  following  limiting  value  of  the  parameters  on  the 
two  sides  of  the  interface  at  z=d: 
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Pi  =  P(z  =  d+£) 

(14.17) 

p2=limE^p(z  =  d-£) 

(14.18) 

k{  -  lim£_,0  k{z  =  d  +  £) 

(14.19) 

k2  =  lim£^0  k(z  =  d-£) 

(14.20) 

An  alternative  method  of  accounting  for  small  wave-height  roughness  is  to  treat  the 
effects  of  surface  roughness  as  a  perturbation  of  the  normal  mode  equations.  Let  F(z :  qn ) 
denote  the  depth  functions  of  the  unperturbed  waveguide  with  flat  boundaries. 

These  functions  are  solutions  of  the  differential  equation: 

— F{-\ '  g)  +  (k2(z)  -q2)F(z  :q)  =  0  (14.21) 

dz 

and  satisfy  the  following  normalization  condition. 


+00 

jdzF(z:q)F(z:q)/ p(z)  =  l  (14.22) 

o 

Here,  one  must  make  the  assumption  the  density  is  piece-wise  constant.  The  effect  of  surface 
roughness  is  to  perturb  the  above  eigenvalues  and  eigenfunctions. 

Let  F{z  :qn)  +  SF(z :  qn )  denote  the  perturbed  depth  function  and  qn+8qn  the  perturbed 

eigenvalue  due  to  roughness  of  the  interface  at  z=d.  The  perturbed  eigenfunction  satisfies  the 
following  boundary  condition  at  z=d: 


limf_>0(^F(J  +  £ q)  -  SF(d -e:q))  =<  a 2  >  (bxF(z  =  d\q)  +  b2 


1  dF(z  =  d:q) 


dz 


■) 


lim£_0( 


1 


p(d  +  e)  dz 


d  SF(d  +  e:q)-  1 


p{d  -  £ )  dz 


SFV-S-.q))  = 


<a 2  >(c]F(z  =  d  :q)  +  c: 


1  dF(z  =  d:q) 


dz 


) 


(14.23) 


(14.24) 


The  coefficients  on  the  right  hand  side  of  the  boundary  conditions  are  given  by  either 
Equations  13.102  through  13.105  for  a  general  surface  roughness  spectrum,  or  by  Equations 
13.106  through  13.109  for  an  isotropic  surface  roughness  spectrum. 
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The  perturbed  depth  function  satisfies  the  differential  equation: 


d\F\z :  q)  +  5F{z :  q) 
d? 


+  (k2(z)-q2  —8q2)(F(z :  q)+SF(z  :q))  =  0 


(14.25) 


where  Sq2  is  the  perturbation  of  the  horizontal  wavenumber  squared.  Keeping  only  first  order 
terms,  one  arrives  at  the  differential  equation: 

dl5F£2:q)  +  tt\z)  -q‘)SF(z :  q)  =  Sq1F(z:q)  (14.26) 

relating  the  perturbed  depth  functions  and  the  perturbed  eigenvalue.  Here  the  perturbed 
eigenvalue  acts  as  a  source  term  for  the  perturbed  depth  function.  This  result  is  typical  of  first 
order  perturbation  theory. 


Combining  Equations  14.21  and  14.26  one  arrives  at  the  following  result  for  the  depth 
derivative  of  the  Wronskian  of  the  perturbed  and  unperturbed  depth  functions: 

F(z :  q)-——^-  —  —  SF(z :  =  Sq2F(z :  q)F(z :  q)  (14.27) 

clz  dz 


Dividing  the  above  product  by  the  density  and  integrating  the  resulting  product  over  the 
entire  waveguide,  one  arrives  at  the  following  relation. 


Sa2  -  8a 2 "f dz  F(Z :  *)F(Z :  &  -  "f A  F(Z :  g)<S?r,,(z :  g)  ~  F"(z :  g)<SF(z : 
oJ  PM  l  P(z) 


(14.28) 


Here,  one  has  made  use  of  the  orthogonality  condition  of  the  unperturbed  depth  functions. 
Assuming  the  density  is  piece-wise  constant,  the  above  integrand  is  a  total  derivative  in  each 
layer  of  the  waveguide  except  at  the  interface  z=d,  where  the  perturbed  depth  function  is 
discontinuous.  In  this  case,  the  above  integral  is  equal  to  the  following  limit: 

Sq2  =  -lim£_,0(F(z :  q)SF'{z :  q)-F\z :  q)SF(z :  q))l  p(z)  \dzted_e  (14.29) 

Using  the  boundary  conditions  given  by  Equations  14.23  and  14.24  for  the  discontinuity 
of  the  perturbed  field  at  the  boundary,  one  arrives  at  the  following  expression  for  the  perturbed 
eigenvalue: 

Sq2  =  - < a2  >  {c, F(d : q)2  +  (c2  - 6, )F( d : q)F\d  :q)/ p-b2(F'(d :q)/ p)2}  (1 4.30) 
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To  first  order  in  perturbation  theory,  the  perturbed  depth  functions  are  equal  to  the 
unperturbed  depth  functions  evaluated  at  the  perturbed  eigenvalue  q  +  Sq .  In  particular,  the 
perturbed  depth  function  8F{z :  q)  is  given  by  the  difference 

SF(z:q)  =  F(z:q  +  Sq)-F(z:q)  (14.30) 

Note,  in  the  above  first  order  perturbation  theory,  the  effect  of  multiple  rough  surfaces  is 
additive  to  first  order.  Thus,  one  can  represent  the  effect  of  a  rough  surface  and  bottom  by  the 
sum  of  the  perturbations  of  each  surface  on  the  unperturbed  eigenvalues. 

In  the  case  of  a  sound  soft  surface  at  z  =  0,  the  perturbation  of  the  eigenvalue  is  given  by 
the  following  relationship: 

Sq2  =  +  <a2  >b2(F'(z  =  0:q)/p(0))2  (14.31) 

The  coefficient  b2  is  given  by  the  following  limit  p2  — >  0  of  Equation  13.107  in  the  case 
of  an  isotropic  surface  roughness: 


OO 

=  +ipi  \q'dq'S(q'-q%,z  ( 14.32) 

0 

In  the  case  of  a  rigid  rough  surface  at  z=D,  the  perturbed  eigenvalue  is  given  by  the 
following  expression: 

Sq2=-<a2>cx(F(z  =  0-.q))2  (14.33) 

The  coefficient  c,  is  given  by  the  following  limit  p,  — >  +®°  of  Equation  13.108  in  the 
case  of  an  isotropic  rough  surface: 


Sjq'-q) 

Su 


(14.34) 
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15.  ASYMPTOTIC  EVALUATION  OF  THE  TIME  DOMAIN  SOLUTION 


This  section  presents  an  asymptotic  analysis  of  the  propagation  of  a  Gaussian  pulse  in  a 
waveguide  in  the  time  domain.  This  analysis  is  based  on  a  saddle  point  evaluation  of  the 
asymptotic  evaluation  of  the  normal  mode  expression  for  the  Green’s  Function. 

The  Green’s  Function  for  the  propagation  of  a  monochromatic  signal  in  a 
waveguide  is  of  the  following  form. 


G(r,zs,z ) 


i 

4  p(zs) 


(15.1) 


Substituting  the  asymptotic  expression 


H"\qr)~\—ei(qr-*,4)  (15.2) 

\nqr 

for  the  Hankel  Function  into  Equation  15.1,  one  obtains  the  following  asymptotic  expansion  of 
the  Green’s  Function: 


G(r,zs,z)  ~ 


+i{qnr-7t!4) 


4 p(zs)  „  \mqnr 


F(zs:qn)F(z:qn) 


(15.3) 


Propagation  of  a  pulse  is  described  by  the  Fourier  Transform: 


+oo 

G(r,zs,z :  t )  =  jS(co)G(r,zs,z  :aj)e~,M da) 


(15.4) 


where  S{co)  is  the  spectrum  of  the  incident  signal.  In  the  following  one  may  assume  a  spectrum 
of  the  following  form: 

S((o)  =  exp(-  )  (15-5) 

2A  a 

where  6)0  is  the  center  frequency,  and  A  to  is  the  bandwidth  of  the  pulse.  The  incident  signal  is 
of  the  form 

+oo  __ 

s(t)=  ^S(co)e~“0'  da)  =  yf2aA(oexp(-t2A6)2  l2-icoQt)  (15.6) 
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where 


r 


1 

A  co 


(15.7) 


represents  the  pulse  length  of  the  incident  signal. 


Substituting  Equation  15.3  and  15.5  into  Equation  15.4,  one  obtains  the  following 
expression  for  the  propagation  of  a  Gaussian  pulse  in  a  waveguide: 


G(r,zs,z:t)  = 

.  +oo 

7-7—-Z  J^POW -6X-7r/4)-(6)-6)0)2  / 2Aco2)F(zs  :  q„ )F(z :  qn )d(0 

**P\ZS)  n  —00 


(15.8) 


One  must  now  assume  that  this  integral  is  dominated  by  the  exponential  term  and  the 
depth  functions  inside  the  integral  are  approximately  equal  to  their  value  at  the  center  frequency 
of  the  pulse.  Under  this  assumption,  evaluation  of  Equation  15.8  reduces  to  the  evaluation  of 
integrals  of  the  following  form: 


/„  (r,  zs,z:t)=  jexp(z  (qnr  -  cot  -  n /  4)  -  (eo  -  co0  )2  /  2  Ary2) dco 


(15.9) 


where  the  horizontal  wavenumber  qn(a))  is  an  implicit  function  of  frequency.  Expanding  the 

argument  of  the  exponential  function  about  the  center  frequency,  one  obtains  the  following 
expression  for  the  integral: 


-too 

/„ (r , z5 , z : 0  - exp(/(<7„ (a)0)r - o)0t - /r / 4)  Jexp {-a2{co-o)0f -i{(o-co0){t-tn))dco  (15.10) 


where  the  reciprocal  of  the  group  velocity  of  the  mode  is: 


q\  (®)  = 


dco 


the  rate  of  change  of  the  reciprocal  of  the  group  velocity  is 


*"„(<»)  = 


d\Scd) 

dco 2 


the  travel  time  for  the  mode  is 


t*  =rq\  (<»o) 


(15.11) 


(15.12) 


(15.13) 
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and  the  complex  width  squared  of  the  effective  Gaussian  is 
1 


a2  = 


2A*y 


2  2 


(15.14) 


Equation  15.10  can  now  be  computed  by  completing  the  square  in  the  argument  of  the 
exponential  function  to  recast  the  integral  into  a  simple  Gaussian  of  the  form: 


jexp(-#2x2)<s!r 


(15.15) 


The  resulting  integral  is  given  by  the  following  expression: 

I„(r,zsiz:t)~—  exp(i(qn(a)0)r-a)0t-n/4)-(t-t„)2 /4a2)  (15.16) 

a 

Substituting  Equation  15.16  into  Equation  15.8,  one  obtains  the  following  asymptotic 
expression  for  the  time  domain  representation  of  a  Gaussian  pulse  propagating  in  a  waveguide: 

G(r,zs,z:t)  = 

(15.17) 

- 4 :  <ln)F(z :  tfjexp (i(q„r -(0Qt-nlA)-{t- tnf  / 4a2) 

4p(zs)  q„r7C  a 

where  the  normal  modes  are  evaluated  at  the  center  frequency. 
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16.  REPRESENTATIONS  OF  CYLINDRICAL  WAVE  FUNCTIONS 


This  section  focuses  on  different  representations  of  cylindrical  wave  functions. 
The  regular  cylindrical  wave  functions  are  given  by  the  following  expressions: 


Re*£  (r :  a)  =  (kp sin (a))e+iM lcos(m<p)<7  el 


8;r 


\sm.{m(p)(J  =  o , 


(16.1) 


(16.2) 


Here,  the  function  e(m)  is  defined  by  the  following  relationship: 


e(m)  = 


1, m  =  0] 

2, m*  0| 


(16.3) 


The  irregular  cylindrical  wave  functions  corresponding  to  outgoing  waves  are  given  by 
the  following  expressions: 


X%(r:cc)  =  , 


-Hnm(kpsin(a)y  . 

&7U  [sm(m(p)<T  -  o  j 


(16.4) 


_  l£(m)  u(D  /t/,c;r>r/v^fl_,fcC0S(af)/cos^W3^<7 


ir '<*)=- 


8# 


-//'  m(£psin(tfr))e~ 


[sin(m9>)<7  =  o 


(16.5) 


The  variable  a  represents  the  angle  of  the  wave  with  respect  to  the  z-axis.  The  horizontal 
wavenumber  q,  and  vertical  wavenumber  h  are  given  by  the  following  expressions  in  terms  of 
this  angle: 

q  =  ksm(a)  (16.6) 

h  =  kcos(a )  (16.7) 

The  downward-going  cylindrical  wave  functions  and  upward-going  wave  functions  are 
related  by  the  following  relationship: 


Re  (r :  a)  =  Re  z™  (r :  n  -  a) 


(16.8) 
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The  following  expansion 

exp[;xcos(or)]  =  £  inJn(x)e+ina=J0(x)  +  2j^inJn(x)cos(na) 

n~~ 00  n=l 

follows  by  substituting  the  expression 
t  =  z'exp(7a] 

into  the  generating  function  for  the  regular  Bessel  Functions: 
exp[*(f -/-')/ 2]  = 


(16.9) 


(16.10) 


(16.11) 


Equation  16.8  serves  as  the  basis  for  the  following  expansion  of  a  plane  wave  in  terms  of 
the  regular  cylindrical  wave  functions: 


— 

exp[±/£  •  r]  =  Y,  e(m)(±i)nJm  (kp  sin(tf))  exp[±ifc  cos(a)]  cos{m{cp  -  /?)) 

m=0 


=  4*E£(±0"Re;r“(r:a) 

m=0  C 


cos(m/3),cr  =  e 
sin(m/3),a  =  o 


where  the  components  of  the  wave  vector  are  given  by  the  following  expressions: 
kx  =  ksin(a)cos(fl) 

ky  =  £sin(#)sin(/?) 


k,  =kcos(a) 


(16.12) 


(16.13) 

(16.14) 

(16.15) 


Using  the  orthogonality  of  the  trigonometric  functions  and  the  above  expansion  of  the 
plane  wave  in  terms  of  regular  cylindrical  wave  functions,  one  arrives  at  the  following  plane 
wave  expansion  of  the  regular  cylindrical  wave  functions: 


**Z™{r:a)  = 


1 

4  K(±i)m 


cos(/m/?),  a  =  e 
s‘m(mft),CF  =  o 


exp[±z'£  •  r]d{3 


(16.16) 


The  above  integral  representation  of  the  regular  cylindrical  functions  may  also  be  inferred 
from  the  following  integral  representation  of  the  regular  cylindrical  Bessel  Function: 
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Jn (x)  =  -^^^exp[±ixcos(fi)]cos(nfi)dfi  =  2^+jy  Jexp[±zxcos(/?)]cos(«/?)J/? 


^■(±0 

=  2x(+iy  Jexp [±ixcos(/3)±infi]dp 


Making  the  following  substitution: 
k*r  =kp  sin(or)  cos(/?  -<p)  +  kz  cos  (a) 

of  the  dot  product  in  Equation  16.16,  one  arrives  at  the  following  integral: 
R  eZ%(r:cc)  = 


— -j-—  j/M|C0S^^  e|exp[±/(k/?sin(a)cos(y5  -  (p)  +  kzcos(a)]dJ3 
An(+i)  o'  In  [sin(m/?),<7  =  o) 

which  is  of  the  form  of  Equation  16.17,  and  leads  to  the  result  listed  in  Equation  16.16. 
A  similar  plane  wave  expansion  of  the  outgoing  cylindrical  waves: 

Zi fr : «)  — - T-rf * ?]  Ii«  lcos(mf  '  ■  eU 

™  2 mm  1  2k  [sin(m/?),<7  =  oj  H 

can  be  obtained  from  the  following  integral  representation  of  the  Hankel  Functions: 


(16.17) 


(16.18) 


(16.19) 


(16.20) 


1  tr/2+i°° 

H(nu  (x)  =  Jexp[+/xcos(/?)  -  in/3~\dp 


K  /  2-/oo 


(16.21) 


The  expansions: 

G(r, r’)  =  exp [ik  \  r  - r'|] IAk  \  r -  r'\=  2/k]T  ^m(a)da Re xlm (r :  a)  Re  (r':  a)  ( 1 6.22) 

am  ct 

G(r, r')  =  exp[zk  |  r  -  r’|]  /4k  \  r -  r'\=  2 ik^  Jsin(a)</fla£,  (n,  :  a) RexL  (r<  :a)  0  6.23) 

am  c 


of  the  free  field  Green’s  Function  follow  from  the  integral  representation  of  the  free  field  Green’s 
Function.  The  contours  C  and  C+  are  the  following  contours: 
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C+  =[0,^/2]U[^/2,^/2-/co) 

(16.24) 

C_  =  (n  /  2  +  i°o,nl  2]  U  [nl2,n] 

(16.25) 

C  =  (+/°°,0]  U  [0,zr]  U  \n,it  -  /oo) 

(16.26) 

The  contour  C+  is  used  in  the  case  (z  -  z')  >  0  and  the  contour  C_  is  used  in  the  case 
(z  -  z')  <  0  in  Equation  1 6.22.  Conversely,  the  contour  C_  is  used  in  the  case  (z  -  z')  >  0 ,  and 
the  contour  C+  is  used  in  the  case  (z  -  z')  <  0  in  Equation  1 6.23. 

A  derivation  of  Equations  16.22  and  16.23  follows.  Start  with  the  following  integral 
representation  of  the  free  field  Green’s  Function: 

C?(r,r*)=  J- ^r^p[+ip»(r-r'Mp2 -k2 -iO)  (16.27) 

J(2n) 

Introduce  the  variable 
M  =  4p2:~k2 

and  express  the  denominator  of  Equation  16.27  as  the  following  product: 

(p2  ~k2  —iO)  =  (py  +  i^p]+M2  - iO)(py  -  iJp2x+M2  - /0)  ( 1 6.28) 

Equation  16.27  takes  on  the  following  form: 

G(r,r')  = 

J  +°°  +00  +oo 

— — j-  fdp:  exp[+(p.z"]  \dpx  exp[+i>,x"]  \dpy  exp [+ipyy"]/(py  +  /-/ p]+M2  -  iO)(py  -  i^p]  +  M2 

V  /  — oo  — oo  — OO 

(16.29) 

The  integral  over  py  can  be  performed  by  the  method  of  residues,  where  one  either 

closes  the  contour  in  the  upper  half  or  lower  half  plane,  depending  upon  whether  the  variable 
y,=  (y~y')  is  positive  or  negative.  Performing  the  integration  over  the  variable  pv  one  obtains 
the  following  expression  for  Equation  16.29: 

G(r,r')  = 

2  j  +°°  +<*> 

r±2i  \^Pz  exp[+(P:z"]  jdpx  exp[+ipxx”]exp[+y"Jp2  +  M2  ]/jp2  +  M2 


(16.30) 
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where  the  upper  sign  is  chosen  in  the  case  /' =  (y  -  y' )  >  0  is  positive,  and  the  lower  sign  is 
chosen  in  the  case  it  is  negative.  Introducing  the  following  transformation: 

px=iMsin(P)  (16.31) 


in  Equation  16.30,  one  arrives  at  the  following  contour  integral: 
G(r,r')  = 

.  +oo  n  /  2 +A» 

— l—  | dp ,  exp[+/p7z"]  ^d/3 exp[-x"A/sin(/?)  +  y"M cos( j3)\ 

jr/2-/«> 


(16.32) 


Expressing  the  argument  of  the  exponential  in  polar  coordinates: 
x"Msin(/?)±/'Mcos(/?)  =  Mp"cos(p±<p )  (16.33) 

Equation  16.32  takes  on  the  following  form: 

+oo  7t!2+i<x> 

G(r,r')  =  — —  jdpz  exp[+zp2z"]  \dP exp[-p" M cos(/? ± cpj\  (16.34) 


The  integral  over  the  variable  p  is  equal  to  the  following  expression: 

1tl2+i<x>  +°o 

J<7/? exp[-Mp"cos(/?  ±  (p)}  =  — /  jd/ l exp[-Mp"cosh(/l)  =  kH^  {iMp")]  ( 1 6.35) 

K  /  2-/oo  -o© 

Substituting  the  above  expression  into  Equation  16.34,  one  arrives  at  the  following  expression 
for  the  free  field  Green’s  Function: 

G(r,r')  =  f dp,  exp [+ip,zn]H^)(i^f^kIp")  (16.36) 

Introducing  the  following  coordinate  transformation: 
p,=kcos(a)  (16.37) 

iM -iyjpl  ~k2  =k sin(cr)  (16.38) 

into  Equation  16.36,  one  arrives  at  the  following  integral  representation  of  the  free  field  Green’s 
Function.: 


G{r,r')  =  —  |<^orsin(a)exp[/fe"cos(a)]//ol)(/:/9"sin(a)) 


8 k  ,J 

+/=» 


(16.39) 
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Replacing  the  integral  over  the  complex  angle  by  an  integration  over  the  horizontal 
wavenumber: 

?  =  £sin(a)  (16.40) 

one  arrives  at  the  following  version  of  Equation  1 6.39: 

G(r’r')  =  ijr  «f  (16.41) 

where  the  integral  is  around  the  branch  cut  of  the  vertical  wave  number: 

h(q)  =  i^jq2 -k2  =  kcos(a)  (16.42) 

that  is,  the  integral  is  along  the  following  contour: 

(+i° °  -e,0-  ie]  U  [0  -  is,  k  -  ie]  U[£  +  ie, 0 + is]  U  [0 + ie,  i°° +£]  (1 6.43) 

where  £  is  an  infinitesimal  displacement. 

One  can  deform  the  contour  around  the  branch  cut  of  the  vertical  wave  number  into  one 
along  the  real  axis,  provided  one  chooses  the  sign  of  the  vertical  wave  number  such  that  the 
imaginary  part  of  the  product  Im(fo")  >  0  is  positive.  In  this  case,  Equation  16.41  takes  on  the 
following  form: 

.1  +OO-/0  I 

=  f  P[+ih\z"\]Ho)(<iP'')  (16.44) 

— OO-/0 


Using  the  relationships: 

(16.45) 

+  to)  (16,46) 

one  can  express  Equation  16.44  into  the  following  integral  over  the  regular  Bessel  Function: 

=  ^  1  ^exp[+/7i|z"|]y0(^)  (16.47) 

W  0-/0  n 

Equations  16.44  and  16.47  serve  as  a  basis  for  generating  an  integral  representation  of  the 
free  field  Green’s  Function  in  terms  of  the  cylindrical  wave  functions. 
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Substituting  the  expansion: 

H<"  {qp")  =  fje{m)H^(qP>  )Jm  (qP< )  cos(m(p  -  (p'))  (1 6.48) 

m= 0 

into  Equation  16.44,  one  arrives  at  the  following  expression  for  the  free  field  Green’s  Function: 

G(r,r')  =  iYp^cos{m((p-(p'))  \H"\qP> )Jm (qP< ) exp[+i/i \z-z'\]qdq/h  (16.49) 

^  8/r  _i 

Equation  16.49  can  be  expressed  in  terms  of  the  outgoing  and  regular  cylindrical  wave 
function: 


G(r,r')  =  /£ 

am 


(r,  :^r)Re^)(r<  :q) 


(16.50) 


where  r>-r  and  r<  =  F  if  P>  P\  and  r>  =  F  and  r<-r  otherwise.  The  upper  sign  is  chosen 
if  the  quantity  (z>-z<)>  0  is  positive,  and  the  lower  sign  is  chosen  if  it  is  negative.  Converting 
Equation  16.50  into  an  integral  over  the  complex  angle  a: 

q  =  ksin(a)  (16.51) 

one  arrives  at  the  following  expression  for  the  free  field  Green’s  Function: 

G(r,  r')  =  ikY,  \x™  0> : «)  Re  X™  (r<  ■ «)  sin(a)da  ( 1 6.52) 

am  c 


where  C  is  the  contour  defined  by  Equation  16.26. 

Substituting  the  expansion: 

JQ(qP")  =  (qP')cos(m(<p  -  q>'))  (1 6.53) 

m= 0 

into  Equation  1 6.47,  one  arrives  at  the  following  expression: 

G(r,r')  =  2iY~~^cos(m(<P ~ <P'))  U m(qp)J  MP')e*P[+ih \z-z'S\qdqlh  (16.54) 

Expressing  Equation  16.54  in  terms  of  the  regular  cylindrical  wave  functions,  one  arrives 
at  the  following  expression  for  the  free  field  Green’s  Function  in  terms  of  the  cylindrical  wave 
functions: 
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+oo 

G(r,r')  =  2/ £  jRej&fr :  ^)Re^(r’: 9) 

0 


qdq_ 

h 


(16.55) 


The  upper  sign  is  used  in  Equation  16.55  if  the  quantity  ( z-z ')  >  0  is  positive,  and  the 
lower  sign  is  used  if  it  is  negative.  Transforming  the  integral  in  Equation  16.55  into  one  over  the 
complex  angle  a,  one  obtains  the  expression: 

G(r,r’)  =  lik^  jRe^(r :  a)Rez£(r':a)sm(a)da  (16.56) 

om  C± 

where  the  contour  C+  is  used  if  the  quantity  (z-z')  >  0  is  positive,  and  the  contour  C_ 
is  chosen  if  the  quantity  is  negative. 
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17.  REPRESENTATIONS  OF  SPHERICAL  WAVE  FUNCTIONS 


This  section  describes  various  representations  of  the  spherical  wave  functions. 

The  regular  and  outgoing  spherical  wave  functions  are  defined  by  the  following 
expressions: 


g("0(2/+l)  (;-'”)!p.(cos(g))|cos(m^’).<7  =  el 

4n  (/  +  /«)!  Isin(m^’),<7  =o\ 


(17.1) 


V„M  =  hm,(.kr\ 


eJjnM + 1)  (/^j!  =  ‘\ 

Ak  ( l  +  m)\  [sin(m^),cr  -o\ 


(17.2) 


where  j,{kr )  and  h\X)(Jcr)  are  the  regular  spherical  Bessel  Function  and  Hankel  Function  of  the 
first  kind,  respectively.  The  functions: 


2  \m/2 


1)' 


dx 


i'n 


dx 


(17.3) 


are  the  real  associated  Legendre  Functions  of  the  first  kind.  These  functions  can  be  extended  to 
the  complex  plane  in  two  ways.  The  first  method  is  to  define  the  function 


/>m(cos(0))  =  (-l)msin(0)w 


im 

■^rPl(x  =  cos(9)) 


(17.4) 


in  terms  of  a  complex  angle  0 .  The  second  method  is  to  define  the  associated  Legendre 
polynomial  as  the  analytic  function 


pr(Z)=(z2- 1) 


-»  d"  !}(*)-' 'z!-'r 


im+l 


dx 


2'1\  dz 


m+l 


(z2-  iy 


(17.5) 


of  the  variable  z,  where  the  function  P,m  (z)  is  related  to  the  real  associated  Legendre  Functions 
on  the  branch  cut  [-1,+1]  by  the  following  expressions: 


um„0pr(x+i£)=(-i)mpr(x) 

(17.6) 

lim£_,0  P,m{x-ie)  =  (+/)m  P,m  (x) 

(17.6) 
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In  the  following  one  may  use  the  extension  of  the  spherical  wave  functions  by  the 
extension  of  the  associated  Legendre  polynomials  defined  by  Equation  17.4  by  the  notation 
Re^,  and  y/^ .  Similarly,  one  may  use  the  notation: 


Re  ft*  O') =.//(*»■) 


g(m)(2/  +  l)  (/-m)! 
\  4k  (/  +  w)! 


p,m(  corny 


co  s{m<p),o  =  e 
sin  {m(p),cr  =  o 


(17.7) 


e)){cQS('"g),g  =  *1 

4#  (/  +  /«)!  1  [sin(w9>),cr  =  oJ 


(17.8) 


to  denote  the  extension  of  the  spherical  wave  functions  to  complex  variables  defined  by 
Equation  17.5.  The  two  representations  are  related  by  the  following  relationships  along  the 
branch  cut  [-1,+1]  of  the  associated  Legendre  polynomials: 


lim^0  W^ix  +  ie)  =  (17.9) 

lim^o  Ymd  (*  ~ i£)  =  (+0"  Von,,  (*)  (17.10) 


The  complex  associated  Legendre  polynomials  satisfy  the  following  addition  theorem: 

P,(cos(y))  =  £e(m)  P,m  (cos(t?))  P"  (cos(  t?’  ))cos(m(<p-(p' )) 

m=o  (l  +  m)l 

+i  (l—  V  0^* 

=  X^(m)(-l)m^— - ^j^”(cos(t?))^'"(cos(t?,))cos(/M(^-^)) 

m-0  (l  +  my. 


The  complex  angle  y  is  defined  by  the  following  expression: 
cos(y)  =  cos(t?)  cos(i?')  +  sin(r?)sin(tf')cos(^>  -  <p')  ( 1 7. 1 2) 

The  associated  Legendre  polynomials  satisfy  the  following  orthogonality  relationships: 


\pr{x)p?{X)dx= 


2  (l  +  m)ly 
2l  +  \(l-m)\  1 


f P,m  (x  ±  iO)P.r  (x  ±  i0)dx  =  (-l)m  — 2  (/  +  w)V 
i  2/  +  1  (l-m)\  ‘ 


(17.13) 

(17.14) 


\pr{x)Pr\x) 


-l 


1  (l  +  m)\ 
m  ( l-m)\ 


S”',m±0 


(17.15) 
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Jy(^±<0)r'(x±/0)7^?=(  'r  0 

_i  VI  — JC2  m 

The  spherical  harmonics  are  defined  by  the  following  expressions: 

=  tomW+l)  (l-mjl  Pr(c0s(j™(^  =  °\ 

V  \  Ak  (l  +  m)\  '  \sm(mtp),a  =  oj 

~  1 
V  4#  (l+m)\  [sin(w^),<T  =  oJ 

They  satisfy  the  following  orthogonality  conditions: 


2*  +JT 


0  -*■ 

2*  +* 


0  -;r 


Using  the  following  identity  from  Reference  9  evaluated  at  v  =  1/2 , 


*=o 


exp[/xcos(a)]  =  r(v)(v/2rX(v'  +  ^'i4+*WCr(cos(a)) 
and  the  identities: 

J i  (x) =  6+1/2 

cr\x)=pk(x) 


(17.16) 


(17.17) 


(17.18) 


(17.19) 


(17.20) 


(17.21) 


(17.22) 


(17.23) 


one  arrives  at  the  following  expansion  of  a  plane  wave  in  terms  of  regular  spherical  Bessel 
Functions: 


exp[+/&  •  r]  =  exp[//7-cos(y)]  =  ^(2/  + 1  )i'j,  (kr)P,(cos(y)) 


(17.24) 


1=0 


The  angle  y  is  the  angle  between  the  two  vectors  k  and  r  .  Let  {$,(p)  and  (a,  p) 
be  the  polar  coordinates  of  the  direction  of  the  position  and  wavenumber  vectors,  respectively. 
Then  the  angle  y  is  defined  by  the  following  relationship: 

cos(y)  =  cos(tt)  cos(z9)  +  sin(«)  sin(^)  cos(]3  -  (p)  (17.25) 
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Substituting  the  expansion  of  the  Legendre  polynomial  given  by  Equation  17.1 1  into 
Equation  17.24,  one  arrives  at  the  following  expansion  of  a  plane  wave: 


expD'fc  •  r )  =  £(2/  + 1)/*  j, (kr)^T e(m)  m^P,m (cos (t?))P,m (cos(ar))cos(w(£ - (p)) 

1=0  m=0  (<  +  W)! 

Z  (2/  +  ty'Ji  (^)E(-l)mf  (cos(t?))/>m  (cos(a))cos(m(fi  -  <p)) 

1=0  m=0  ( l  +  m)\ 


(17.26) 


Expressing  this  equation  in  terms  of  the  regular  spherical  wave  functions  and  spherical 
harmonics,  one  arrives  at  the  following  expansion  of  the  plane  wave: 


exp[+/£  •  r  ]  = 


/=0  m=0(T=e 


/= 0  m~0cr=e 


(17.27) 


Using  the  orthogonality  condition  of  the  spherical  harmonics  and  the  expansion  of  a  plane 
wave  in  terms  of  spherical  wave  functions  given  by  Equation  17.27,  one  arrives  at  the  following 
plane  wave  representation  of  the  regular  spherical  wave  functions: 

|  2  K  +x 

Re  Vami  0)  =  t -J  \dP  jsin(a)da exp[+/£  •  ?]Yaml  (a,  p)  (1 7.28) 

0  -/r 

j  2  n  +* 

Re  ft*  00 = — r  Jsin(or) Ja  exp[+/&  •  r]!'*  (a,  /?)  (1 7.29) 

0 


The  irregular  (outgoing)  spherical  wave  functions  have  the  following  integral 
representation: 

\  lK 

Von, OO  =  T~T  \dP  Jsin(«)<7«exp[+/^  •  r]Yml (a, p)  (1 7.30) 

l7U  0  c+ 

|  In 

Von, (r)  =  — T  \dp  Jsin(or) daexp[+ik  •  rtf'* (a, /?)  ( 1 7.3 1 ) 

0  c± 

The  contour  C+  is  used  in  the  case  of  the  inequality  z  >  0,  and  the  contour  C_  is  used  in 
the  case  z  <  0 . 

The  proof  of  Equation  17.30  follows.  From  the  addition  theorem  of  the  spherical 
harmonics  and  the  definition  of  the  angle  y  in  Equation  17.25,  one  arrives  at  the  following 
transformation: 
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P,(cos(y))  = 


4£_yy 

)]  J_  1 

*«Tl  nj-D  rr=r 


(17.32) 


Equation  17.32  rotates  the  vector  {d,(p)  to  (0,0) .  Using  the  fact  that  the  spherical 
harmonics  form  a  linear  representation  of  the  rotation  group,  one  may  represent  the  effects  of  the 
rotation: 


^,(p)  =  Rz{Y)Ry{l5')Rz{a'){  0,0) 

on  the  spherical  harmonics  by  the  following  unitary  transformation: 

/  o  (0 

*/(£>*;  («')=  Z  HD  {CC\P\Y)Y.W, 

m'-Ocr’-e  a'm'.om 


(17.33) 


(17.34) 


£><'? , 
om.om 


The  rotation  coefficients  in  Equation  17.34  are 

(.or,?,?) = wifely uw 


given  by  the  following  expression: 


(17.35) 


From  Equations  17.32  through  17.35  one  arrives  at  the  following  result: 

?<*(<*,  p)  =  Y^i^P^co&iy))  (17.36) 

Substituting  Equation  17.36  into  Equation  17.30,  one  arrives  at  the  following  expression: 


1  in 

—~f  \dp  \sm(a)daexv[+ik»r]Yam,(a,p)  = 

0  c± 


|  2x 

— -Yml(0,<p)  \d/3  Jsin(a)t/aexp[+/^cos(y)]P/(cos(7)) 
2-m  o  c+ 


(17.37) 


Making  a  change  of  integration  form  the  variables  (cz,fi)  to  the  variables  (y,fi")  one 
obtains  the  result: 


1 


2m 


in 

7  \dp  Jsin(«)Jaexp[+/^  •  r]?^,  (a,  p)  = 


0  ct 


1  2;r 

—jY^i^cp)  jdP”  Jsin(7)^7exp[+z/rcos(7)]P/(cos(7)) 

a  r 


(17.38) 


o  c± 


The  above  integrand  is  independent  of  the  variable  P"  and  the  integration  over  this  variable  is 
trivial: 
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\  2" 

\dp  Jsin(ar)<foexp[+/£  •r]Yaml(a,/3)  = 

^7  Yami  (#>  <P)  Jsin(y)^7exp[+/^  cos(7)]P,  (cos(7)) 

Making  the  change  of  variables 
Z  =  cos(7) 

in  Equation  17.39,  one  arrives  at  the  following  integral: 

—j  J dp  Jsin(ff)<taexp[+/£  •  ?]Yml(a,p)  = 

0  c± 

^rYml^,V)ldZexV[+ikrZ]P,(,Z) 


(17.39) 


(17.40) 


(17.41) 


Substituting  the  Rodriguez  Formula  for  the  Legendre  polynomials  from  Reference  9  into 
Equation  17.41: 


P,(Z)  = 


_j_  d[_ 

2 '/!  dZ‘ 


(Z2-\)' 


(17.42) 


and  integrating  by  parts  one  arrives  at  the  following  expression  for  the  integral: 
1  2* 

—r  \dp  \sm{a)dazxv[+ik  •  r]Yoml(a,p)  = 

0  c± 

^•^-n./W^j^rapI+lfe'ZKl-Z’)' 


(1743) 


Using  the  following  integral  representation  of  the  spherical  Hankel  Function: 

4V)  =  ^  \dZ exp[rxZ](l  -Z2)1  (17.44) 

*  (OO 

one  arrives  at  the  result: 


\  l7t 

Vomir)  =  hf'OcryY^ip,?)  =  —  J dp  Jsin(a)^orexp[+/A  •  rjF^O*,/?)  (17.45) 

0  ct 

A  similar  derivation  follows  for  Equation  1 7.3 1 . 
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The  free  field  Greens’  Function  has  the  following  expansion  in  terms  of  the  spherical 
wave  functions: 

G(r, r':k)  =  exp (ik  |  r  -  F'|)  / 4?r  \  r  -  r' '|=  “ (k  \  r  - r ’|)  =  ik£ Wami  (r> ) Re Vomt  (r< )  ( 17-46) 

Z! 


This  expansion  follows  from  the  following  addition  theorem  for  Bessel  Functions  from 
Reference  9  evaluated  at  v  =  1  /  2 : 

H™(w)lwv  =2vr(v)f(v  +  k)H^k}u)  (17.47) 

£ 0  u  u 


w  =  ^/w2  +  w'2  -2uu'cos(a)  (17.48) 

the  following  identities  and  the  addition  theorem  for  spherical  harmonics  given  by  Equation 
17.32  are: 


(17.49) 

(17.50) 


C<I/2)(*)  =  i>(jc) 


(17.51) 
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18.  TRANSFORMATION  BETWEEN  SPHERICAL  AND  CYLINDRICAL  WAVE 

FUNCTIONS 


This  section  discusses  the  transformations  between  spherical  and  cylindrical  wave 
functions  of  the  free  field  Helmholtz  Equation. 

The  transformations  from  spherical  to  cylindrical  wave  functions  are  given  by  the 
following  expressions: 


Re  *®  {a)  =  X  C*  (<*)  Re  O8-1) 

l~m 

=  B<«(x-a)  =  (±i)'-J2/2+1  |'~”|y(cos(«))  (18.2) 

Equation  18.1  follows  from  the  expansion  of  the  free  field  Greens’  Function  in  terms  of  spherical 
and  cylindrical  wave  functions. 

The  transformations  from  cylindrical  to  spherical  wave  functions  are  given  by  the 
following  expressions.  The  upper  contour  is  used  in  Equations  18.4  and  18.5  if  the  z-coordinate 
is  positive,  and  the  lower  contour  is  used  if  it  is  negative: 


Re  Wami  =  Jsin  (a)daB™  (a)  Re*®  (a)  ( 1 8.3) 

0 

rm,=2|sin(a)rfas:-;,(a)Re^>(a)  (18.4) 

c± 

Womi  =  2  \sm{oc)daB^ (a) Re *£  (a)  (18.5) 

C- 

Worm  =  Jsin  (a)daB™  (a)*®  (a)  (18.6) 

c 

Equation  18.3  follows  from  Equations  18.1  and  18.2  and  the  orthogonality  of  the 
associated  Legendre  polynomials.  The  remaining  transformations  follow  from  the  expansion  of 
the  Greens’  Function.  Note:  Equation  18.5  follows  from  Equation  18.4  under  the 
transformation : 

a^n-a  (18.7) 
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in  the  integrand.  A  discussion  of  these  transformations  may  be  found  in  References  10  through 
12  by  the  author. 

The  proof  of  Equation  18.1  follows.  Equations  16.12  and  17.27  give  the  expansion  of  a 
plane  wave  in  terms  of  cylindrical  and  spherical  wave  functions.  Equating  these  two  expansions, 
one  arrives  at  the  following  relationship  between  cylindrical  and  spherical  wave  functions: 


+PQ  +/  O 


exp[+z*  •  r]  =  4;r£ /'  £ X  Re y/^,  (r)^ (a,  fi)  = 

/=0  m=0  o~e 

4^I(±0"Re^(r:a) 


m=0  a 


e(jn )  [cos(/»yg),<7  =  e] 
2 n  [sin(m^),cr  =  oJ 


(18.8) 


Using  the  orthogonality  of  the  trigonometric  basis  set, 

e{m)  cos{mp),o  =  e 
2k  sin  {mp),a  =  o 


(18.9) 


Equation  18.1  is  obtained  upon  multiplying  both  sides  of  Equation  1 8.8  by  the  basis  given 
in  Equation  18.9  and  integrating  with  respect  to  /?  from  0  to  2n . 


Equation  18.3  follows  from  Equation  18.1  and  the  orthonormality  condition. 


Jsin  (a)dccB™  (a)B™  (a)  =  Sf 
0 


Equations  1 8.4  and  1 8.5  follow  from  the  expansions  provided  by  Equations  1 6.56  and 
17.46  of  the  free  field  Greens’  Function  in  terms  of  cylindrical  and  spherical  wave  functions: 

G(r,  F:  k)  =  ik £  yfm,  (r> )  Re  y/ml  (r< ) 

am I 

=  likY,  jRe^’ (n. : a) Re^ (^  : a)sm(a)da  (1 8. 1 1 ) 

am  c± 

=  2ik'Z  :«)Re^,(r<  :a)sin(a)da 

am  r 

The  upper  contour  is  used  in  the  case  the  quantity  (z>  -z<)>  0  is  positive,  and  the  lower 
contour  is  used  if  the  quantity  is  negative.  Substituting  the  expansion  of  the  regular  cylindrical 
wave  function  in  terms  spherical  wave  functions  provided  by  Equation  18.1,  one  arrives  at  the 
following  expression: 
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G(?,P:t)  =  /*2;r„,(r>)  Re  !/„,(>•«) 

oml 

=  jRe  (r, :  («)  Re  if/^,  (rK )  sin(oO<ta  (18.12) 

om/  c± 

=  2z*£  jRe^ (r, : «’)^> (or) Re ^ (rjsin^or 

om/  c_ 

Equating  the  coefficients  of  the  regular  spherical  wave  function  in  Equation  18.12,  one 
arrives  at  the  following  expressions  for  the  representation  of  the  outgoing  spherical  wave 
function  in  terms  of  regular  cylindrical  wave  functions: 

Vonlir) 

=  2jReX^(r-.cx)B(-)(a)sm(a)da  (18.13) 

c± 

=  2  jRe^  (r :  a)B^  (a)sm(a)da 
thus  proving  Equations  18.4  and  18.5. 

Equation  18.6  follows  from  Equations  16.52  and  17.46,  which  provide  the  following 
expansion  of  the  free  field  Greens’  Function  in  terms  of  cylindrical  and  spherical  wave  functions: 

G(r, ?':k)  =  xff^ (r> ) Re y/^ (r< )  =  ik^  :a)Re^)(r<  :a)sin(a)da  (18.14) 

ami  am  c 


Substituting  the  expansion  provided  by  Equation  18.1  of  the  regular  cylindrical  wave 
functions  in  terms  of  spherical  wave  functions  one  arrives  at  the  following  expansion  of  the  free 
field  Greens’  Function: 

G(r,r':  k )  =  ^^(rJRe^rJ  =  jz£(r>  :  a)B$ (a) Re ^ (rK ) sin(a)da  (18.15) 

ami  am  c 


Equating  the  coefficients  of  the  regular  spherical  wave  functions,  one  arrives  at  the 
following  expression: 

Won,,  to  =  JjJT«  ir :  a)B™  (a)  sin  (a)da  (18.16) 

c 

thus  proving  Equation  18.6. 
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19.  EXPANSION  OF  NORMAL  MODE  IN  TERMS  OF  A  SPHERICAL  BASIS 


This  section  describes  the  expansion  of  the  normal  mode  representation  of  propagation  in 
a  waveguide  in  terms  of  a  spherical  basis  set  about  the  source  and  receiver  points.  Throughout 
this  section,  one  will  assume  the  density  and  sound  speed  are  piece-wise  constant  as  described  in 
Sections  8  and  9. 

From  Equations  9.18  and  9.19  the  normal  mode  representation  of  the  Greens’  Function  in 
a  homogeneously  layered  waveguide  with  a  homogeneous  half-space  as  a  basement  is  given  by 
the  following  normal  mode  and  branch  cut  contributions: 


G(r,zs,z ) 


i 

4 pM 


YiH"\q„r)F(z,-.q,)F(z-.q„)  +  Gc„(r,l,a) 


(19.1) 


1 

Gc^r’Zs,z)  =  —  \hNdhNH"\qr){Gz{z,zs  :+hN)-G:(z,zs : -hN )} 


(19.2) 


The  following  function  is  the  depth  function  of  the  normal  mode  in  the  n’th  layer: 


F(z, q )  =  A+F;  (z  :  q)  +  A~„F~  (z :  q), zn  <  z  <  z„+I 


(19.3) 


The  functions: 

F;  (z:q)  =  exp(±ihn  (z  -  zn ))  (19.4) 

are  the  basis  set  of  upward  and  downward-going  waves  in  the  n’th  layer,  and  the  coefficients 
{A*  :  n  =  0,1... TV}  are  the  coefficients  of  the  depth  functions  of  a  normal  mode  as  described  in 
Sections  8  and  9.  The  function  G;(z,zs )  is  the  depth  dependent  Greens’  Function  defined  in 
Section  8  and  9.  In  the  case  of  a  rigid  basement  the  cut  contribution  vanishes. 

Using  the  expansion: 

«“’(»■)=  f,HZ(u)JtW)cosm  (19.5) 

k=^<> 

w  =  ^u2  +«'2—  2mw'cos(^>)  (19.6) 
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from  Reference  9,  and  the  following  properties  for  integral  order  Bessel  Functions: 

^(x)  =  H)‘^°(x)  (19.7) 

J-k(x)  =  (-»kJk(x)  (19.8) 

one  obtains  the  following  expansion  of  the  normal  mode  contribution  to  the  Greens’  Function 
about  the  source  point: 


4  p{zs) 


H"\qa\r-*?s  | )F(z : qn )F(zs  +  Az s )  = 


i  00 

JZ7Z  )  (Cln^s  )F(Z  '  )F(Z,  +  Az,  )COS(tfjA^  ) 

** P\ZS  ) 


(19.9) 


Here,  (Ars,Azs,A(ps)  are  the  polar  coordinates  of  the  displacement  of  the  point  about  the  source 
point  (0,z5,0) ,  where  the  orientation  of  our  coordinate  system  is  such  that  the  source  and 
receiver  points  lie  in  the  (p  =  0  plane  of  the  coordinate  system.  Expanding  the  depth  functions  in 
terms  of  upward  and  downward-going  plane  waves,  one  may  express  Equation  19.9  in  terms  of 
the  product  cylindrical  wave  functions  about  the  source  and  receive  points: 


2m 


H"\qn\?-AFs  |  )F(z:qn)F(zs  +A zs)  = 


P(zs) 


SEE^rcv(^  ’  9n)Gsrc(Zs  •(ln)Xom,Nro,(.r$$  •  (A^  ,  Az^  ,  ’q„) 


(19.10) 


s'  om  .v=+j'=+ 


The  following  functions  are  the  coefficients  of  the  expansion  of  the  depth  functions  in 
terms  of  upward  and  downward-going  plane  waves  about  the  source  and  receiver  points,  where 
Nsrc  and  Nrcv  are  the  indices  of  the  layers  within  which  the  source  and  receiver  are  located. 

:  $„)  =  A±N„FL  ( Zs-<ln )  (19.1  1) 

GL  (z-qn)  =  an„.  (z:qn)  (19.12) 


Define  the  following  generalization  of  the  expansion  coefficients  of  the  regular 
cylindrical  wave  functions  in  terms  of  regular  spherical  wave  functions  in  the  N’th  layer: 


5m,jv(?)  =  (±  0 


l-m 


2/  +  l(/-m)1,  ,,  dm 

K)  *rPm 


x=hy(q)/ks 


(19.13) 


Equation  19.13  is  equivalent  to  Equation  18.2  evaluated  at  the  complex  angle  aN  defined  by 
Equation  19.14,  where  kN  is  the  wave  number  in  the  N’th  layer: 
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q  =  kNsm(aN) 

In  this  case  the  vertical  wave  number  is  given  by  the  following  expression: 

K  (q)  =  =  kN  c°s  (aN ) 


(19.14) 


(19.15) 


Substituting  the  following  expansion  of  the  regular  cylindrical  wave  functions  in  terms  of 
regular  spherical  wave  functions  into  Equation  19.12: 


R-e  X am, Nsrc  ~  ml,N„c  0? )  V ami,  N,r< 


(19.16) 


one  arrives  at  the  following  expansion  of  the  normal  mode  contribution  to  the  Greens’  Function 
in  terms  of  a  spherical  basis  set  around  the  source  point: 


.  ,  (g.lf-AF,  \)F(z:q„)F(zs+Azs)  = 
4  pizs) 

^aml,Nsrc  (r,  zs,z:qn)  Re  Wami,Nm  >  A<Ps ) 


(19.17) 


am I 


A^.NJr,z„z:q,)  =  ^f-Yt G : «■  >G« <■*. : ». W>.® : «. WU. («. )  < 1  »• 1 8) 

P\Zs)  S=+i'=+ 

Simplifying  Equation  19.18,  one  arrives  at  the  following  expression: 

Aaml,N(r,zs,z:qn)  = 


(19.19) 


The  odd  parity  terms  vanish  in  Equations  19.18  and  19.19  because  one  has  chosen  the  orientation 
of  our  coordinate  system  such  that  the  source  and  receiver  points  lie  in  the  (p-  0  plane. 

The  expression  for  expanding  the  normal  mode  contribution  about  the  receiver  point 
instead  of  the  source  point  is  made  by  exchanging  the  source  and  receiver  depths  in 
Equations  19.18  and  19.19,  except  in  the  density  term  in  these  equations. 

The  expression  for  expanding  the  normal  mode  contribution  about  both  the  source  and 
receiver  points,  where  (Ars,Azs,A<ps)  is  the  displacement  of  the  source  point  and  (Ar,Az,A<p)  is 

the  displacement  of  the  receiver  point  follows  from  expanding  Equation  19.19  about  the  receiver 
point. 
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(\r+tf\,zs,z  +  Az\qn)  = 


r  +  Ar\)FN^(z+Az:qn) 


iALFL (*s  ■  <in)K,,»„  (*. : qn)B-al„  } 


(19.20) 


Using  Equations  19.3,  19.4,  and  19.5  to  expand  the  Hankel  Function  and  depth  functions 
in  Equation  19.20  about  the  receiver  point,  one  arrives  at  the  following  expression: 


A aml.N „  (|r  +  Ar\,zt,z  +  Az:qn)  = 

K  ^nrVMn^)  COS (k(7C  ~  A<p)) 

iALFL  (z :  ^)exp(+/^Az)  +  ^F^(z :  qn)txp{-ihNm  Az)} 

lALFL  (z*  '■  +  a'n„fs„  (z,  ■  qn)Ki.Nm ) 


(19.21) 


The  factor  of  n  in  the  cosine  function  arises  from  the  choice  of  orienting  our  coordinate 
system  such  that  the  direction  of  the  line  from  the  source  point  to  the  receiver  point  is  in  the 
direction  <p  =  0 .  One  can  make  the  following  rearrangement  of  the  summation  in  Equation  19.21 : 

+oo 

Z  Hmlk  fa*V*(?Ar)  cos(Ar(;r  -  A  <p))  = 

*=-« o 

mL (W) + (- l)m  H (qr ))J„. (qAr) cos(/n'  (;r  -  A<p))  =  (19.22) 

#«‘= 0  * 

Z  (~  1) m  (#£+«•  (qr) + (-  l)m  H^n.  (qr  ))Jm.  (qAr )  cos(m'  A  <p) 

m'=0  z 


Upon  substitution  of  Equation  19.22  into  Equation  19.21,  one  obtains  the  following 
expression  for  the  expansion  about  the  receiver  point: 


(\r+Ar\,zs,z  +  Az:qn)  = 

(ALFL (z : q»)exp(+i\mAz)  +  A-^ (z :  q„)txp(-ihN 
{An„Fn„  (zj  •  qn)Bml^  +  An^Fn^  (zs  :  q„)BmlN  } 


(qnr))Jm.  (q„Ar)  cos {m\A(p)) 


(19.23) 


Making  the  following  substitution  into  Equation  19.23: 


Re^.m.W(CT(Ar,Az,A^) 


J  m-  (q^f )  exp(±//?Az)  { 


cos(m'A(p),o'=e 
sin(w'A^),  cr'=o 


(19.24) 
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one  arrives  at  an  expansion  in  terms  of  the  regular  cylindrical  wave  functions  about  the  receiver 
point: 


WOT(IF+A?l’z^z+Az:^)= 

-  O  -H»  f 

it 


s'=+o'=em'= 0 


2m  \e(m) 


Hr  4*. 


p(z,)Af  8/r  '  '  V  8?r 
M*  /?(z,  +A-sirF~JZs  :qn)B-ml'„j4mFZjz:qn) Re^JAr.Az.A*) 


(19.25) 


Substituting  the  expansion  of  the  regular  cylindrical  wave  function  in  terms  of  the  regular 
spherical  wave  functions  about  the  receiver  point  one  obtains  the  following  expansion: 


Aoml,NJ\r  +  Ar\>zs,z  +  Az:qn)  = 


aml,Nsrc  ,a'm'l\Nrc 


( r ,  z, ,  z)  Re  WmTiA,  (Ar,  Az,  Ap)  (19.26) 


G'm'V 


•  oml  ,Nsrc<a'm'l\Nrc 


P(Zs) 


(z, : <?„)i^  +  (*.  :  } 

(z : 9.  (z  •  9.  > 


(19.27) 


The  overall  factor  of  (-l)m*  in  Equation  19.27  is  due  to  the  fact  that  the  orientation  of  the 
coordinate  systems  at  the  source  and  receiver  points  are  oriented  such  that  the  line  from  the 
source  to  the  receiver  point  is  in  the  direction  (p  =  0 ,  that  is  the  orientation  of  the  local  coordinate 

systems  at  the  source  and  receiver  points  are  parallel. 

Substituting  Equation  19.26  into  19.17,  one  arrives  at  the  following  expression  for  the 
expansion  of  the  normal  mode  component  about  both  the  source  and  receiver  points: 

■  7  ~ H ^ (qn\F  +  Ar-Ars  \)F(z  +  A z:qn )F(zs  +  A z, )  = 

*P(zs)  (19.28) 

(rtz„z : q„) Re^  (Ars , Az^A^Re^,^  (Ar,Az,A<p) 

ami  a'mT 

Here,  (Ar^Az^A^)  is  the  displacement  of  the  source  point,  and  (Ar,Az,A<p)  is  the 
displacement  of  the  receiver  point. 

A  similar  expansion  may  be  made  for  the  integrand  of  the  branch  cut  contribution,  where 
the  coefficients  {A^}  of  the  depth  functions  are  replaced  by  the  expansion  coefficients  {4,4} 
of  the  depth  dependent  Green’s  Function  in  the  integrand. 
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The  generalization  of  Equations  19.19  and  19.27  to  local  coordinate  systems  whose  z-axis 
are  vertical  but  whose  orientation  of  the  x-axis  is  not  parallel  to  the  line  from  the  source  to  the 
receiver  are  given  by  the  following  expressions: 


+a^fnJ2s } 


(19.29) 


m  sin  (m<ps),(X  =  o 

(*  1.  )«;,»„ + 4i.fi.  <* :  9.  w*. }  e } 

sm(wV),<r'=o 


(19.30) 


Here  p,  is  the  angle  between  the  x-axis  of  the  local  coordinate  system  at  the  source  and  the 
location  of  the  receiver  as  seen  from  the  source,  and  q>  is  the  angle  between  the  x-axis  of  the 
local  coordinate  system  at  the  receiver  and  the  location  of  the  source  as  seen  from  the  receiver. 
Note,  in  Equations  19.19  and  19.27  these  two  angles  have  the  following  values:  (ps  =  0 ,  and 

(p  =  Jt . 
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20.  METHOD  FOR  INCLUDING  DIRECTIVITY  OF  THE  SOURCE  AND  THE 
RECEIVER  IN  THE  NORMAL  MODE  TERM 


This  section  uses  the  results  of  Section  19  to  include  the  effects  of  the  directivity  of  the 
source  and  the  receiver  in  the  normal  mode  description  of  propagation. 

Let  DSrc(&,<p)  denote  the  beam  pattern  for  the  source,  and  DRcv(&,<p )  denote  the  beam 

pattern  of  the  receiver.  Using  the  orthogonality  of  the  spherical  harmonics,  one  can  project  the 
beam  patterns  for  the  source  and  receiver  onto  the  spherical  harmonics  by  performing  the 
following  surface  integrals: 

ds^i  =  \d(p\sm(p)d  (d,  (p)DSrc  (d,  (p)  (20.1) 

0  0 

(20.2) 

0  0 

These  integrals  can  easily  be  performed  for  arbitrary  beam  patterns  by  using  Gauss 
Legendre  integration  of  sufficiently  high  order.  Using  the  expansion  of  the  normal  mode 
contribution  about  the  source  and  receiver  points  given  by  Equation  19.28: 


— — r # o°  (?  J  f + AF  -  A rs  \)F(z  +  A z:qn  )F(zs  +  A zs)  = 

4p(zs)  (20.3) 

(r,zs,z,(ps  ,(p :  (ZjRe^^JAr,,  ,  A^Re^^JA^Az.Ap) 

ami  a'm'l' 

one  may  project  the  above  normal  mode  contribution  onto  the  directivity  pattern  of  the  source 
and  receiver  by  performing  the  following  summations: 


X  X^om/.v,,,.  (r,zs,z,(ps,(p:q)dSrccmld 

Rcv.a'm'l' 

ami  a'm'F 


(20.4) 


The  coefficients  {T^itNm  ^mr>Nrcv(r,zs,z,q)s,<p:q)}  are  defined  in  Equation  19.30. 

In  the  case  the  receiver  is  omni  directional,  the  projection  of  the  normal  mode  onto  the 
directivity  of  the  source  is  given  by  the  following  summation: 

X zs’z><Ps :  ^dSrcaml  (20.5) 

ami 
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The  coefficients  {Aml  ( r,zs,z,(ps :  q)}  are  the  expansion  coefficients  of  the  normal  mode 

contribution  about  the  source  point  defined  in  Equation  19.29. 

Equations  20.4  and  20.5  assume  the  coordinate  system  used  in  computing  the  projection 
of  the  source  and  receiver  beam  pattern  onto  a  spherical  basis  are  parallel  to  the  coordinate 
system  used  in  computing  the  coefficients  in  Equations  19.29,  and  19.30.  In  the  event  these 
coordinate  systems  are  not  parallel,  the  computation  of  projection  of  the  source  and  receiver 
beam  patterns  onto  a  spherical  basis  set  will  require  rotating  the  origin  of  the  argument  {d,(p)  of 
the  directivity  functions. 
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21.  ROTATION  OF  COORDINATE  SYSTEMS 


This  section  describes  the  procedure  for  rotating  the  local  coordinate  systems  at  the 
source  and  the  receiver  for  an  arbitrary  orientation.  A  practical  introduction  to  the  theory  of 
representations  of  the  rotation  group  “SO(3)”  is  provided  by  Tinkham13.  More  advanced 
treatises  on  representation  theory  are  provided  in  texts  by  Helgason14  and  Weyl15. 

The  conventions  of  Rose  are  used  to  describe  an  arbitrary  rotation  as  a  series  of  three 
rotations  about  the  fixed  axes  (x,y,z) : 

R(a,  f5,Y)  =  Rz(a)oRy{p)oRz{y)  (21.1) 

The  equivalent  rotation  about  relative  axes  is  given  by  the  following  set  of  three  rotations 
about  the  rotated  axes: 

R(a,p,Y)  =  R2{a)°Ry{p)oRz{Y)  =  R:„{Y>Rym°RM  (21-2) 

Here,  (x\y',z')  are  the  axes  obtained  by  rotating  the  axes  (x,y,z)  by  the  rotation  R,  (y)  about 
the  z-axis,  (x",y",z")  are  the  axes  obtained  by  further  rotating  the  axes  by  the  rotation  Ry(fi) 
about  the  y'  -axis,  and  (x"',/",z"')  are  the  axes  obtained  by  further  rotating  these  axes  by  the 
rotation  R,» (cc)  about  the  z”-axis. 

The  complex  spherical  harmonics  are  defined  as  follows: 

Yml(&,<p)  =  J2/a+1  77  Qxp(+im<p)Plm (cos(ff)), m  =  -/,-(/ - 1)... + /  (21.3) 

\  An  ( l  +  my. 

Assume  the  Condon-Shortley  phase  convention  for  negative  indices  given  below: 
p-m  (cos(tf))  =  (-l)m  P,m  (cos(r?))  (2 1 .4) 

Y.ml(#,(p)  =  (-\)mY'mlW,(p)  (21.5) 

The  spherical  harmonics  (Fm/(z?,^>) :m  =  -/,-(/ -1)...+/}  forma  (2/ + 1) -dimensional 
irreducible  representation  of  the  rotation  group  “SO(3)”  on  the  homogeneous  manifold 
SO(3)/SO(2)~S2 .  The  spherical  harmonics  obey  the  following  transformation  rules  under  an 
arbitrary  rotation: 


j?(«,  >5,  r)  *  )(^, ^  r) 


m  — -/ 


(21.6) 
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The  Wigner  rotation  coefficients  D(J?m,(a,0,y)  are  given  by  the  following  relationships 
from  Reference  13: 

Dm 2»  («.  P>  r)  =  exp(-im'a)  exp  (-im  y)d(J?m  (0)  (2 1 .7) 


d%m{P)  = 

y  (-i)* 


(cos{0/2))2J-n-m'+m 


(-sm(/3/2))2k+m'-m 


(21.8) 


The  summation  in  Equation  21 .8  extends  overall  possible  values  of  the  index  k  for  which 
the  denominator  is  finite,  that  is,  k  satisfies  the  following  inequalities: 


k>  0 

j+m-k> 0 
j  -  m'-k  >  0 
k  +  m'-m  >  0 


(21.9a) 

(21.9b) 

(21.9c) 

(21.9d) 


Equation  21 .9  results  in  the  following  limits  for  the  index  k: 
min(y  +m,j-m')>k  > max(m - m' ,0)  (21.10) 

The  Wigner  coefficients  form  an  irreducible,  unitary  representation  of  the  group 
Spin(3)~SU(2)  corresponding  to  the  covering  group  of  the  rotation  group  “SO(3)’\  “SU(2)”  is 
the  Z2 -valued  covering  group  of  the  rotation  group  “SO(3)”  associated  with  the  spinor 
representation  of  the  rotation  group.  One  of  the  properties  of  the  group  “SU(2)”  is  that  the 
rotation  of  a  spin  1/2  representation  by  2k  results  in  a  multiplication  of  the  spinor  by  —  1 . 

The  relationship  between  group  elements  in  “SU(2)”  and  “SO(3)”  is  given  by  the 
following  correspondence: 

u(a,  /J,y)eSU (2)  u(a,  p,  y)atu  (a,  0,  y)  =  Ru  (a, 0,  y)aj  (21.11) 

Here,  u(oc,0,y)G  SU( 2)  is  the  following  group  element: 

u(a,  0,  y)  =  f +  exp(+l(or  +  y)^)  cos(0  /  2)  +exp(+i(y-a)/2)sm(0/2)\ 

V  ~  exp(-/(y  -a)  12)  sm(0  /  2)  +exp(-/(or+y)/2)cos(/?/2)J  (  '  ’ 
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The  Hermitian  matrices  {<J, :  /  =  1,2,3}  are  the  following  Pauli  spin  matrices: 


r0 

(21.13a) 

oj 

'  0 

(21.13b) 

<+i 

0. 

r+l 

0) 

(21.13c) 

,0 

-ij 

The  matrix  RiJ(a,0,y)  is  the  three-dimensional  rotation  matrix  with  the  following 


elements: 

Rii(a,J3,y)  =  +cos(a)cos(/?)cos(y)-sin(ar)sin(}')  (21.14a) 

Rl2(a,/3,y)  =  + sin(a)  cos(/?)  cos(  y)  +  cos(or)  sin(  y)  (21.14b) 

Rl2(a,fi,y)  =  -sin(/?)cos(y)  (21.14c) 

R2](a,fi,y)  =  -cos(a)cos(/?)sin(y)-sin(a)cos(y)  (21.14d) 

R22(a,f5,y)  =  -sin(«)cos(/?)sin(y)  +  cos(a)cos(y)  (21.14e) 

R2  3  (a,  ft,  y)  =  +sin(fi)  sin  (y)  (2 1 . 1 4f) 

R3 ,  ( a,p,y)  =  +  cos  (a)  sin(J3)  (21.1 4g) 

R22(a,p,y)  =  +sin(«)sin(/?)  (21.14h) 

RXi(a,fi,y)  =  +cos(fi)  (21.1 4i) 


This  matrix  is  given  by  the  following  rotations  about  the  z-axis  and  y-axis: 


Rziy)Ry{p)Rz{a)  = 

cos(y)  +  sin(y)  0Y+cos(/?)  0  -sin(/?)Y  +cos(a) 
-sin(y)  +  cos(y)  0  0  10  -  sin(or) 

0  0  ljl^+sin(/?)  0  +cos(/?)J[  0 


+  sin(or) 
+  cos(or) 
0 


0^ 

0 

1 


(21.15) 
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The  correspondence  between  the  Wigner  coefficients  and  the  matrix  representation  of  the 
corresponding  group  element  is  given  by  the  following  relationship: 

Dl:l(«:0,y)=rJLUK^r)Kr(-P)KA-a))=Ul%<(R/a)RJP)RJy)y')  (21.16) 

Here,  ((/".( g)  :  m,  m'=  -/,...+ /}  are  the  matrix  elements  of  a  (21  + 1)  -dimensional  unitary 
representation  of  the  group  element  g  e  50(3)  c  SU (2) .  Strictly  speaking,  this  matrix  is  a 
representation  of  the  covering  group  “SU(2)’\  The  group  “SU(2)”  is  a  non-trivial  Z2  -valued 
principal  fiber  bundle  over  the  group  “S0(3)’\  However,  in  the  case  of  integral  angular 
momentum,  the  above  matrix  elements  are  constant  along  the  fiber  of  this  principal  bundle,  and 
the  matrix  elements  may  be  considered  to  be  functions  of  the  base  space  “S0(3)’\  The  constant 
nature  of  this  matrix  along  the  fiber  follows  from  the  fact  that  for  integral  angular  momentum  a 
rotation  by  multiples  of  2 n  leaves  the  representation  unchanged. 

The  action  of  the  group  element  ge  50(3)  acting  on  the  spherical  harmonics  can  be 
expressed  in  the  following  form  in  terms  of  this  unitary  representation: 

V0-)=I 

m' 

The  function: 

Lg  :  50(3)/ 50(2)  50(3)/50(2) 

Lg\.S'] H [5 •  5']  =  {g* g'*h :  he  50(2)} 
denotes  the  natural  left  action 
50(3)x50(3)  50(3)  ->  50(3)/ 50(2)  -  52 

(g,g')  g  •[£’]  =  [g  •  g'] 

of  the  group  element  g  e  50(3)  acting  on  the  homogeneous  space 
it :  50(3)  50(3)/ 50(2)  *  52 

g  ~>n(g)  =  [g]  ={5  •h-.he  50(2)} 
and  Lg  *  denotes  the  pullback  of  this  function. 

One  may  introduce  the  following  transformation  from  the  complex  spherical  harmonics: 

=  -/,...+/}  to  the  real  spherical  harmonics  :  o-e,o\m  =  0,..i} : 


(21.17) 

(21.18) 

(21.19) 

(21.20) 
(21.21) 

(21.22) 

(21.23) 
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i '_(*«>) = <*.«>) +(-d*  ««>)) 


(21.24) 


n-  o>,« = (X,  «>,«>)-  (-lr  z_  o>,  p» 


(21.25) 


The  effect  of  the  rotation  R(a,(3,y)  on  the  real  spherical  harmonics  is  of  the  form: 


R(aj,y)‘(.r„,x0,<p) = £ 


(21.26) 


The  matrix  elements  of  the  real  Wigner  rotation  coefficients  are  given  by  the  following 
relationships: 


(/)  _  1  je(m)  ,  a)  ,i\m'nu)  .+.  r-um+m' nU) 


2  X  £(/»•) 


{DZ  +  (-l r  DZ,m  +(-l)mX>:U  +(-l)ffl+m 


(21.27) 


(/)  +i_  £(m)_  (/)  _(-_1n«'£)(/)  } 


(21.28) 


D(/)  =—i  mw  +(_nm'£)(/>,  -(-lY”//'?  £)<'),  } 

LJem',om  -i  ,  ...  V  V  ^ m',-m  V 

2  \  £(m  ) 


(21.29) 


(/>  _  1  (e(m) 


(/)  ±  |D(/)  _<_iyDv)  } 

^om'.om  2  \j  f  (w')  '~LJ’n  'm  '  V  V  Lym\-m  V  V 


(21.30) 


From  Equations  21 .7  and  21 .8  one  may  derive  the  following  relationships  by  rearranging 
the  terms  in  the  summation  in  Equation  21.8: 


D%(a,p,y)  =  exp(-z'm'  a )  exp(-z'ro  y)d(J?m  (ft) 


(21.31) 


dZ,-JP)  =  (-W+n'dZ,+m(fi) 


(21.32) 


(21.33) 


Substituting  Equations  21.31  through  21.33  into  Equations  21.27  through  21.30,  one 
arrives  at  the  following  representation  of  the  real  Wigner  rotation  coefficients: 


D(Zem  =  ^^{d%(P)cos(m'a  +  my)  +  (-\)mPZm{p)cos(m'a-rny)} 

V  £(m ) 


(21.34) 
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Dom',em  =  J~^  Wi>  (/?)sin(m’tf  +  m  y)  +  (-  l)m'  m  (y?)  sinO’a  -my)} 


(21.35) 


WK.  (£)  sin(m'a + my)  -  (-l)m'^.m  (/?)  sin(m'a-  my)} 


(21.36) 


Dl<on,  =  (/?)cos(m’a  +  my)  -  (-l)m  m  (p)cos(m'a  -  my)} 


(21.37) 


Consider  the  case  j  -  1.  The  spherical  harmonics  are  given  by  the  following 
relationships: 


"  n*S  (*)e  ~  V&r  r 


(21.38) 


n,=+jAcos(,)=+JXi 


(21.39) 


r-=+fe  si„We-=J±£z£ 


(21.40) 


In  the  case  of  the  rotation  /?(<*, 0,0)  =  /?2(a)  about  the  z-axis,  the  Wigner  rotation 
coefficients  in  the  case  j=l  are  given  by  the  following  relationships: 

rexp(-/a)  0  0  ^ 

Z>i'?m(a,0,0)  =  0  +1  0  (21.41) 

,  0  0  exp(+za) 


In  this  case,  the  rotated  spherical  harmonics  are  given  by  the  following  linear  combination 
of  the  spherical  harmonics: 


YuoR:(a)  =  e-iaYu 
Y0loR;(a)  =  Y0] 


(21.42a) 

(21.42b) 


Y_u°R:(a)  =  e+iaY_u 


(21.42c) 


In  the  case  of  the  rotation  /?( 0,/?,0)  =  Ry(fi)  the  Wigner  rotation  coefficients  in  the  case 
j=l  are  given  by  the  following  relationships: 
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DZ(o,m= 

f  cos(/?/2)2  -  V2  cos(y(?  /  2)  sin(yS  /  2) 

+V2cos(/?/2)sin(/?/2)  cos(/?/2)2  -sin(/?/2)2 

^  sin(y5  /  2) 2  +  V2  cos(/?/2)sin(/?/2) 

r(l  +  cos(/?))/2  -V2sin(yS)/2  (l-cos(y5))/2N' 
+V2sin(/?)/2  cos(yff)  -V2sin(/?)/2 
(l-cos(/?))/2  +  V2sin(/?)/2  (l  +  cos(/?))/2 


sin(/?  /  2)2 

-  V2  cos(/?  /  2)  sin(/?  /  2) 
cos(/?/2  )2 


(21.43) 


The  rotated  spherical  harmonics  are  given  by  the  following  linear  combination  of  the 
spherical  harmonics: 


(xcvM-zitoU) )+»)  = 

+^*c«r.  +  , 


(21.44a) 


y  (xsin(/?)  +  zcos(/?))_ 

-^sin(«y„  +cos(«r„  +^ySin(/))Y_ll 


(21.44b) 


Y_uoRy^)  =  +y 


3  (*cos(/?)-zsin(/?)-*>) 


(21.44c) 
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22.  HELMHOLTZ  INTEGRAL  EQUATION 


This  section  reviews  the  Helmholtz  Integral  Equation  in  anticipation  of  the  next  section 
on  the  T-matrix  description  of  scattering  from  a  target. 


The  free-field  Greens’  Function  for  the  Helmholtz  differential  equation  is  given  by  the 
following  expression: 


G(r,r':k) 


exp(+/fc  |  r  —  r'|) 
Ait  \r-r'\ 


(22.1) 


This  function  satisfies  the  following  inhomogeneous  partial  differential  equation: 


(V  •  V  +  k2  )G(r,F:  k )  =  - S 3  (r-F) 


(22.2) 


Let  S„  u  S0  denote  the  bounding  surface  of  the  volume  V ,  where  the  surface  S„  is  the 
surface  at  spatial  infinity,  and  S0  is  the  union  of  the  surfaces  of  all  scatterers  in  the  volume.  Let 
the  function 


¥  =  V Incident  +  V Scattered 

denote  the  solution  of  the  Helmholtz  Equation  in  the  volume  V  resulting  from  the  scattering  of 
the  incident  field  ¥  incident  from  the  scatterers  interior  to  this  volume.  The  incident  field  is  given 
by  the  following  surface  integral  over  the  sphere  at  spatial  infinity: 

¥ incident (^" )  ~  •  {G(r,F:  *)V y(F) - G(r ,F: £)}  (22.4) 

Using  Greens’  Theorem,  one  obtains  the  following  integral  equation  for  the  total  field: 
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J< dA  •  {G(r,r':  k)V  \j/{r')  -VG(F,f:  k)y/(r')} 

-  jdA  •  {G(r, k)V  y/(r')  - VG{r,r':  k)\ff{r')} 

So 

=  ^dV{G{r,r':k)V •  V^(r')-V  •VG(r,r':  k)y/(r')}  (22.5) 

V 

=  ~^dV{W  •VG(r,F:  k)  +  k2G(r,r':  k)}y/(F) 
v 

y/(r),re  V 
^(r),redV 
0,riVudV 

Substituting  Equation  22.4  for  the  surface  integral  over  the  sphere  at  spatial  infinity,  one 
obtains  the  following  integral  equation  for  the  field  in  terms  of  the  incident  field  and  the  surface 
integral  over  the  surface  of  the  scatterers: 

VincidenXr)-  fdA •  {G(r  ,r': k)V y/(r') -VG(r,r':  k)y/(F)}  =  <  ^y/(r),redV  (22.6) 

{0,?eVvdV 

Equation  22.6  is  the  Helmholtz  Integral  Equation  for  the  acoustic  field. 
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23.  SCATTERING  FROM  A  RIGID  TARGET 


This  section  gives  a  basic  description  of  the  spherical  T-matrix  approach  to  scattering 
from  a  rigid  target.  One  of  the  first  articles  to  appear  on  the  use  of  the  spherical  T-matrix 
description  of  acoustic  scattering  is  the  article  by  Waterman16.  Hackman17  provides  a  review  of 
the  spherical  and  spheroidal  T-matrix  methods. 

In  the  following  derivation,  one  may  make  use  of  the  following  bilinear  expansion  of  the 
free-field  Greens’  Function  in  terms  of  spherical  wave  functions: 

G(r,  F:  k)  =  ik^Womi  (^>)Re  ¥<mi  </<)  (23- 1 ) 

amt 

Let  the  function 


¥  ¥  Incident  ¥ , Scattered 


(23.2) 


denote  a  solution  of  the  Helmholtz  Equation  representing  the  acoustic  field  in  the  presence  of  a 
rigid  target.  The  acoustic  field  is  required  to  satisfy  the  following  boundary  condition  on  the 
surface  of  the  rigid  target: 

n*V\j/  |s=0  (23.3) 


The  incident  and  scattered  fields  satisfy  the  following  Helmholtz  Integral  Equations: 

W incident^)  =  •  {G(r  ,F:  k)V  y/{F)  -  yr(F)V  G(r  ,F:  k)}  (23.4) 


¥ Incident (?) "  fa •  {G(r , F:  k)V yr(F) -  V G(r ,F:  k)y/(F)}  = 

So 


y/(r),reV 

^y/(r),redV 

0,r£VvdV 


(23.5) 


Substituting  Equation  23.3  for  the  normal  gradient  of  the  acoustic  field  on  the  surface  of 
the  target  into  Equation  23.5,  one  arrives  at  the  following  integral  equation: 


¥  Incident^)  +  fa 
So 


•  VG(F,F':)fc)^(P)H 


y/(r),reV 

^(F),Fe9F 

0,FeFu3F 


(23.6) 
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The  incident  field  and  scattered  field  may  be  expanded  in  terms  of  the  regular  and 


outgoing  spherical  wave  functions: 

^ Incident  ^  \  ®oml  R®  W and 

(23.7) 

am! 

W Scattered  “  ^  j  ^aml  W am! 

(23.8) 

am! 


The  scattered  field  on  the  surface  of  the  target  may  be  expanded  in  terms  of  the  regular 
spherical  wave  functions: 

^ Scattered  =  Z  PonlVontl  (23.9) 

am! 


Substituting  Equations  23.7  and  23.9  into  Equation  23.2,  one  arrives  at  the  following 
expansion  of  the  acoustic  field  on  the  surface  of  the  target: 

^ls=Z<a™/+A™/)Re^A  (23.10) 

am! 


Substituting  Equations  23.1  and  23.10  into  Equation  23.6  evaluated  at  a  point  interior  to 
the  target,  one  arrives  at  the  following  linear  equation: 


Reyu +XRe^/  +&•„•/•) = o 

am!  am!  o'm't' 


(23.11) 


The  matrix  Q  is  the  infinite  dimensional  matrix  defined  by  the  following  surface  integral  over  the 
surface  of  the  target: 

Qomi.a-m'r  =  ~  \\dA  •  V  y/ml  Re^ff.mV.  (23. 1 2) 

5 


Substituting  Equations  23.1  and  23.10  into  Equation  23.6  evaluated  at  a  point  exterior  to 
the  target,  one  arrives  at  the  following  linear  equation: 

"  2>-  X  Re  Qaml.a’mT  (aa'mT  +  PoWr)  =  0  (23. 1  3) 

am!  am!  a'mT 

The  matrix  ReQ  is  the  infinite  dimensional  matrix  defined  by  the  following  surface  integral  over 
the  surface  of  the  target: 

Re  Qomiw  =  “  \\dA  •  V  Re  y/ml  Re^w/,  (23.14) 

s 

The  solution  of  this  pair  of  infinite  dimensional  linear  equations  is  given  by  the  following 
expression: 
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^aml  '^'oml.o'm'vQo'm'V 


(23.15) 


a'm'V 


Here,  the  matrix  T  is  the  T-matrix  of  the  scattering  from  the  target  defined  by  the  following 
matrix  product 


^aml.a'm'V  ^ReQoml,o'm’rQ  lo"mT,a'mT 


(23.16) 


Gnm"r 


In  the  case  of  a  rigid  sphere,  the  T-matrix  is  diagonal  and  has  the  following  form  in  terms 
of  the  derivative  of  the  regular  spherical  Bessel  Function  and  Hankel  Function: 


rp  _ _ ocr'mT  j I  (M 

^aml.a'm'V  ami  ^(1),  (Jca) 


(23.17) 
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24.  SCATTERING  FROM  A  SOLID  ELASTIC  TARGET 


This  section  describes  the  scattering  from  a  solid,  elastic  target  using  the  spherical 
T-matrix  methods  described  by  Hackman17. 

In  order  to  describe  the  scattering  from  an  elastic  target  it  is  necessary  to  introduce  the 
vector  Helmholtz  Equation  for  an  elastic  solid.  Let  us  begin  by  describing  the  Helmholtz 
Equation  for  a  fluid  in  terms  of  Newton’s’  Second  Law  of  Motion  and  the  stress  tensor  for  a 
homogeneous,  isotropic  fluid: 


d2“  =v«f 

(24.1) 

dt 2 

=  -SyP 

(24.2) 

Equation  24.1  is  Newton’s’  Second  Law  of  Motion  for  a  fluid  in  terms  of  the  displacement 
(m.)  of  the  fluid  and  the  stress  tensor  of  the  fluid  (Ttj ).  Equation  24.2  is  the  stress  tensor  of  a 

homogeneous  isotropic  fluid,  where  P  is  the  pressure  of  the  fluid,  and  gy  is  the  metric  tensor  of 

the  coordinate  system.  Substituting  Equation  24.2  into  Equation  24.1  and  assuming  a  harmonic 
dependence  on  time,  one  arrives  at  the  following  expression  for  the  displacement  of  the  fluid  in 
terms  of  the  gradient  of  the  pressure: 

J  =  -L-V/>  (24.3) 

poo 

Analogous  equations  of  motion  for  a  homogeneous,  isotropic,  elastic  solid  are  given 

below: 


P 


=  V«f 


(24.4) 


Ty  =  gyXV  •  u  +  JU(V.Uj  +  VjU, )  (24.5) 

Here,  the  parameters  X  and  p  are  the  compression  and  shear  components  of  the  Lame’  constants 
of  the  material.  The  longitudinal  and  shear  velocities  are  given  by  the  following  expressions  in 
terms  of  the  Lame’  constants  and  the  density  of  the  material: 


A  +  2// 
V  P 


(24.6) 
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(24.7) 


Substituting  Equation  24.5  into  Equation  24.4,  one  arrives  at  the  following  equation  of 
motion  for  the  elastic  solid: 

=  (A + 2p)V (V  •  u) -//(V  x  (Vx  w)),. 

In  the  case  of  an  elastic  solid,  it  is  often  convenient  to  express  the  displacement  vector  in 
terms  of  the  Debye  potentials  O  and  A : 

u  =  VO+Vx.4  (24.9) 

Substituting  Equation  24.9  into  Equation  24.8,  one  arrives  at  the  following  equation  of 
motion  for  the  Debye  potentials: 

P^y(VO  +  VxA)  =  a  +  2  p)V(V  •  VO)  -  p(V  x  (V  x  (V  x  ^)))  (24. 1 0) 


This  equation  of  motion  is  equivalent  to  the  following  pair  of  differential  equations,  since 
the  longitudinal  and  transverse  components  of  the  vector  in  Equation  24.10  are  uncoupled: 


p-^VO)  =  (A + 2p)V(V  •  VO) 
at 

(24.11) 

j2  __ 

p— y(Vx^)  =  -p(Vx(Vx(Vx^))) 
at 

(24.12) 

Extracting  an  overall  gradient  from  Equation  24.1 1,  and  curl  from  Equation  24.12,  one 
may  express  this  pair  of  equations  in  the  following  form: 


j2 

p— yO  =  (A+2p)V»VO 
at 


(24.13) 


p^y^  =  -//(Vx(Vx^)) 


(24.14) 


Let  u  and  vv  be  a  pair  of  solutions  of  the  equations  of  motions  for  an  elastic  solid 
occupying  the  volume  V .  The  following  surface  integral  over  the  boundary  of  this  volume  is 
equal  to  zero: 


24-2 


CSS/TR-02/10 


WldViufljr  (w)  -  w,VjTiJ  («)}  =  \\\dWJ  {uiTij  (w)  -  w.p  (w)}  = 

V  V 

fjdArij  {u,T9  (w)  -  w,Tv  (u)}  =  \\dA{u  •t(w)-w*l{u)}  =  § 


(24.15) 


Here,  the  vector 
t,(u)  =  TfWnJ 


(24.16) 


is  the  surface  traction  at  the  boundary  of  the  volume  V,  where  the  vector  n  is  the  outgoing 
normal  to  the  surface  dV .  Equation  24.15  is  the  third  Betti  Identity  for  the  equations  of  motion 
for  the  elastic  solid.  This  identity  is  used  extensively  in  the  construction  of  the  T-matrix 
description  of  an  elastic  target. 

The  boundary  conditions  between  a  pair  of  elastic  solids  are  that  the  displacements  and 
surface  tractions  be  continuous  at  the  boundary. 


u+  =  u_ 


t(u+)  =  t(u_) 


(24.17) 

(24.18) 


Here,  u+  and  u_  are  the  displacement  vectors  on  the  two  sides  of  the  boundary  evaluated  at  the 
boundary.  In  the  limit  the  shear  modulus  of  both  solids  approach  zero,  the  above  boundary 
conditions  become  the  following  boundary  conditions  between  two  fluids,  respectively. 


— n»VP+  =  —fi»VP_ 
P+  P- 


P+=P_ 


(24.19) 

(24.20) 


The  boundary  conditions  between  a  fluid  and  a  solid  are  that  the  normal  components  of 
the  displacement  and  surface  traction  are  continuous,  and  the  tangential  surface  traction  of  the 
solid  vanishes  at  the  boundary. 


m+  •  n  =  u_  •  n 

(24.21) 

t(u+)»n  =  ?(«_)•« 

(24.22) 

?(w+)xn  —  f(«_)xn  —  0 

(24.23) 

The  tangential  surface  traction  of  an  inviscid  fluid  is  zero.  Equations  24.3  and  24.2  give 
the  displacement  vector  and  stress  tensor  of  an  inviscid  fluid  in  terms  of  pressure. 
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The  vector  spherical  wave  functions: 

^w=Vx(r>l((r)) 

=-rVx^ 

kT 

are  the  analogues  of  the  scalar  spherical  wave  functions.  The  functions 


(24.24a) 

(24.24b) 

(24.24c) 


(24.25a) 


VLV)  =  hf\  k^Y^cp)  (24.25b) 

are  the  scalar  spherical  harmonics  for  the  longitudinal  and  transverse  components  of  the  acoustic 
field.  They  are  solutions  of  the  following  scalar  Helmholtz  Equations: 

V*V^L+  =  0  (24.26a) 

V • ' V  WL, + =  0  (24.26b) 

The  first  two  vector  spherical  wave  functions  are  solutions  of  the  transverse  wave  equation,  and 
the  third  vector  spherical  wave  function  is  a  solution  of  the  longitudinal  wave  equation. 

The  incident  and  scattered  fields  may  be  expanded  in  terms  of  the  regular  and  outgoing 
spherical  wave  functions,  respectively: 

Wtnciden,  =  £  Re  ¥cml  (24.27) 

ami 

¥  Scattered  =Y,b am, ¥mt  (24.28) 

ami 


The  displacements  and  surface  tractions  interior  to  the  target  may  be  expanded  in  terms  of 
the  regular  vector  spherical  wave  functions: 

"-=ZSAr.o»/ReFriOM/  (24.29) 

t-\  ami 

3 

L-ZIA  rm,t{  Rerrow/)  (24.30) 

x- 1  ami 
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The  displacements  and  surface  tractions  on  the  fluid  side  of  the  interface  may  be 
expressed  in  the  following  form: 


^  +  2  ^  (¥ Incident  ¥ Scattered  ) 

py 


(24.31) 


Incident  +V Scattered)* 


(24.32) 


The  displacements  and  surface  tractions  on  the  fluid  and  elastic  solid  sides  of  the 
boundary  are  related  by  the  following  boundary  conditions: 


w+  •  «  =  w_  •  n 


(24.33) 


t+  •  h  =  t_  •  n 


(24.34) 


t_xn  =  0 


(24.35) 


Substituting  Equations  24.31  and  24.32  into  Equations  24.33  and  24.34,  one  obtains  the 
following  expressions  for  the  acoustic  field  and  its  normal  gradient  at  the  boundary  in  terms  of 
the  elastic  field  in  the  interior  of  the  target: 

y/\s=7+  •h  =  F  •fi  (24.36) 

n»Vi//\s=  p+ct)2u+»fi  =  p+0)2u_»n  (24.37) 


The  scattered  field  satisfies  the  following  integral  equation: 


| dA  •  {G(r,F:  k)V  y/(r')-VG(r,r':  k)y/(F)}  = 


+  ¥  Scattered  V 

-'/'incident  (rift  V 


(24.38) 


Substituting  Equations  24.36  and  24.37  into  Equation  24.38,  one  obtains  the  following 
integral  equation  for  the  scattered  field: 


jdA{G(P,r’:  k)p+(02u_  •n-fi*VG{r,r':  k)t_  •«}  = 

5 


+  ¥ Scattered  (f  If  i  V 

~¥ Incident  (fife  V 


(24.39) 


Substituting  the  bilinear  expansion  of  the  free  field  Greens’  Function 
G(r,  r':k)  =  ik £  y/ml  (F, )  Re  y/^,  (rK )  (24-40) 

ami 

and  the  expansions  for  the  normal  displacements  and  surface  traction  given  by  Equations  24.29 
and  24.30  into  Equation  24.39  evaluated  at  a  point  exterior  to  the  scatterer,  one  obtains  the 
following  linear  equation  for  the  scattered  field: 
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boml  R-e  Qamt Xc'm'l' .a' m'V  (24.41) 

r'=l  <x'm7' 

Refi-f.rwr  =*  jJ^{p+6>2  Re^,«.ReKr.ffV/. -«•  VRe^n-^ReK^.,.)}  (24.42) 


Evaluating  Equation  24.39  at  a  point  interior  to  the  scatterer,  one  obtains  the  following 
linear  equation  for  the  incident  field: 


@om!  1  ^  Qoml  j'a’mT  ^ T\o'm'l ' 

T'=l  a'm'l' 

&-.««•  =  *  V„,»  •  R«Kt  w  -  «  •  V  •  ?(Re  Kw,.)( 

s 

The  interior  fields  satisfy  Betti’s  Third  Identity: 

JJt£4{Re  VT  aml  • t_-u_  •  F(Re  Fr  om, )}  =  0 

s 

Using  the  identity 

m_  =  hx(hxu_)-(u_  •  n)n  =  hx(nxiu_)-{u+  •«)« 
one  can  express  Equation  23.45  in  the  following  form: 
jjdA{ReVr<m/  •t_  - («  •  w+ )(« •  / (Re VTCml ))+(«x(«x u_ )) •  1  (Re VTgml)}  =  0 


(24.43) 

(24.44) 


(24.45) 


(24.46) 


(24.47) 


Assume  the  acoustic  field  on  the  surface  of  the  target  has  the  following  expansion  in 
terms  of  the  regular  spherical  wave  functions: 

^=5X,Re^,  (24.48) 

ami 


In  this  case,  the  displacements  on  the  fluid  side  of  the  interface  have  the  following  expansion: 

«♦  =  — (24.49) 

P+CO  ami 


Substituting  the  above  expansion  into  Equation  24.47,  one  arrives  at  the  following  linear 
equation  between  the  expansion  coefficients  of  the  fields  on  the  fluid  and  target  sides  of  the 
interface: 
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E  Wftw,  - 1 1  *««r4wr  =  0  <24'50> 

cr'mV  r-1  ff'mV' 

P„w,  =  /|<i({(n-F(ReK„,)X«*VRert,..r)l  <24-51> 

5 

*„w»r  =  Ij^lReRnw,,  •F(ReRr-,T-B,r)-^(Re^n»i/)*(”x(”xRe^T’1T'*T))}  (24.52) 

Combining  Equations  24.41, 24.43,  and  24.50,  one  obtains  the  following  representation 
for  the  T-matrix: 


T  =  ReQR~'P(QR~'Py 1 


(24.53) 


The  expansion  coefficients  for  the  interior  fields  are  given  by  the  following  linear 
equation  in  terms  of  the  expansion  coefficients  of  the  incident  field: 

A  =  R-'P(QR~'Py'a  (24.54) 


An  alternative  formulation  of  the  T-matrix  for  a  solid  target  is  obtained  by  using  the 
boundary  conditions  to  express  Equations  24.39  and  24.45  in  the  following  form: 


^dA{G{r, r'\  k ) p+(02u_  •  n  -  h  •  V G(r ,  r ':  kjt+  •n}  = 
s 


+  VscallereAr),?e  V 

-V Incidents  *  V 


(24.55) 


\\dA{RtVr<am,  *t+  -u_*K Ref7>om/)}  =  0  (24.56) 


Suppose  the  surface  traction  at  the  surface  of  the  target  has  the  following  expansion  in 
terms  of  the  spherical  harmonics: 

Z/vn-  (24*57) 

ami 


Evaluating  Equation  24.55  at  a  point  exterior  to  the  target,  one  obtains  the  following 
linear  equation  for  the  scattered  field: 


l^aml  Qcmt ,za'mT^T,a'm'l'  Re  M aml.o'm'l' Po'mT 

r= 1  c'm'V 

ReOm/.rvxwr  =  ikpy  ReRrW 


(24.58) 

(24.59) 
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Re  Mon,i.*mr  = ik  \\dAh  •  V  Re  (24.60) 

Evaluating  Equation  24.55  at  a  point  interior  to  the  target  one  obtains  the  following  linear 
equation  for  the  incident  field: 


3 


& Oml  ^  \Q<Pnl.m*mT^T.g%mT  ^ rani  n'm'r  ft fr'm'P 

T~]  o'm'V 

(24.61) 

Qamixo'm'v  =  tip,®2  \\dA  ¥anil  ReVwr 
s 

(24.62) 

Mmi.cwr  = ik  \\dAn  •  V  yJY^ 

s 

(24.63) 

Evaluating  Equation  24.56,  one  obtains  the  following  linear  equation  between  the  interior 
and  exterior  fields: 

3 

^TOrnt^'m'I'P  o' ^TOml  ,T'o'mT^t',o'mT 
a'm'l'  r'-l  o'm'V 

(24.64) 

W-r  =  JJ^(«*ReF,„,)W 

(24.65) 

5 

(24.66) 

Combining  Equations  24.58, 24.61,  and  24.64  one  obtains  the  following  expression  for 
the  T-matrix: 

T  =  -(Re  M  -  Re  Q1C'P)(M  -  QR-'P)-' 

(24.67) 
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25.  SCATTERING  FROM  AN  ELASTIC  TARGET  IN  A  WAVEGUIDE 


This  section  uses  the  results  of  Sections  19  through  24  to  describe  the  scattering  from  an 
elastic  target  in  a  waveguide.  This  section  assumes  the  scattering  from  the  target  is  described  by 
the  spherical  T-matrix  of  the  target. 

The  free-field  scattering  from  the  target  is  of  the  following  form: 

W Scattered  W aml^’aml ,o'mT®o'mT  (25 . 1 ) 

oml  a'/wV 

The  coefficients  a ml  are  the  expansion  coefficients  of  the  incident  wave  in  terms  of 
regular  spherical  wave  functions 

!^j Incident  ^ oml  <25-2> 

ami 


In  Section  1 9,  the  following  expansion  of  the  normal  mode  contribution  to  the  waveguide 
Greens’  Function  was  derived  in  the  following  form: 


4  p(zs) 


Ho  (<ln  l^  +  AF-Ar;.  \)F(z+Az:qn)F(zs  +  Azs)  = 


Z  Z ,*•«•/•, AU  (>*>  Zs’Z'-<ln) Re  Wond,^  ,  Az, ,  A<PS ) Re y/and.Nr,.  < ^ A (p) 


(25.3) 


ami  o'mT 


The  coefficients  of  this  expansion  are  defined  as  follows: 


cos  (m(p\(7  =  e^ 

iANmFNm  (Zv  :  <ln)Bn,l,Nm  +  ANsreFN„  (Zv  :  <ln)Bml,Nsrc  }  ^ 

,  _  _  cos(w'ffl),rr'=e 

(Z  :  <In)Bn,l,N„.  +  ANn.FNm  (Z  1  ?n)Bt mlji„  }  {sin(m< p)>£j-’=  J 


(25.4) 


Here  (ps  is  the  angle  between  the  x-axis  of  the  local  coordinate  system  at  the  source  and  the 
location  of  the  receiver  as  seen  from  the  source,  and  (p  is  the  angle  between  the  x-axis  of  the 
local  coordinate  system  at  the  receiver  and  the  location  of  the  source  as  seen  from  the  receiver. 
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In  Section  20,  the  following  expansion  of  the  directivity  functions  of  the  source  and 
receiver  in  terms  of  spherical  harmonics  was  introduced: 

In  n 

dsrc.^  =  \d<p  Jsin  ( (0,  <p)DSn  (0,  <p)  (25.5) 

0  0 

In  n 

dRcv,omt  =  jd<P  Jsin  (i?)d#Yml  (0,  <p)DRcv  (0,  (p)  (25.6) 

0  0 

Neglecting  multiple  scattering  terms,  the  normal  mode  contribution  of  the  scattering  from 
an  elastic  target  has  the  form: 


Scattered '  4)  = 

S  S  /C  dSrc,am!^aml,Nw  \a”m"l",NTtl  ^Src-Tgl  ’  ZSrc  >  ZTgl  >  *Psrc-Tgt  >  *pTgt-Src  -4) 

ami  a'mT  a"'muTu 

To'm'l\o''m''V'^G'm'l\NTRl  \a'''m'''r\NRc,.  (^Rcv-Tgt^Tgt^Rcv^Tgt-Rcv^Rcv-Tgt  *  9  ^ Rev, 


(25.7) 


The  addition  of  multiple  scattering  terms  to  the  calculations  replaces  the  free  field 
T-matrix  in  the  above  expression  with  the  waveguide  T-matrix  of  the  following  form: 

Tw  =(T~'  +iR)~'  (25.8) 

Here,  the  matrix  R  is  the  rescattering  matrix  of  References  10  through  12.  Generally,  the 
rescattering  matrix  only  affects  the  scattered  field  by  a  few  decibels. 
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26.  NORMAL  MODE  REPRESENTATION  OF  SCATTERING  FROM  A  ROUGH 

INTERFACE 


This  section  presents  the  normal  mode  description  of  the  scattering  from  a  rough 
interface.  This  formalism  can  be  utilized  to  describe  both  surface  and  bottom  reverberation  from 
a  rough  interface.  The  topic  of  bottom  reverberation  due  to  volume  inhomogeneities  in  the 
sediment  will  be  addressed  in  a  separate  section  on  the  normal  mode  representation  of  volume 
reverberation. 

This  section  begins  by  recalling  Equation  19.29  representing  the  expansion  of  the  normal 
mode  contribution  to  the  Greens’  Function  in  terms  of  a  spherical  basis  about  the  source: 

,  2m  ,  cos(m(ps),o  =  e 

Aon,i.Njr,zs,z,<ps:qn)  p(z  ^  %Jt  H">  N^Z ‘qn\m{m(ps),a  =  o  (26.1) 

F,  (zs .  q„  +  AN  Fn*"  (zs  .  qn  )Bml_Nm } 


Expand  the  normal  mode  contribution  to  the  Greens’  Function  in  terms  of  a  cylindrical 
basis  set  about  the  field  point: 


A^Jlr  +  Arlz^z  +  Az,^  :qn)  = 

£ £  5  (r, zs , z, (ps , (pr  :qn) Re/  (Ar,Az,A«p) 


(26.2) 


One  may  make  use  of  the  following  expansion  of  the  Hankel  Function: 


(26.3) 


w  =  Ju2  +u'2-2uu'cos((p) 


(26.4) 


One  may  also  make  use  of  the  following  representation  of  the  depth  function  in  terms  of 
upward  and  downward-going  plane  waves: 


Fn  (z  +  Az :  q)  =  A+  (q)Fj  (z :  q)  exp(+ih„Az)  +  A~  ( q)F„  (z :  q)  exp(-ihNAz) 


(26.5) 


Substituting  Equations  26.3  and  26.5  into  Equation  26.1,  one  obtains  the  following 
expression  for  the  expansion  coefficients  in  Equation  26.2. 
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BS oml,N>rc;o-m\N(r ,Z%,Z,<Ps,<Pr  q„  )  = 

+  (-1  THtAlS)) 

^N,rc(Zs  •  ^n)^ml,Nlrc  +  ^Nm.  ^Nm  (Zv  •  <Jn)^ml.Nm^Nr„^N,cv(Z  •  Qn  ) 

cos(m<ps),cr  =  cos(m<pr),<j'= 
sin(w0>1),<7  =  o  sm(w^r),<T'=  o* 


(26.6) 


Here  ^  is  the  angle  between  the  x-axis  of  the  local  coordinate  system  at  the  source  and  the 
location  of  the  receiver  as  seen  from  the  source,  and  <pr  is  the  angle  between  the  x-axis  of  the 
local  coordinate  system  at  the  receiver  and  the  location  of  the  source  as  seen  from  the  receiver. 

One  may  introduce  the  following  expansion  of  the  directivity  functions  of  the  source  and 
receiver  in  terms  of  spherical  harmonics: 

2  7t  It 

d src,omi  =  \d(p  Jsin {p)d^Yml  (tf,  (p)DSrc  (0,  <p)  (26.7) 

0  0 

2  k  it 

d =  jd<p  Jsin  (&)d#Yml  (tf,  <p)DRcv  (tf,  <p)  (26.8) 

0  0 

The  above  expansion  is  used  to  include  the  effects  of  the  beam  pattern  of  the  source  and  receiver 
on  the  scattered  field. 


Suppose  scattering  from  the  interface  is  described  statistically  by  a  scattering  function 
S{d,d',(p)  dependent  upon  the  frequency  and  the  incident  and  scattered  angles.  The  scattering 
from  the  rough  interface  is  equal  to  the  surface  integral  of  the  product  of  the  incident  field,  the 
scattered  field,  and  the  scattering  function  over  the  rough  surface.  Decompose  the  rough 
interface  into  rectangular  facets,  whose  areas  are  small  in  comparison  to  the  range  to  and  from 
the  facet.  In  this  case,  the  scattering  function  is  approximately  constant  over  the  area  of  the  facet, 
and  the  surface  integral  over  the  facet  requires  the  integration  of  the  following  product  of 
cylindrical  wave  functions: 

+1/2  +hl  2 

jdx  \dy Re zL (r,0, <p:q)Re xL- (r,0, q>:q')  (26.9) 

-// 2  -h!  2 


In  general,  the  above  integral  is  restricted  to  one  over  either  the  upward  or 
downward-going  components  of  the  cylindrical  wave  functions.  Equation  26.9  may  be  expressed 
as  follows: 


c  ,  ^£(m)e{m')  +"f  +hlr2  cos (m(p),a  =  e  cos(m'(p),o'=e 

Soni.a'Ar-‘l')  =  2L—r- -  J  dx  J  dyJm{rq)Ja,{rq'){  .  v  r  ’  }{  .  V  V)' 

A7r  "  «  cm/  vm/AI  /r  —  /\  rm/  ivi'/At  _ 


sin  (mq>),cr  =  o  sin  {m'(p),c'=o 


}  (26.10) 
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One  may  simplify  the  calculation  by  replacing  the  rectangular  facet  by  a  circular  facet  of 
the  same  area.  In  this  case,  Equation  26.10  takes  on  the  following  simplified  form: 


s„.„  W<« :  <t\ -  7WKW.  (?» 


(26.11) 


In  the  limit  R  -»  +°°  ,  the  integral  approaches  the  following  limit: 
linw  Smrm.(g :  q') = -Uf  -  ?') 

Here,  one  has  used  the  fact  the  following  integral  is  equal  to  the  Dirac  delta  distribution: 
'\xdxJm{qx)JMx)  =  -S(.q-q') 

J  n 


(26.12) 


(26.13) 


In  order  to  simplify  Equation  26.1 1  and  obtain  a  closed  form  solution,  replace  the  above 
integral  by  the  following  Gaussian  weighted  integral: 


\  K sf\rdrJ,  (qr)Jm  (?V)  ex^-a  V ) 


(26.14) 


The  parameter  “a  ”  is  chosen  such  that  the  following  integral  is  equal  to  the  area  of  the  facet: 


2 n  jrdr  exp (~a2r2 )  = — a  2  =JtR2  =lh 


(26.15) 


4ijtR2  4m 


(26.16) 


One  may  regard  this  approximation  as  replacing  the  incident  field  by  a  Gaussian  beam 
with  the  appropriate  beam  width. 


The  following  integral  is  found  on  page  718  of  Reference  18: 


00  1 
J xdx  exp  (~a2x2  )Jm  ( OCx)Jm  (fix)  =  — —  exp(- 


4a2  m  2a2 


(26.17) 


Substituting  this  integral  into  Equation  26.14,  one  arrives  at  the  following  approximation 
of  Equation  26.11. 


■  4)  “  W  exP(-(<72  +  q'2  W l2)Im(qq'KR2) 


(26.18) 
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The  integral  given  in  Equation  26.18  is  proportional  to  the  area  of  the  facet.  The  ratio  of 
Equation  26. 1 8  and  the  area  of  the  facet  may  be  regarded  as  the  directivity  function  of  a  flat  facet. 
The  normal  mode  contribution  to  the  scattering  from  the  facet  is  given  by  the  following 
expression: 


V* Scattered » Kb  •  ^l^Src,oml^oml,NSK;<T,m,,b/b(rsb’Zs>Zb’(Psb’(Phs  •  #) 

r\^b)  oml  s=+  o'm' 

)^a"mmr,N „„\a'm\blb  (rbr  >Zr’Zb>  Prh  >  <Phr  ‘‘  <l')d  Rev., 


ami 


(26.19) 


where  rsh  is  the  horizontal  range  from  the  source  to  the  facet,  rhr  is  the  horizontal  range  from  the 
facet  to  the  receiver,  zs  is  the  depth  of  the  source,  zr  is  the  depth  of  the  receiver,  and  zh  is  the 
depth  of  the  facet.  Other  variables  are  the  angle  between  the  x-axis  at  the  source  and  the  location 
of  the  facet  ( (psh ),  the  angle  between  the  x-axis  at  the  facet  and  the  location  of  the  source  ( (phs ), 

the  angle  between  the  x-axis  at  the  receiver  and  the  location  of  the  facet  {(prh ),  and  the  angle 
between  the  x-axis  at  the  facet  and  the  location  of  the  receiver  ( <phr ). 


The  function  Ssa1m,a.m.(a,a\p :  q,q')  is  given  by  the  following  expression  in  terms  of  the 
scattering  function  S(d,ti\(p)  and  the  function  S^^q-.q')  defined  in  Equation  26.18. 


ss£a.m\a,a\p:  q,q')  =  ■  ?') 


S(+a,-a',fi),s  =  +,s'=  + 
S(+a,+a',/}),s  =  +,s'=- 
S(-a-a\p),s  =  -,s'=  + 
[S  (~a,+a'  ,P),s  =  -,s'=- 


(26.20) 


The  angles  a ,  a'  are  the  grazing  angles  of  the  incident  and  scattered  field  at  the  interface 
defined  by  the  following  relationship,  and  the  angle  P  is  the  azimuthal  angle  between  the  source 
and  receiver  as  seen  from  the  facet: 

q  =  kNhCos(a)  (26.21a) 

q'=  kNh  cos(ar’)  (26.21b) 

In  the  case  one  is  given  a  statistical  representation  of  the  scattering  from  a  randomly 
rough  facet  in  terms  of  a  scattering  strength,  Equation  26.20  is  replaced  by  the  following 
function,  where  the  function  S(+a,+ct',P)  is  the  square  root  of  the  scattering  strength: 
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f»^j' 

^  om,G'm' 


(a,OC',p:q,q') 


=  S 


Random 

om.o'm' 


(q  ■  q') 


S(+a-a',p),s  =  +,s'=  + 
S(+a,+a',J3),s  =  +,s'=  - 
]S(-a-a\fi),s  =  -,s'=  + 


[S(-a,+a' ,  J3),s  =  s' =- 


(26.22) 


O* : «’)  = 

4nR 


8aa8mm  exp (-(q2  +  q'2  )7tR2  /  2)1  m  (qq'xR2  )^/-21og(Ml)  exp(+2;hu2) 


(26.23) 


Here  w,  and  w2  are  a  pair  of  uniformly  distributed  random  numbers  in  the  interval  from  0  to  1 
used  to  generate  a  single  realization  of  a  randomly  rough  surface. 
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27.  NORMAL  MODE  REPRESENTATION  OF  SCATTERING  FROM  VOLUME 

INHOMOGENEITIES 


This  section  presents  the  normal  mode  description  of  the  scattering  from  volume 
inhomogeneities.  This  formalism  can  be  utilized  to  describe  scattering  from  volume 
inhomogeneities  in  the  case  of  both  volume  and  bottom  reverberation.  For  simplicity,  attention 
will  be  restricted  to  the  case  of  isotropic  scattering  from  volume  inhomogeneities. 

This  section  begins  by  recalling  Equation  1 9.29  representing  the  expansion  of  the  normal 
mode  contribution  to  the  Greens’  Function  in  terms  of  a  spherical  basis  about  the  source: 


Aml,Nsn(r,zs,z,<ps:qn)  = 
(zs  '■  9n)BmI.N5re 


2m 


(q- 


P(z,) 

+  ANsrr  F,  (Zs  :  qn)Bml,Nm  } 


r)FN(z:qn){ 


cos(m<ps),<T  =  e 
sin  (m<ps),<r  =  o 


(27.1) 


Express  Equation  27.1  in  terms  of  upward  and  downward-going  plane  waves  about  the 
field  point: 

Aaml,Nsrc  (riZs>Z’(Ps'-Cln)~  q*]  ' 

A+ami,Nm.N(r,zs,<ps :  q„ ) exp(+ihN (z-zN))+  ir,zs,(ps :  q„)e\p(-ihN(z-zN)) 


J+  ,  ,  2m  fani)  „(1)/  ,cos (m(ps),<j  =  e 

(A^F^(zs:q„)B^Nm  +  A^F-Jzs:q„)B^„JA±»(qn) 


(27.3) 


Here,  one  makes  use  of  the  following  expansion  of  the  depth  function  in  terms  of  the  propagator 
matrix: 

Fn ( z:q)  =  A+n (q)exp(+ihN (z  -zN))  +  A~  (q) exp (~ihN (z - zN ))  (27.4) 


Introduce  the  following  expansion  of  the  directivity  functions  of  the  source  and  receiver 
in  terms  of  spherical  harmonics: 
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2n  it 

dsrc.am,  =  j d<P  (0,  <p)DSn.  (0,  <p)  (27.5) 

0  0 

2/r  n 

d Rev, ami  ~  \d(p  jsinCtf)^}^, (tf, (p)DRcv(P (p)  (27.6) 

0  0 

Use  the  above  expansion  to  include  the  effects  of  the  beam  pattern  of  the  source  and 
receiver  on  the  scattered  field. 


Consider  the  isotropic  scattering  due  to  volume  inhomogeneities  from  a  given  volume 
element  whose  horizontal  cross  section  is  constant.  Scattering  from  such  volume 
inhomogeneities  can  be  approximated  by  the  following  integral  over  depth: 


■“  -  z+d  /  \ 

^ijLi^jdsrc.omi  d Rcv  a.m,rArea  j  dz-———SyolumeAgmi  N  (rs,zs,z,<ps  .q)Ac.m,i,N  (rr,zr,z,(pr 

j=+  $'=+  ami  o'mT  .  P\Z) 


(27.7) 

Here,  Area  is  the  area  of  the  horizontal  cross-section  of  the  volume  element  and  SVolume  is  the 

linear  scattering  strength  per  unit  volume,  that  is,  the  scattered  field  at  unit  distance  is  equal  to 
the  product  of  the  linear  scattering  strength,  the  volume  of  the  scattering  volume,  and  the  incident 
field.  Assume  the  scattering  strength  is  constant  and  isotropic  within  each  layer  of  the 
waveguide.  In  this  case,  the  integral  over  depth  in  each  of  the  layers  in  the  waveguide  contains 
an  integral  of  the  following  form: 


jdz  exp[+ishN  (z  z  ^ )  +  is  h  N  (z  z ^ )]  — 


1 


i(shN+s'h'N) 


{exp [+ishN(zN+]  -zN)+is'h'N  (zN+]  -z„)]-l} 


(27.8) 


In  the  case  the  denominator  i(sh  +  s'h')  vanishes,  this  integral  is  equal  to  the  depth  of  the  layer 
dN  =(zn+i  ~zn) •  The  depth  integrated  scattering  from  volume  inhomogeneities  with  the  given 
horizontal  cross-sectional  area  may  therefore  be  expressed  in  terms  of  the  following  summation: 


{exp[(+/5^  +is'h'N  )(zN+]  -z„)]-l} 
(+ishN  +  is'h'N  ) 

Svolume.NJ^reaj4oml.Ns„,N(t's’Zs’<Ps  *  Q)^a'mT,N,n,N  Or  >  Zr  ’Vr  "  Q  ) 


d Src.omld Rcv.a'mT 

N  s=+  s'=+  ami  a'm'l' 

P^r) 

Pn 


(27.9) 


Here,  pN  is  the  density  of  the  N’th  layer,  and  SVo,ume  N  is  the  linear  scattering  strength  of  the  N’th 
layer. 
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Consider  the  case  were  the  scattering  from  a  random  collection  of  volume 
inhomogeneities  is  described  by  a  depth  dependent  volume  scattering  strength.  The  energy  in  the 
scattered  field  at  a  unit  distance  is  equal  to  the  product  of  the  volume  scattering  strength,  the 
volume  of  the  volume  containing  the  scatterers,  and  the  energy  of  the  incident  signal.  In  this 
case,  replace  the  product  of  the  linear  scattering  strength  and  the  cross-sectional  area  of  the 
volume  in  Equation  27.9  by  the  following  expression: 

Svoiun,e,NArea  ^^"^-^-210 g(», )  exp(+2 Jtiu2)  (27.10) 

Here,  VSSN  is  the  volume  scattering  strength  in  the  N’th  layer,  dN  is  the  thickness  of  the  N’th 
layer,  and  w,  and  u2  are  a  pair  of  uniformly  distributed  random  variables  between  0  and  1 .  In 
this  manner,  one  produces  a  single  realization  of  the  scattering  from  random  volume 
inhomogeneities.  Equation  27.10  describes  a  zero  mean  Gaussian  random  process,  with  standard 
deviation  squared  proportional  the  product  of  the  volume  and  the  volume  scattering  strength. 
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28.  TIME  DOMAIN  REPRESENTATION  OF  SCATTERING  FROM  A  ROUGH 
INTERFACE  AND  VOLUME  INHOMOGENEITIES 


This  section  combines  the  results  of  Sections  15, 26  and  27  to  create  a  time  domain 
representation  of  the  normal  mode  contribution  to  the  scattering  from  a  rough  interface  and 
volume  inhomogeneities.  Our  attention  will  be  limited  to  the  propagation  of  a  simple  Gaussian 
pulse  with  the  following  frequency  spectrum: 


S(o))  =  exp[- 


2A  O)2  1 


(28.1) 


Here,  Q)0  =  27tfQ  is  the  angular  frequency  of  the  center  frequency,  and  Act)  is  equal  to  2 n  times 
the  bandwidth  of  the  incident  signal.  The  time  domain  representation  of  the  incident  signal  is  of 
the  following  form: 


+*> 

s(t )  =  J<&aS(ty)exp[-/fUf]  =  ^^27^Act)e\p[-t2Act)2  /  2  -  ico0t] 


(28.2) 


The  approach  adopted  in  this  section  is  to  approximate  the  scattered  field  by  a  Gaussian 
integral  of  the  following  form: 


+oo 

^ Scattered (0  =  \dCt)F  {CO)eX^[~a2  ((O  ~  0)Q  f  +  fi(0)~ fi>0)] 


(28.3) 


One  approximates  this  integral  by  either  the  saddle  point  approximation  or  the  method  of 
stationary  phase.  In  the  case  of  a  broadband  signal,  the  bandwidth  is  sub-divided  into  a 
collection  of  sub-bands  for  which  the  above  technique  is  applied: 


400  'J~7T 

sScattered  (t)  -  F(0)0 )  f dm  exp  [-a2  (at -  co0  )2  +  p(co -  co0  )]  =  F(co0)— exp[+(/?  /2a)2]  (28.4) 

oc 

— oo 

This  technique  is  extended  to  an  arbitrary  band-limited  signal  by  convolving  the  scattered 
signal  from  a  broadband  Gaussian  pulse  with  the  incident  signal.  This  technique  approximates 
the  spectrum  of  the  incident  signal  by  the  following  Gaussian  weighted  spectmm,  where  the 
bandwidth  of  the  Gaussian  pulse  is  sufficiently  large  this  approximation  has  adequate  accuracy: 

S(a>)  -  S(ft>)exp[-(a>-6>0)2  /2A  co2]  (28.5) 
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From  Equations  26.6,  26.18, 26.19,  and  26.20  one  obtains  the  following  representation  of 
the  scattering  of  a  monochromatic  wave  by  a  facet  on  a  rough  interface,  where  the  function 
S(tf,  Td\(p)  is  the  bistatic  scattering  function  of  the  rough  surface: 


VscaUerAVrb  '  = 

PiA 

p(zb)  otnl  j=+  a,m’a"mV=+<T,,,m"7l 
?s,s' 


ZZZZI  2X.  ,  ami  ^Cfml ,  N  w ;  &  m Nb  (r5h,zs,zh,(psh,<phs :  q) 


Sa'm\o"m"(a>a  iP  :  QiQ  (hr  >Zr’Zh>  Wrb >  Vbr  '■ 


Rcv,a"'m"'l" 


B  am!, N„\o'm',N(r ,Zs,Z,(ps,(pr  qn)  — 


m 


P(Zs) 


j£(m)£(m')  {H("+m.{qnr)  +  H^_Jqnr)} 


(Zs  :  <ln)Kl,Nm  +  AN„FN„  (ZS  ••  9n)Kl,Nm  )AN„FNm  (Z  '  ?„) 

cos(m(ps),G  =  cos{ m<pr),o'=  e 
si n(m(ps),cr  =  o  sin (nupr),o'=o 


s^.m.(a,a\ P :  q,q')  =  S^iq  :  q')\ 


S(+a-a',p),s  =  +,s'=  + 
S(+a,+a',j3),s  =  +,s'=  - 
S(-a-a\fi),s  =  -,s'=  + 
{S(-a,+a',p),s  =  -,s'=- 


jtR2 

SowAq :  q')=—8%8*  exp (-(q2  +  q'2  )nR2  /2)Im(qq'nR2) 


(28.6) 


(28.7) 


(28.8) 


(28.9) 


The  problem  of  deriving  the  time  domain  representation  of  the  scattering  from  a  rough 
surface  entails  expressing  the  following  integral  in  the  form  of  Equation  28.3: 


¥ Scone Arsh,rrh,:  q„,q„. :  t )  = 

+  O0 

|(iryexp[-/6Jt]exp[-(fi)-r»0)2/2A6)2]^Wm,(/(FtA,^,:  qn,qn. :  co) 


(28.10) 


One  may  begin  by  deriving  an  asymptotic  expression  for  the  coefficients  defined  in 
Equation  28.7  by  substituting  the  following  asymptotic  expression  for  the  Hankel  Functions: 


~J—  exp[ix - imn 1 2- in / 4] 
V^c 


(28.11) 


Upon  substitution  of  Equation  28.1 1  into  Equation  28.7,  one  obtains  the  following 
asymptotic  expression: 
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Zs > Z' <Ps ’ tyr  •  Qn) 


exp [iqnr -in  1 4}BJ  N^.a.mJzs,z,(ps,(pr :  q„) 


I  ms 


Bs™w„v'n,',Njzs,z,<ps,<pr  :qn)  =  -^-(-/)'"+m'V^(w)f(m') 

P(Zs) 

Jn^N„  (Zs  •  Qn)Bml,Ntrc  +  ^N„Fn„  (Zs  ’  Bn)^ml,N!n:  (Z  *  #«) 

^cos(m^),<T  =  ^cos(w^r),cr'= 
sin(m^),<7  =  o  sin(/w^r),<r'=o 


(28.12) 


(28.13) 


Next,  one  derives  an  asymptotic  expression  for  the  directivity  function  of  the  facet  based 
on  the  asymptotic  expansion  for  the  modified  Bessel  Function  in  Equation  28.9.  Make  the 
following  approximation  for  Equation  28.9: 


:  ?’)  “ 


exp[-(?  -  q')27tR 2  /2]exp[-^'^R2  ]Im  (qq'nR2) 


(28.14) 


Define  the  following  asymptotic  expansion  of  Equation  28.8: 
8'X,a'm'^Z,cc' ,  0  •  Q’Q')  ~  exp[— iq~ q  )  / 2]S.m  a,m<(cc,(X , 0 .  q,q  ) 


S"Jt>,c(,P:q,q')  = 


7rR 


8*8*  exp[-qq'nR2]Im(qq'xR2) 


\S(+a,-a',0),s  =  +,s'=  + 
S(+a,+a',0),s  =  +,s'=  - 
S(-a,-a',0),s  =  -,s'=  + 
S(-a,+a',0),s  =  -,s'=  - 


(28.15) 


(28.16) 


Substituting  Equations  28.12  and  28.15  into  Equation  28.6,  one  obtains  the  following 
asymptotic  expansion  for  the  scattering  from  a  facet  on  a  rough  interface: 


V  Scatter  Jr shXrh  ‘  ^')  = 


—J—r - exp [i(qrsh+q'rrh)-{q-q')2KR2  f2]y/ScallerJrsh,rrh :  q,q ') 

m\qq  rshrrb 


V scone  Jr sh,rrh  ■  q,q ')  = 

P\^b)  Gml  .«=+  {T'm'<TwmV=+<y"'mM7" 

So'm\eT"m"(a’a'>0  •  Q’<l<)Ba"’m'"r",Nl"V;o"m'',Nt(Zr>Zb’(Prh’<Pbr  •  Q  Wtfcv.i 


;o" 


\a'm\Nh  (zs’zb’<Psh’<Pbs  ■  q) 


( y'"w"7" 


(28.17) 


(28.18) 
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Substituting  Equation  28.17  for  the  asymptotic  expression  for  the  scattered  field  from  a 
monochromatic  source  into  Equation  28.10  one  obtains  the  following  expression: 

Vsca«erejl?sh,rrh,:qn,qn.  :t)  = 

+oo 

J^yexp[-/^]exp[-(<y-6>0)2 /2&a)2  +  Kq„rsh+q„.rrh)-(qn -qn.)27rR2 12}  (28.19) 


•  -I  I  V Scattered  (^tA  »  ^ rh  ^q  n  ’  q  n'  •  ®) 

m  V  <lnqn'rshrrh 

Expand  the  integral  about  the  center  frequency  by  making  the  following  substitutions  and 
keeping  only  terms  of  quadratic  order  in  the  exponential  evaluated  about  the  center  frequency: 


<ln  «  q „0  +  q'n 0  (®  -  ©0  )  +  <l''nO  (O)  ~  0)Q)2  /  2 

(28.20a) 

<i i ,■  =  qn- o  +  q'n- o  ) + ?"„•<>  )2  /  2 

(28.20b) 

(28.20c) 

_dqn(0)  =  0) o) 

H  n0  “  , 

d(o 

(28.20d) 

q  n0~  ,2 

day 

(28.20e) 

V'scanered(rsh,r,h,:qn,qn.:t)  = 

2  1  1  ~ 

_•  *1  /  ^  Scattered  ^sh  »  ^rb  »  *7n0  ’  9«'0  *  ^0  ) 

**  V  <1*4*'**** 

(28.21) 

+oo 

eXP[^nor*A  +  q  n'O^rb  )  “  “  (tf„0  “  q  n'tyf  ^  /  2]  Jt/tU  eXp[-flT2  (<0  -  6>0  )2 

-oo 

-iy?(<y-©0)] 

«2  =  2Affl2  W'-O^A  +^0^a)  +  (^O~^o)2^2/2 

(28.22) 

P  =  t~  rsh  +  q'n-0  rrh  )  -  Kq„  o  -  ?„'<)  )(q'n0~q'n'0  )^2 

(28.23) 

Upon  completing  the  Gaussian  in  Equation  28.21,  one  arrives  at  the  following  expression 
for  the  field  scattered  from  the  facet  on  the  rough  interface: 
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¥ Scattered  (Kb  ’Kb’'  Qn’Qn'  :0 


1 


-Ysca„ered(Kb’rrb’<lnO’<ln' 0  :®0) 


(28.24) 


,>p2m 


^PU(Cl„Orsb+(In-Orrb)~i0}OtO  ~  (?«0  ~  ^  n'of  ^  I  2}—— GXp[~{t  ~  t  I  It*  lO)Q(t  tQ  )] 

c 

t0  =  (q\o  rsb  +q\'o  rrb)-Kqn  0  )^2 


2_n2 


f2  =  2a2  =  — — ^  -  2i(q"„0  rsh  +  q\,0  rrh )  +  (q'n0-q'n. 0  )  kR 
Ad) 


(28.25) 

(28.26) 


The  parameters  t0  and  r  are  the  complex  arrival  time  and  pulse-length  of  the  scattered  signal, 
respectively.  Neglecting  terms  proportional  to  the  area  of  the  facet,  the  travel  time  reduces  to  the 
sum  of  the  travel  times  of  the  two  modes. 

From  Equations  27.3  and  27.9,  one  obtains  the  following  expression  for  the  depth 
integrated  scattering  from  volume  inhomogeneities  with  horizontal  cross-sectional  area  Area : 


¥  Scattered  (.W  ,Z  S,Z  r,(ps  ,Qr  q  q  I  0})  = 


^  j  ^  \  ^  i  ^  ^  1  ^ Src.oml^ Rcv.a'mT 


{Gxp[(+ish„  +  is'ti  N  )(zA,+,  -Zjy)]-!} 


(+ishN  +is'h'N  ) 

^  Volume,  N  AreaAml  Ns"  N (rs , zs ,<ps  .  qn )Aa,m,VN^N (rr ,zr,(pr  .  qn- ) 


(28.27) 


N  s=+  5  — +  ami  a'm'l 
Pi£r) 

PN 


.  ,  .  2m  \£{m)uW.  cos  (m<ps),(T  =  e 

A  ami,N^N{r,zs,(ps.qn)  p(zj^  ”  (q"r){sin(m<pf),(T  =  o} 

Mal, (zs  ‘.qn)Bml^  +  A^ (zs  -q„)BmljN^}A  N(qn) 


(28.28) 


The  time  domain  representation  of  the  scattered  field  is  given  by  the  following  equation: 

¥  Scattered  0,  .  f,  „ Zs  ,Zr  ,<PS  ,<Pr  l  q„  ,  „•  it)  = 

+oo 

\dcoexp[-i(tX]exp[-((0-coQ)2 l2Aco2]y/Scatlered{rs,rr„zs,zr,(ps,(pr : qn,q„r.Q)) 


(28.29) 


Substituting  the  asymptotic  expansion  for  the  Hankel  Function  into  Equation  28.28,  one 
arrives  at  the  following  asymptotic  expansion  of  Equation  28.27: 


¥  Scattered  >  F  ’  Zs  ’Zr’VX’<Pr  ’^n’  ?  „•  •  = 

2  1 

—  - exp [i(qnrs  +qn-rr)]¥Scallered(zs,zr, (ps, (pr  :qn,qnr.O)) 

m  \  qnqn^srr 


(28.30) 
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W  Scattered  ’  ^Ps  ’  Vr  ' 

yyyyv,  d  {exp[(+/^  ±isjh\  -zN)]~  1} 

La  La  La  La  La  U Src,omla Rcvto'mT  7  •  /  ^Tm  \ 

Af  5=+  J'=+  ami  a'm'i'  (+tf  A?^  +  W  AT  ^  ) 

S volume^ AreaA(^1N^  N{zx,(ps  •  q„)Aa,mT N(zr9<pr  :qn,) 


(28.31) 


5 : *.)  ~^H>- -  *) 

p(zj  V  8;r  sm {m(ps\a  =  o 
(ALF»Jzs  :^n)K,.Nm^A-mF;jzs:qn)B-}A\{qn) 


(28.32) 


Substituting  Equation  28.30  into  Equation  28.29  and  expanding  the  argument  about  the 
center  frequency,  keeping  terms  of  quadratic  order  in  the  expansion,  one  arrives  at  the  following 
approximation  of  the  scattered  field: 


V Scattered  O's » »  ^r > >  ^Pr  •  fin  >  9n'  •  0  — 


**'  \Mn0Wri 


W Scattered  (^j  ’  ^r’^Ps  >  ^P r  •  7  nO  ">Pl n'O  •  ^0  ) 


(28.33) 


Tw 

+^orrA)-'6>0/]  |^6>exp[-a2(«-<y0)2  -//?(&>- 6;0 )] 


^  =  77  ~  +  *"«’o  ) 


(28.34) 


P  =  *-(<!' »0rsh+q'n-0rrh) 


(28.35) 


Upon  completing  the  Gaussian  in  the  integrand,  one  obtains  the  following  asymptotic 
representation  of  the  scattered  field: 

V Scattered  (?, ,  T,  „  Z, ,  Z,  ,  ^  ^  ,  0„.  :  /)  = 


®1w,0 


V Scattered  (^s  >^r  >  1 Ps  »  ^Pr  •  9n0  »  *7 n’0  •  ^0  ) 


(28.36) 


eXP[%n0ri  +  Q n'O^r  )  -  ^fl]- - eXp[-(/  -  /0  )2  /  2T2  -  /6>0  (/  -  /„  )] 


'o  +?7>'V) 


(28.37) 


r2=2«2=^r-W„o^+^o^) 


(28.38) 


The  parameters  /0  and  r  are  the  effective  arrival  time  and  pulse-length  of  the  scattered  signal. 
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29.  PLANE  WAVE  APPROXIMATION  OF  SIGNAL-TO-NOISE  RATIO 

CALCULATIONS 


This  section  uses  the  plane  wave  approximation  of  the  normal  modes  to  calculate 
signal-to-noise  ratio  (SNR).  This  approximation  decomposes  the  normal  mode  contribution  to 
the  Greens’  Function  at  the  source  and  field  point  into  upward  and  downward-going  plane  waves 
with  a  grazing  angle  given  by  the  following  relationship  between  vertical  wavenumber  and 
acoustic  wavenumber  at  the  source  and  field  points: 

Re(^c(g))  =  sin(0  )  (29. 1 ) 

kSrc 


The  parameter  &Src  denotes  the  grazing  angle  of  the  normal  mode  at  the  source.  A  similar 

expression  is  used  for  the  grazing  angle  of  the  normal  mode  at  the  field  point.  The  vertical 
wavenumber  is  positive  for  downward  going  plane  waves,  and  negative  for  upward  going  plane 
waves.  Similarly,  the  grazing  angle  is  positive  for  downward-going  plane  waves  and  negative  for 
upward-going  plane  waves.  This  sign  convention  is  related  to  the  fact  that  the  z-axis  is  directed 
downwards  in  PC  SWAT. 

One  begins  with  the  time  domain  representation  of  a  normal  mode.  Propagation  of  a 
pulse  is  described  by  the  Fourier  Transform: 

Gn(r,z,,z:t)  =  jS(Q))G„(r,zs,z:o))e~lMdO)  (29.2) 


Here  S(a>)  is  the  spectrum  of  the  incident  signal.  In  the  following  instance,  one  may  assume  a 
spectrum  of  the  following  form: 

(29-3) 

2A  nr 

Here,  co0  is  the  center  angular  frequency,  and  A  co  is  the  angular  bandwidth  of  the  pulse. 

Propagation  of  an  arbitrary  band-limited  pulse  can  be  obtained  by  convolving  the  scattered  signal 
from  a  broadband  Gaussian  pulse  with  the  incident  signal.  The  incident  signal  is  of  the  following 
form: 


+oo 


s(t)  =  J S(co)e  '°*d(Q  -  y/lizAeoexpi-t2 An)2  /  2  -  ico0t ) 

— oo 


(29.4) 


The  parameter  x  represents  the  pulse  length  of  the  incident  signal. 


29-1 


CSS/TR-02/10 


(29.5) 


Equation  15.17  describes  the  time  domain  representation  of  the  Greens’  Function,  as 
shown  below: 


Gn  (r,  zs,z:t)  -  :  )F(Z :  )  exp(/(^„r  -  co0t  -  it  /  4)  -  (t  -  t„  )2  /  4a2 ) 


Here,  the  reciprocal  of  the  group  velocity  of  the  mode  is: 


*’.(«) 


_dqn{0)) 


(29.6) 


(29.7) 


the  rate  of  change  of  the  reciprocal  of  the  group  velocity  is: 

„  ,  v  d2qAco) 
da 

the  travel  time  for  the  mode  is: 

K  =rq'n(Q) o) 

and,  the  complex  width  squared  of  the  effective  Gaussian  is: 


(29.8) 


(29.9) 


(29.10) 


The  normal  mode  terms  in  Equation  29.6  can  be  decomposed  into  an  upward  and 
downward-going  plane  wave: 


Gin'±(r,zs,z:t)  = 
i  12  ‘Jk 
4P(2s)\qn r*  « 


F±(ZS'-  q„ )F± (z'-qn) exp (i(q„r  -  o)0t  -  n / 4)  -  (/  -  tn  f  / 4a2 ) 


(29.11) 


Here,  the  function  F  (z :  q)  represents  the  decomposition  of  the  depth  function  into  upward  and 
downward-going  plane  waves,  as  defined  below  in  terms  of  the  components  of  the  propagator 
matrix  in  the  N’th  layer  containing  the  depth  z: 


F±  ( z:q)  =  A\  ( q)F *  (z:q)  =  A*  (q)  exp [±ihN  ( q)(z  -  zN  )] 


(29.12) 
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The  normal  mode  contribution  to  the  signal  from  the  target  is  of  the  following  form: 


y/s’s'’s”’s'" Scattered  (rs , rr , zt ,  Zr ,  2, <ps ,  <pr  :  qn ,  qn.  :  t  + 1')  =  (4 K)1  ASrcDSrc  (s $Sn. , K  -  <ps ) 
DRcv  (s'  "&Rcv,7t  -  <pr  )GS/  (rs,zs,z:  t)GsJ's'"  (rr,z,zr:t')ST„get(s'dn,n-(ps,s"  ,  <pr ) 


(29.13) 


The  parameter^  =  10Si/2°  is  the  amplitude  of  the  incident  field  at  a  unit  distance  from  the 
source,  where  SL  is  the  source  level  of  the  projector  in  decibels.  The  vectors  (rs  ,zs,(ps )  and  s 
(r,  ,zr,<pr)  are  the  cylindrical  coordinates  of  the  source  and  receiver  relative  to  the  target,  that  is, 
rs  is  the  range  from  the  target  to  the  source,  zs  is  the  depth  of  the  source,  and  <ps  is  the  bearing 
of  the  source  as  seen  from  the  target.  Similarly,  rr  is  the  range  from  the  target  to  the  receiver,  zr 
is  the  depth  of  the  receiver,  and  <pr  is  the  bearing  of  the  receiver  as  seen  from  the  target.  The 
functions  Gs/  (r,z,z'\  t )  are  the  decomposition  of  the  n’th  normal  mode  contribution  to  the 

Greens’  Function  in  terms  of  upward  and  downward  going  plane  waves  at  the  two  endpoints. 

The  angle: 


&n  =  sin-1  (Re(-^-^))  (29.14) 

k 

is  the  grazing  angle  of  the  n’th  normal  mode  at  the  target.  The  angles  &Src  and  tfRcv  are  the 
grazing  angles  for  the  outgoing  and  incoming  normal  mode  at  the  source  and  receiver, 
respectively.  The  function  STarget($,(p,d',<p')  is  the  bistatic  scattering  function  for  the  target.  In 

the  case  of  an  omni-directional  point  target,  this  function  is  equal  to  the  constant  10 75/20 ,  where 
TS  is  the  target  strength  in  decibels. 

In  general,  one  estimates  the  signal  from  the  target  by  incoherently  adding  the  magnitude 
squared  of  all  the  normal  mode  contributions  given  by  Equation  29.13  that  arrive  at  the  receiver 
within  a  pulse  length  of  each  other. 

The  next  step  in  evaluating  the  signal-to-noise  ratio  of  a  system  is  to  estimate  the  amount 
of  surface,  bottom,  and  volume  reverberation  arriving  at  the  receiver  as  a  function  of  travel  time. 


First,  consider  the  case  of  calculating  the  mean  square  surface  reverberation.  The 
magnitude  squared  of  the  normal  mode  contribution  to  surface  reverberation  is  of  the  following 
form: 


|  y/S'S' Surface  Sr  ,  Z,  ,Z,  , Z  =  0,  Q  q„  ,q„.  I  t  +  f)  |2  = 

(4 *)4 1  ASrc  |2|  DSrc(s0Src, 0)DRcv(s'#Rcv,<p)  |2|  G?-(r„z„z 


=  0:t)G+/(rr,z  =  0,zr:O\ 2 


S Surface  (dn-’dn-.K  ~  <P)(V  SrS* 


dqn 

dco 


A  #„orAR) 


(29.14) 


Here,  the  quantity  SSurface(d, is  the  bistatic  surface  scattering  strength. 
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The  term: 
do 

Area  =  v^r,  -^-AdHorAR  (29. 1 5) 

is  the  area  ensonified  by  the  incident  signal,  where  AR  is  the  range  resolution  of  the  signal,  and 
At? H0r  is  the  horizontal  beam- width  of  the  system.  The  parameter  rs  is  the  range  from  the 

projector  to  the  ensonified  area  on  the  surface,  and  zs  is  the  depth  of  the  source.  Similarly,  rr  is 
the  range  from  the  receiver  to  the  ensonified  area,  and  zr  is  the  depth  of  the  receiver.  The 
parameter  z  is  the  depth  of  the  surface.  The  angle  (p  is  the  azimuthal  angle  between  the  projector 

and  the  receiver,  as  seen  from  the  center  of  the  ensonified  area.  For  a  monostatic,  active  sonar, 
this  angle  is  approximately  zero. 

In  the  discussion  of  bottom  reverberation,  one  may  decompose  the  bottom  reverberation 
into  two  parts.  The  first  describes  bottom  reverberation  due  to  the  roughness  of  the  interface. 

The  second  describes  bottom  reverberation  due  to  volume  inhomogeneities  in  the  sediment. 

High  frequency  models  of  bottom  reverberation  usually  combine  these  two  contributions  into  a 
single  scattering  function,  since  the  sound  is  generally  limited  to  a  narrow  layer  about  the 
interface  owing  to  the  attenuation  of  sound  in  the  sediment.  In  the  case  of  low  frequency 
propagation,  the  sound  has  significant  penetration  into  the  sediment,  and  scattering  from  volume 
inhomogeneities  cannot  be  considered  to  arise  from  scattering  from  a  thin  layer  about  the 
interface. 

|  y/s-s' bo,, on, (rs,rr,zs,zr,zh,q>:qn,qn. :/  +  /')  |2  = 

(4*)4 1  Asrc  \2\ DSrc (s #5rc ,0)DRiv (s' dRcv , <p)  |2|  Gsn’+(rs,zs,zh :  t)G~/(rr,zh,zr :  /')  |2  (29.16) 

SBo„oM,t}n,,<pr  -<ps  +7C)(vSrcrs^Ai)„nrAR) 

act) 

Here,  the  quantity  SBotlom (i?,  is  the  bistatic  bottom  scattering  strength  for  the  rough 
interface.  The  term: 

da 

Area  =  -£-A  #H„rAR  (29. 1 7) 

is  the  area  ensonified  by  the  incident  signal,  where  A R  is  the  range  resolution  of  the  signal,  and 
A*?//or  is  the  horizontal  beam-width  of  the  system.  The  parameter  rs  is  the  range  from  the 

projector  to  the  ensonified  area,  and  zs  is  the  depth  of  the  source.  Similarly,  rr  is  the  range  from 
the  receiver  to  the  ensonified  area,  and  zr  is  the  depth  of  the  receiver.  The  angle  <p  is  the 
azimuthal  angle  between  the  projector  and  the  receiver,  as  seen  from  the  center  of  the  ensonified 
area.  The  parameter  z  is  the  depth  of  the  bottom  interface.  In  the  case  of  a  multi-layered  bottom, 
one  can  easily  sum  over  the  scattering  from  each  interface. 
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The  magnitude  squared  of  the  mean  normal  mode  contribution  due  to  the  scattering  from 
volume  inhomogeneities  in  the  N’th  layer  is  of  the  following  form: 


ll^' 


Volume 


(rt,rr,z„z„z,<p:q„q,.:t+t')\2= 


my  |  ASrc  |J|  DSrc(s#Src,0)D,Js"'A„,<P) I’l  P^rVPs  l!|  G’S (r„z,  :r)Gf'>r,z,  :C) |! 


(W, 

dco 


exp [{+is'hN(qH)  +  is'h'N  (gJHf^  zlullll 
(+is"hN(qn)  +  is"h'N  (<?„.)) 


(29.18) 

The  function  G^'±(r,zs  \t)  is  the  following  function,  where  the  coefficients  A±N{q)  are 
the  coefficients  of  the  depth  function  in  the  N’th  layer: 


G°/{r,zsU)  = 


2 


4 p(zs)Xq„rx  a 


Fs(zs  :  qn)ANs'(qn)exp(i(qnr  -  O)0t  t„f  /  4a2) 


(29.19) 


These  functions  arise  by  integrating  the  product  of  the  two  Greens’  Functions  over  the 
depth  of  the  scatterer: 


zN+\ 


jdzG„s’s’(rf,zs,z :  t)Gf's'\rr,z,zr :  t')  = 


\dzGns  s  (rs,zs,z  :  t)Gnf"'s"(rr,zr,z :  /') 

Pn  l 


(29.20) 


P(Zr) 

Pn 


G°/{rs,zs:t)G^\rr,zr:t') 


exp  [{+/5,^(gw)  +  /V,/;V  (g„.)}(^+1  -Zy)]-1 
(+is'hN(qn)  +  is"h'N  (qn,)) 


One  may  use  the  following  expression  for  the  upward  and  downward-going  components 
of  the  depth  function  in  terms  of  exponentials  and  the  coefficients  of  the  propagator  matrix  in 
evaluating  the  above  integral. 

F±(z:q)  =  A*  (q)  exp[±ihN  (q)(z  -  zN )]  (29.2 1 ) 


An  alternative  expression  for  the  scattering  from  volume  inhomogeneities  is  to  replace 
Equation  29. 1 8  by  the  following  expression: 
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|  y/sy': volume  (rs,rr,zs,zr,z,<p:qn,qn.:t  +  t')\2= 

(4*)4  I  Asrc  \1\DsrC(s#Sr,mRJs'''K,,<P)  \2\p(Zr)^N  f|  l  t)G™  [t,  ,  Z,  if)  |2 

(v  r  4-  a  AffiV  —7  vie  (1  —  exp[-(25  Im(/t)  +  25  Im(^'))(zAf+l  —  zN )]) 

V  Src*  s 


da) 


a*  ,»(.  .  „«.  (l-exp[-(2S'to(A)  +  2S"Im(A1))(z„.,-; 

L*UHorlM'-\ZN+ 1  Z N  Volume  TT  '  - “ 

(25  Im(/j)  +  25  Im(A  ))(zAr+1  -zN) 


(29.22) 

In  deriving  this  expression  one  replaces  Equation  29.20  by  the  following  incoherent  sum  of  the 
magnitude  squared  of  the  Greens’  Function  over  depth,  that  is,  one  performs  a  summation  over 
the  energy  rather  than  pressure  as  a  function  of  depth: 


r*+l 

\dz\Gn  '  (rs,zs,z :  t)Gj  '*\rr,z,zr :  t') 

l  =  f  dz 

G„  (rs,zs,z't)Gn.  (r  ,z  ,z: /') 

*S 

J 

2N 

Pn 

^G,"'(r„z,  :r)G"(r„z,  :C)| 
Pn 

P(Zr) 


2  Z.V.I 


Ps 


Gs/(rs,zs:t)Gsny(rr,zr:t') 


\dz  exp[-(2s'  Im  (h)  +  25"  Im (h'  ))(z -z„)]  = 

:S 

'  _z  (1  - exp[-(25'Im(/»)  +  25"Im(^'))(zJV+,  - z„)]) 
N+'  N  (2s']m(h)  +  2s"lm(h'))(zN+l -zN) 


(29.23) 

In  the  limit  the  imaginary  component  of  the  vertical  wavenumber  vanishes  the  following 
limit  is  obtained: 


(1  -  exp[-(25'Im(/Q  +  25"Im(/f  ))(zv+1  -  zv )]) 
(2s’lm(h)  +  2s"]m(h'))(zfUi  -zN) 


(29.24) 


The  quantity  SVo,ume  is  the  volume  scattering  strength  for  volume  inhomogeneities  in  the  N’th 
layer.  The  term: 


Volume  =  vSrcrs 


do) 


^t?WorAR(zAf+1  zN) 


(29.25) 


is  the  volume  ensonified  by  the  incident  signal.  The  parameter  r,  is  the  range  from  the  projector 
to  the  ensonified  volume,  and  zs  is  the  depth  of  the  source.  Similarly,  rr  is  the  range  from  the 
receiver  to  the  ensonified  volume,  and  zr  is  the  depth  of  the  receiver.  The  angle  (p  is  the 
azimuthal  angle  between  the  projector  and  the  receiver,  as  seen  from  the  center  of  the  ensonified 
volume.  In  the  case  of  a  multi-layered  bottom,  one  can  easily  sum  over  the  scattering  of  the 
volume  inhomogeneities  in  each  layer.  Scattering  from  volume  inhomogeneities  in  the  water 
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column  (volume  reverberation)  is  similarly  treated.  One  must  implicitly  assume  the  directivity  of 
the  beam  is  approximately  constant  over  the  depth  of  the  layer,  thus  application  of  Equation 
29.18  or  29.22  may  require  sub-dividing  the  layers  into  smaller  layers. 

Calculation  of  the  signal-to-noise  ratio  begins  with  calculating  the  return  from  the  target 
at  a  fixed  range  by  adding  those  returns  from  the  target  that  return  within  a  pulse  length  of  each 
other  using  Equation  29.13.  This  process  produces  the  envelope  of  the  time  series  representing 
the  different  multipath  returns  from  the  target.  Next  one  calculates  the  scattering  from  the 
surface  and  the  bottom  by  integrating  Equations  29.14  and  29.16  with  respect  to  the  range  from 
the  projector  to  the  ensonified  area  on  the  surface  and  bottom  respectively.  Again,  one  adds 
contributions  to  the  envelope  of  the  surface  and  bottom  reverberation  that  lie  within  a  pulse 
length  of  each  other.  This  process  produces  the  envelope  of  a  time  series  representing  the 
scattering  from  the  surface  and  the  bottom.  Next,  one  uses  either  Equation  29.18  or  29.22  to 
compute  the  scattering  from  volume  inhomogeneities  in  the  sediment  by  integrating  these 
equations  with  respect  to  the  range  from  the  projector  to  the  ensonified  volume  in  the  sediment. 
Similarly,  one  integrates  either  Equation  29.18  or  29.22  with  respect  to  range  from  the  projector 
to  the  ensonified  volume  to  obtain  the  scattering  from  volume  inhomogeneities  in  the  water 
column.  Finally,  one  adds  the  ambient  noise  term  to  the  sum  of  the  surface,  bottom,  and  volume 
reverberation  levels  to  determine  the  total  noise  in  the  system.  The  signal-to-noise  ratio  for  the 
system  can  be  estimated  by  taking  the  ratio  of  the  largest  return  from  the  target  and  the  total  noise 
arriving  at  the  sonar  at  the  same  time. 
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30.  COMPUTATION  OF  NORMAL  MODES 


Direct  application  of  the  characteristic  equations  in  Sections  8  and  9  lead  to  numerical 
instability  in  the  case  of  large  sound  speed  gradients  in  the  water  column.  This  section  describes 
a  numerically  stable  technique  for  solving  the  characteristic  equation  for  the  normal  modes  in  a 
waveguide. 


Suppose  one  has  a  waveguide  consisting  of  N  homogeneous  layers,  where  the  following 
matrices  transform  the  depth  coefficient  at  a  given  interface  into  the  pressure  and  normal 
displacement  at  that  interface: 


=exp(+/7 iNdN) 

(30.1a) 

M™'N  =cxp(-ihNdN) 

(30.1b) 

M"+',N  =  +I—exp(+ih/.,dN ) 

(30.1c) 

Pn 

w  v+i,  v  =_jKl  exp {-ihNdN ) 

(30.  Id) 

Pn 

K'N  =  +1 

(30.  le) 

M?/  =+\ 

(30.  If) 

mn,n  =+IhL 

Pn 

(30.  lg) 

mnn  =_jhL 

(30.  lh) 

Pn 

Suppose  the  depth  function  in  the  N’th  layer  is  of  the  following  form: 

Fn  (•z)  —  exp (+ihN  (z  zN ))  +  AN  exp(  ihN  (z  zN ))  (30.2) 


The  equations  of  continuity  of  pressure  and  the  normal  displacement  at  the  N’th  interface 
are  of  the  following  form: 
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0  = 


Mnn~'V 
1VI\,2 

MN,N- 
2,2 


J  + 

a-n-\ 


\  “N- 1  J 


tjN.N 

Mu 

\*n,n 

M2,i 


A 'MV  J*  ^ 

A-N 

A 


K.2 

i  jr  NtN 

M2,2  J 


\“N  J 


(30.3) 


In  the  case  of  a  rigid  bottom,  one  can  express  the  equations  of  continuity  and  the 
boundary  condition  at  the  pressure  release  surface  and  the  rigid  bottom  by  the  following  set  of 
linear  equations: 


0  =  CA 

at  =  {a;,a;,a:,a[,..a;,a~n) 


(30.4) 

(30.5) 


The  global  matrix  C  is  the  matrix  whose  first  row  is  given  by  the  following  expression  that 
encapsulates  the  boundary  condition  of  a  pressure  release  surface. 


^\,\...2N  ~  (U,0, 0,...,0,0) 


(30.6) 


The  last  row  of  the  global  matrix  C  takes  on  the  following  for  and  encapsulates  the  boundary 
condition  of  a  rigid  bottom: 


'2N.X...2N 


0,0,...,+/— exp  (+ihNd  N ),-/ — exp  {-ihN  dN ) 
Pn  PN 


(30.7) 


The  intermediary  rows  of  the  global  matrix  C  express  the  equations  of  continuity  of  pressure  and 
normal  displacement  in  the  intermediary  interfaces.  The  non-zero  elements  of  the  global  matrix 
C  in  the  intermediary  rows  are  of  the  following  form: 


^2n+i+l,2n+ j+ 1  —  +Mjj 


'~'2(/j+1  )+/+!,  2«+ y+1  ~~ 


(30.8a) 

(30.8b) 


The  global  matrix  C  is  extremely  sparse  and  has  the  following  quasi-diagonal  form. 


f  1 

i 

0 

0 

0 

0 

K° 

K 

-K 

-MU 

0 

0 

c  = 

mS 

mS 

-K 

-M" 

0 

0  .. 

0 

0 

M$ 

-MS 

,  0 

0 

MS 

MS 

-MS 

-MS 

(30.9) 


In  the  case  the  bottom  layer  is  a  homogeneous  half  space  the  last  row  in  the  global  matrix  is  of 
the  following  form. 


'2N.1...2N 


=  (0,0,.  ..,0,1) 


(30.10) 
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This  equation  encapsulates  the  boundary  condition  that  the  upward-going  component  of  the 
depth  function  in  the  homogeneous  half  space  vanishes. 

Note,  terms  with  a  negative  exponential  in  the  global  matrix  C  may  be  eliminated  by 
rescaling  the  columns  of  the  matrix.  In  this  manner,  the  off-diagonal  terms  are  either  Order  1  or 
they  are  exponentially  damped.  The  characteristic  equation  for  the  normal  modes  is  found  by 
tri-diagonalizing  the  global  matrix  C  and  taking  the  product  of  the  remaining  diagonal  terms  to 
obtain  the  determinant  of  the  global  matrix  C.  The  normal  modes  are  those  values  of  the 
horizontal  wavenumber  for  which  this  determinant  vanishes. 


PC  SWAT  7.0  solves  this  characteristic  equation  for  the  normal  modes  in  the  case  the 
sound  speed  is  real  by  using  a  combination  of  the  method  of  bisection  and  Newton’s  Method  to 
obtain  the  solution  for  the  normal  modes  in  the  absence  of  attenuation.  PC  SWAT  then  uses 
perturbation  theory  of  the  normal  modes  to  estimate  the  complex  eigenvalues  due  to  attenuation 
in  the  water  column.  The  resulting  normal  modes  have  a  small  imaginary  component 
corresponding  to  the  attenuation  of  the  mode.  Attenuation  of  the  normal  mode  in  the  sediment  is 
handled  by  altering  the  vertical  wavenumber  in  the  sediment,  which  generally  has  a  significant 
positive  imaginary  component.  The  effects  of  attenuation  in  the  sediment  on  the  horizontal 
wavenumber  are  ignored,  since  the  attenuation  is  generally  large  enough  that  perturbation  theory 
is  not  applicable. 


In  the  case  of  a  rigid  bottom,  the  coefficients  of  the  depth  functions  are  found  by  solving  a 
linear  equation  of  the  following  form: 


V  =  GB 


(30.11) 


f-K 

-MU 

0 

0 

0 

0  " 

-M'2l 

0 

0 

0 

0 

G  = 

K' 

MU  - 

~M?f 

-Mil 

0 

0  ••• 

(30.12) 

Mil 

Mil  - 

-MU 

Mil 

0 

0 

l  0 

0 

Mil 

MU 

-MU 

-MU  , 

b  =  (a;,a; 

M) 

(30.13) 

'+11 
v-  b 

1). 

(30.14a) 

-(Af10 

r+r 

-b 

)  2 

(30.14b) 

K  = 

0,  i  ^  1,2 

(30.14c) 
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The  matrix  G  is  the  2(N-1)  x  2(N-1)  dimensional  sub-matrix  of  the  global  matrix  C 
obtained  by  eliminating  the  first  two  columns  and  the  first  and  last  rows.  The  vector  V  represents 
the  source  term  due  to  the  depth  coefficients  in  the  top  layer,  where  one  must  assume  the  depth 
coefficients  in  the  top  layer  are  given  by  the  following  expression: 


fA+) 

J 

-K 

PC  SWAT  solves  the  above  equations  by  using  Gaussian  elimination  with  partial 
pivoting  to  find  the  inverse  of  the  matrix  G.  The  coefficients  of  the  depth  function  in  layers 
1,2,...N  are  then  obtained  by  the  following  matrix  multiplication: 

(•^i  >^i  =  G  V  (30.16) 

In  the  case  of  a  homogeneous  half  space,  PC  SWAT  forms  the  matrix  2(N-2)  x  2(N-2) 
dimensional  sub-matrix  of  the  matrix  C  by  eliminating  the  first  two  columns  and  the  last  two 
columns,  the  first  row,  and  the  last  three  rows.  The  vector  B  is  the  2(N-2)  vector  containing  the 
depth  coefficients  for  layers  1,2,...(N-1).  The  depth  coefficient  in  the  homogeneous  layer  is 
obtained  by  using  the  propagator  matrices  to  relate  the  depth  coefficients  in  the  (N-l)’th  layer  to 
the  N’th  layer. 
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31.  SAMPLE  CALCULATIONS 


This  section  describes  a  set  of  test  cases  of  the  low  frequency  propagation  model.  This 
model  is  described  in  Reference  1 . 


Table  31-1  displays  a  comparison  of  the  complex  Green’s  Function  for  a  rigid  waveguide 
of  depth  300  m  calculated  by  PC  SWAT  and  the  exact  normal  mode  solution  at  a  frequency  of 
200  Hz.  The  source  is  located  at  mid  water  column.  The  range  to  the  field  point  is  300  m.  As 
can  be  seen  from  Table  31-1,  PC  SWAT  agrees  with  the  exact  solution  to  six  significant  figures. 

TABLE  31-1.  COMPARISON  OF  EXACT  AND  APPROXIMATE  SOLUTION  OF  THE  GREEN’S 

FUNCTION 


PCSWAT  EXACT 


Depth,  M 

Real 

Imaginary 

Real 

Imaginary 

0 

2.156952E-21 

-1.742372E-48 

0.000000E+00 

0.000000E+00 

3.1 

-7.601606E-06 

-7.848293E-05 

-7.601 606E-06 

-7.848293E-05 

6.2 

-9.370370E-04 

1.848949E-04 

-9.370370E-04 

1.848949E-04 

9.3 

5.9371 88E-04 

5.420785E-05 

5.9371 88E-04 

5.420785E-05 

12.4 

4.005531  E-04 

-2.861 161E-04 

4.005531  E-04 

-2.861 161  E-04 

15.5 

-3.460245E-05 

3.103788E-04 

-3.460245E-05 

3.103788E-04 

18.6 

7.772559E-05 

-1.685571  E-04 

7.772559E-05 

-1.685571  E-04 

21.7 

-7.101 162E-04 

-2.997450E-04 

-7.101 162E-04 

-2.997450E-04 

24.8 

2.578272E-05 

3.199933E-04 

2.578272E-05 

3.199933E-04 

27.9 

3.287343E-04 

-5.274567E-05 

3.287343E-04 

-5.274567E-05 

31 

2.75961 9E-04 

3.6401 88E-04 

2.75961 9E-04 

3.6401 88E-04 

34.1 

-3.68931 3E-06 

6.981 51 6E-06 

-3.689313E-06 

6.981 51 6E-06 

37.2 

-2.695758E-04 

-5.676645E-04 

-2.695758E-04 

-5.676645E-04 

40.3 

1.927261  E-04 

-1.785782E-04 

1.927261  E-04 

-1.785782E-04 

43.4 

-3.423437E-04 

8.77041 4E-05 

-3.423437E-04 

8.770414E-05 

46.5 

-6.500707E-05 

5.725958E-04 

-6.500707E-05 

5.725958E-04 

49.6 

8.186082E-05 

2.823280E-05 

8.186082E-05 

2.823280E-05 

52.7 

-1.378475E-04 

-9.820521  E-05 

-1.378475E-04 

-9.820521  E-05 

55.8 

7.90561 7E-04 

2.8521 05E-05 

7.90561 7E-04 

2.8521 05E-05 

58.9 

7.423926E-05 

-2.768488E-04 

7.423926E-05 

-2.768488E-04 

62 

-6.218280E-04 

2.453624E-05 

-6.218280E-04 

2.453624E-05 

65.1 

8.851 132E-05 

-4.144703E-04 

8.851 132E-05 

-4.144703E-04 

68.2 

-1.450531  E-04 

-2.095526E-05 

-1.450531  E-04 

-2.095526E-05 

71.3 

-3  129984E-04 

4.758295E-04 

3.129984E-04 

4.758295E-04 

74.4 

1.079975E-04 

2.570571  E-04 

1.079975E-04 

2.570571  E-04 

77.5 

-2.536621  E-05 

1. 39091 8E-04 

-2.536621  E-05 

1. 39091 8E-04 

80.6 

2.644454E-04 

-2.3841 32E-04 

2.644454E-04 

-2.3841 32E-04 
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TABLE  31-1.  COMPARISON  OF  EXACT  AND  APPROXIMATE  SOLUTION  OF  THE  GREEN’S 

FUNCTION,  CONTINUED 

PCSWAT  EXACT 


Depth,  M 

Real 

Imaginary 

Real 

Imaginary 

83.7 

5.99901 5E-04 

4.540631  E-04 

5.99901 5E-04 

4.540631  E-04 

86.8 

-1.663070E-04 

-2.418384E-04 

-1.663070E-04 

-2.418384E-04 

89.9 

-3.274314E-05 

-7.096214E-04 

-3.274314E-05 

-7.096214E-04 

93 

3.0781 89E-04 

2.401429E-04 

3.0781 89E-04 

2.401 429E-04 

96.1 

-2.750968E-04 

-4.493027E-04 

-2.750968E-04 

-4.493027E-04 

99.2 

-3.570023E-04 

-2.424585E-05 

-3.570023E-04 

-2.424585E-05 

102.3 

-3.382565E-04 

5.974809E-05 

-3.382565E-04 

5.974809E-05 

105.4 

2.291 967E-05 

-2.152185E-04 

2.291 967E-05 

-2.152185E-04 

108.5 

-1 .661 1 76E-04 

5.873839E-04 

-1.661176E-04 

5.873839E-04 

111.6 

-5.8346 19E-04 

4.81 5931  E-05 

-5.834619E-04 

4.81 593 IE-05 

114.7 

3.460832E-04 

2.040280E-04 

3.460832E-04 

2.040280E-04 

117.8 

1.763182E-04 

2.856471  E-04 

1.763182E-04 

2.856471  E-04 

120.9 

-3.589055E-04 

2.928459E-04 

-3.589055E-04 

2.928459E-04 

124 

3.807956E-04 

3.962045E-04 

3.807956E-04 

3.962045E-04 

127.1 

3.738250E-04 

-3.303400E-04 

3.738250E-04 

-3.303400E-04 

130.2 

-3.827622E-05 

5.404839E-04 

-3.827622E-05 

5.404839E-04 

133.3 

3.019754E-04 

2.317492E-04 

3.019754E-04 

2.317492E-04 

136.4 

3.601 126E-04 

-4.246805E-04 

3.601 126E-04 

-4.246805E-04 

139.5 

1.806243E-04 

3.5961 55E-04 

1.806243E-04 

3.5961 55E-04 

142.6 

3.077695E-04 

3.466247E-05 

3.077695E-04 

3.466247E-05 

145.7 

1. 474091 E-04 

5.509475E-07 

1.474091 E-04 

5.509475E-07 

148.8 

4.343452E-04 

-2.044144E-04 

4.343452E-04 

-2.0441 44E-04 

151.9 

2.670685E-04 

9.472733E-05 

2.670685E-04 

9.472733E-05 

155 

-6.337753E-05 

3.846362E-04 

-6.337753E-05 

3.846362E-04 

158.1 

7.366325E-04 

-5.7941 80E-04 

7.366325E-04 

-5.7941 80E-04 

161.2 

5.727772E-05 

2.611 41 6E-04 

5.727772E-05 

2.611 41 6E-04 

164.3 

-1.872883E-05 

4.375409E-04 

-1.872883E-05 

4.375409E-04 

167.4 

7.384581  E-04 

-3.635707E-04 

7.384581 E-04 

-3.635707E-04 

170.5 

-7.500781  E-05 

2.628381  E-04 

-7.500781  E-05 

2.628381  E-04 

173.6 

1.815121E-04 

3.189365E-04 

1.815121E-04 

3.189365E-04 

176.7 

2.031 955E-04 

2.004866E-04 

2.031 955E-04 

2.004866E-04 

179.8 

9.574840E-05 

1.482104E-04 

9.574840E-05 

1.482104E-04 

182.9 

1.778534E-04 

2.0291 85E-04 

1.778534E-04 

2.0291 85E-04 

186 

-5.309431  E-04 

5.078087E-04 

-5.309431  E-04 

5.078087E-04 

189.1 

2.515085E-04 

7.838585E-05 

2.515085E-04 

7.838585E-05 

192.2 

-2.746774E-05 

6.725184E-05 

-2.746774E-05 

6.7251 84E-05 

195.3 

-8.23651 2E-04 

2.094750E-04 

-8.23651 2E-04 

2.094750E-04 

198.4 

1.61 7321  E-05 

-7.181 784E-05 

1.61 7321  E-05 

-7.181 784E-05 

201.5 

2.697329E-05 

-5.51 9383E-05 

2.697329E-05 

-5.519383E-05 

204.6 

-3.961 550E-04 

-3.287576E-04 

-3.961 550E-04 

-3.287576E-04 

207.7 

-2.016652E-04 

-3.599796E-04 

-2.016652E-04 

-3.599796E-04 

210.8 

4.106582E-04 

5.6581 39E-05 

4.106582E-04 

5.6581 39E-05 
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TABLE  31-1.  COMPARISON  OF  EXACT  AND  APPROXIMATE  SOLUTION  OF  THE  GREEN’S 

FUNCTION,  CONTINUED 


PCSWAT  EXACT 


Depth,  M 

Real 

Imaginary 

Real 

Imaginary 

213.9 

1.827767E-04 

-2.275882E-04 

1.827767E-04 

-2.275882E-04 

217 

9.763994E-05 

-3.040782E-04 

9.763994E-05 

-3.040782E-04 

220.1 

3.429965E-04 

2.782708E-04 

3.429965E-04 

2.782708E-04 

223.2 

-8  857475E-05 

4.557982E-04 

-8.857475E-05 

4.557982E-04 

226.3 

3.194940E-04 

1.128602E-04 

3.194940E-04 

1.128602E-04 

229.4 

-2.8461 70E-04 

-1.218884E-04 

-2.8461 70E-04 

-1 .21 8884E-04 

232.5 

-6.002076E-04 

3.137906E-04 

-6.002076E-04 

3.137906E-04 

235.6 

1. 232101 E-04 

3.605080E-05 

1.232101  E-04 

3.605080E-05 

238.7 

-1. 108531 E-04 

-6.0041 08E-04 

-1.1 08531  E-04 

-6.0041 08E-04 

241.8 

2.773905E-04 

-6.520506E-05 

2.773905E-04 

-6.520506E-05 

244.9 

-1.062994E-04 

-1.729520E-04 

-1.062994E-04 

-1.729520E-04 

248 

1.828837E-04 

1.295178E-04 

1.828837E-04 

1.295178E-04 

251.1 

6.021 326E-04 

6.909265E-04 

6.021 326E-04 

6.909265E-04 

254.2 

-5.457961  E-04 

-3.855391  E-04 

-5.457961  E-04 

-3.855391  E-04 

257.3 

-1.582388E-04 

-1.412799E-05 

-1.582388E-04 

-1.412799E-05 

260.4 

-2.968068E-04 

3.720943E-04 

-2.968068E-04 

3.720943E-04 

263.5 

-1.973936E-05 

-5.479777E-04 

-1.973936E-05 

-5.479777E-04 

266.6 

6.874971  E-04 

-5.917463E-05 

6.874971  E-04 

-5.917463E-05 

269.7 

-1.082978E-04 

1.171910E-05 

-1.082978E-04 

1.171910E-05 

272.8 

-7.321751  E-05 

-1.322922E-06 

-7.321 751  E-05 

-1.322922E-06 

275.9 

-2.318374E-04 

6.503385E-04 

-2.318374E-04 

6.503385E-04 

279 

5.030232E-05 

-4.674538E-05 

5.030232E-05 

-4.674538E-05 

282.1 

4.462274E-05 

-5.669501  E-04 

4.462274E-05 

-5.669501  E-04 

285.2 

-1.375995E-04 

-1.844840E-04 

-1.375995E-04 

-1.844840E-04 

288.3 

3.717545E-04 

9.333937E-05 

3.717545E-04 

9.333937E-05 

291.4 

-3.635994E-04 

4.157779E-04 

-3.635994E-04 

4.157779E-04 

294.5 

2.495043E-04 

2.9921 73E-04 

2.495043E-04 

2.9921 73E-04 

297.6 

9.2351 50E-05 

-3.025371  E-04 

9.2351 50E-05 

-3.025371  E-04 

Figure  31.1  presents  a  comparison  of  the  coherent  transmission  loss  between  the  high 
frequency  and  low  frequency  models  in  PC  SWAT  7.  The  waveguide  is  100  m  deep  with  a  rigid 
bottom  and  pressure  release  surface.  The  source  and  receiver  are  both  located  at  a  depth  of  50  m. 
The  frequency  of  the  projector  is  1000  Hz.  Table  31-2  contains  the  sound  velocity  profile  for  the 
waveguide.  In  the  case  of  the  high  frequency  model  a  user-defined  surface  and  reflection  loss 
of  0  decibels  is  used. 

TABLE  31-2.  SOUND  VELOCITY  PROFILE  OF  TEST  CASE 


Depth  (m) 

Sound  Speed  (m/s) 

0 

1500 

50 

1460 

100 

1500 
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FIGURE  31.1  COMPARISON  OF  TRANSMISSION  LOSS  FOR  HIGH  FREQUENCY  AND  LOW 

FREQUENCY  PROPAGATION  MODELS 

Note,  the  incoherent  transmission  loss  for  the  high  frequency  and  low  frequency 
propagation  models  are  not  directly  comparable,  since  each  ray  trajectory  in  the  high  frequency 
model  represents  a  coherent  sum  of  normal  modes. 
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